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ABSTRACT 

We present the results from a detailed analysis of the 0.6 - 10 keV spectra of 
23 ASCA observations of 18 objects. We find that in most cases the underlying 
continuum can be well- represented by a powerlaw with a photon index F ~ 2. 
However we find strong evidence for photoionized gas in the line-of-sight to 
13/18 objects. We present detailed modelling of this gas based upon the ION 
photoionization code. Other studies have been made of the 'warm absorber' 
phenomenon but this paper contains the first consideration of the importance of 
the covering-fraction of the ionized gas and a direct comparison between models 
of attenuation by ionized versus neutral material. 

We find the X-ray ionization parameter for the ionized material is strongly 
peaked at Ux ~ 0.1. The column densities of ionized material are typically 
in the range Nh,z ~ 10^^-10^^ cm~^, although highly ionized (and hence 
psuedo-transparent) column densities up to lO^'^cm"^ cannot be excluded in 
some cases. We also investigate the importance of the emission-spectrum from 
the ionized gas, finding that it significantly improves the fit to many sources 
with an intensity consistent with material subtending a large solid angle at the 
central source. Allowing a fraction of the continuum to be observed without 
attenuation also improves the fit to many sources, and is definitely required 
in the case of NGC 4151. A deficit of counts is observed at ~ 1 keV in the 
sources exhibiting the strongest absorption features. We suggest this is likely 
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to be the signature of a second zone of (more highly) ionized gas, which might 
have been seen previously in the deep Fc A'-shcll edges observed in some Ginga 
observations. We find evidence that the ionized material in NGC 3227 and 
MCG-6-30-15 contains embedded dust, whilst there is no such evidence in the 
other sources 

We discuss these results in the context of previous studies and briefly explore 
the implications in other wavebands. 

Subject headings: galaxies: active - galaxies:nuclei - galaxies rSeyfert - 
X-rays: galaxies 
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1. INTRODUCTION 

Prior to the launch of ASCA, the paradigm for Seyfert 1 galaxies, was that the 2-10 keV 
regime can be well represented by a powerlaw continuum of (photon) index F ~ 1.9 (e.g. 
Nandra & Pounds 1994). Spectra obtained using the Large Area Counter (LAC) onboard 
Ginga also revealed emission line features superimposed on this underlying 'primary' 
continuum attributable to Fe X-shell fluorescence in the 6.4-7.1 keV band. A flattening of 
the spectra was observed above ~ 10 keV thought to be the 'Compton- reflection hump'. 
Evidence for a X-shell absorption edge in the 7.1-8.9 keV band due to highly-ionized Fe 
was also reported in a number of cases (Nandra & Pounds 1994). The combination of these 
features offers an explanation as to why simple powerlaw models to the data obtained by 
earlier instruments typically revealed flatter (F ~ 1.7) spectra (e.g. Turner & Pounds 1989). 
There are a few exceptions to this rule, most notably NGC 4151 which exhibits a continuum 
noticeably flatter than the average (F ~ 1.5) and absorption within a substantial column of 
gas (~ 10^^ cm~^) along the line of sight covering most (but possibly not all) of the central 
continuum (e.g. Yaqoob, Warwick, Pounds 1993 and references therein). 

This picture has been generally supported by ASCA observations, with the addition 
that the higher spectral resolution afforded by the onboard instruments have shown the Fe 
Ka line to be signiflcantly broadened (Mushotzky et al. 1995; Tanaka et al 1995). However, 
the situation in the soft (< 1 keV) X-ray regime is somewhat less clear. Historically, a 
confusing variety of spectral and temporal properties have been seen and/or implied in 
Seyfert 1 galaxies in the soft X-ray regime. A number of sources seemed to show a steep 
and rapidly variable 'soft excess' which could be characterised by either a low-temperature 
thermal component (e.g. Arnaud et al. 1985) or a steep second powerlaw (e.g. Turner 
& Pounds 1989; Turner, George & Mushotzky 1993). Data obtained using the solid-state 
spectrometer (SSS) onboard the Einstein Observatory showed evidence for a complex 
spectral form in the soft X-ray regime (Turner et al 1991), in some cases this was thought 
to be due to leakage of the primary continuum through a patchy absorber (e.g. Reichert et 
al. 1985). 

Evidence for absorption by ionized material along the line-of-sight to active galactic 
nuclei (AGN) was flrst obtained in observations of the QSO MR225 1-178. Using data 
from the Monitor Proportional Counter (MPC) onboard the Einstein Observatory, Halpern 
(1984) reported the column density in this source to have increased by a factor ^ 4 between 

two observations separated by ~ 1 yr. The excess flux observed in a subsequent observation 
using the High Resolution Imager (HRI) on the Einstein Observatory led Halpern to suggest 
an explanation in terms of material with a column density ~ 10^^ cm~^, but photoionized 
such as to be transparent below the i^-edge of Ovil at 740 eV. Such an interpretation 
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was strongly supported by EXO SAT observations of MR2251-178. Pan, Stewart & Pounds 
(1990) found an inverse correlation between the inferred column density and the source flux, 
and demonstrated that this was consistent with the behaviour expected from photoionized 
material along the line-of-sight. Similar behaviour was observed in NGC 4151 (Fiore, Perola 
& Ramano 1990; Yaqoob & Warwick 1991) and MCG-6-30-15 (Nandra, Pounds & Stewart 
1990), establishing the effect for low and high luminosity sources. Supporting evidence for 
such "warm absorbers" was provided by Ginga observations of emission-line AGN which 
revealed iron i^-shell features in some sources, suggestive of an origin in ionized material 
with columns densities ^ 10^^ cm~^. Nandra & Pounds (1994) find this to be a common 
occurrence, with 12 out of 27 sources in their sample of Seyfert 1 galaxies and Narrow 
Emission Line Galaxies (NELGs) showing such a component. However, the low spectral 
resolution of Ginga meant that the results were likely to be sensitive to the assumptions 
concerning the adjacent Fe i^'-shell emission line. 

Confirmation of the existence of warm absorbers came from ROSAT Position Sensitive 
Proportional Counter (PSPC) observations, which revealed absorption edges attributable 
to ionized oxygen in numerous sources (e.g. MCG-6-30-15, Nandra & Pounds 1992; 
3C351, Fiore et al. 1993; NGC 3783, Turner et al. 1993a). Early ASGA observations 
(e.g MCG-6-30-15, Fabian et al. 1994, Reynolds et al. 1995; NGC 3783, George, Turner 
& Netzer 1995) immediately confirmed the PSPC results and showed evidence for warm 
absorbers in many other sources (e.g. NGC 3227, Ptak et al 1994; NGC 4051, Mihara et al 
1994; Guainazzi et al 1996; EXO 055620-3820.2, Turner, Netzer & George 1996; Mrk 290, 
Turner et al 1996; PG1114-h445, George et al. 1997a). Furthermore Reynolds (1997) finds 
evidence for such a component in at least half of a sample of ASCA observations of 24 
type-1 AGN. The presence of highly-ionized, absorbing gas instrinsic to Seyfert 1 nuclei is 
further confirmed by the presence of 'associated absorbers' seen in the UV in many sources. 
Discussion of these observations and the attempts to relate the UV and X-ray absorbers 
made to-date is postponed to §8.9.2 . 

Comparison of ASGA data with the predictions of detailed photoionization calculations 
have revealed spectroscopic evidence for a second ionized absorber in NGC 3516 (Kriss et 
al. 1996a). In addition the absorption edges due to ionized oxygen appear to vary in a 
way which is inconsistent with the behaviour of a single absorber in at least some sources. 
In MCG-6-30-15, the depth of the Ovill edge changes on a timescale ~ 10^ s with no 
corresponding change in the Ovil edge (Reynolds et al 1995; Otani et al. 1996), whilst the 
converse is true in NGC 4051 (Guainazzi et al. 1996). 

The presence of emission lines at energies ^ 1 keV has been suggested in a number 
of previous observations (e.g. Turner et al 1991). A blend of emission lines due to Ovil 
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(0.568-0.574 keV) has been claimed in the ASCA data from NGC 3783 (George, Turner 
& Netzer 1995) and potentially seen in other sources (e.g. IC 4329A, Cappi et al. 1996). 
However the reality of such features has since been disputed, because of the realization of 
specific inadequacies in the low energy calibration of the ASCA instrumentation (see §0). 



Here we present the third in a series of papers describing the X-ray properties of a 
sample of 18 Seyfert 1 galaxies, using data obtained by ASCA. The sample is comprised 
of 23 ASCA observations of Seyfert class 1.0-1.5 galaxies (as defined by Whittle 1992) 
performed prior to 1994 May 01. In Nandra et al. (1997a, hereafter Paper I), we presented 
imaging and timing data for this sample. In that paper we demonstrate that all of the 
sources show evidence for variability on timescales of minutes-hours, with an amplitude 
anti-correlated with the source luminosity, confirming previous results. We also showed that 
for a number of sources the variability amplitude is greater in soft X-rays (^ 2 keV) than 
in hard X-rays (for 10 of the 18 sources on timescales of ~ 1 hr). In Nandra et al. (1997b, 
hereafter Paper II), we presented the 3-10 keV spectra of the sample sources, and found 14 
out of the 18 sources to contain an iron Ka line which is resolved in the ASCA SIS, with 
mean width a ~ 0.4 ±0.1 keV for a Gaussian profile. However, we found many of the line 
profiles to be asymmetric, suggesting an origin from the innermost regions of the putative 
accretion disk. 

In this paper, we present the 0.6-10.0 keV spectra of the sample sources. The 
extension of the spectral analysis down to 0.6 keV requires consideration of the properties 
of the circumnuclear absorbing material as well as any steepening of the spectrum at low 
energies. We present detailed modelling of absorption by photoionized material and discuss 
likely emission contributions as well as combinations of several emission and absorption 
components. We present the range of properties applicable to our sample, and discuss the 
characteristics and statistics of the ionized absorbers. We also consider the the importance 
of allowing a fraction of the underlying continuum to be unattenuated by ionized gas, and 
the inclusion of the emission features expected from the ionized material. As will become 
apparent below, there is an obvious variation in signal-to-noise ratio between different 
datasets and hence there are variations in detectability of the ionized gas. As in the two 
previous papers our intention is to study the mean properties of the sample and we only 
include detailed notes on individual sources as an Appendix. 

We describe the observations and data reduction in §|^, and the photoionization models 
in §0. Our spectral fitting procedure is introduced in and the results presented in §^ and 
§1^. In §1^ we discuss the objects for which there in more than one observation within our 
sample. Our results are discussed within the general context of the physical properties AGN 
in §1^. We sumarize our overall findings in §^ along with outlining a number of important 
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issues that remain to be resolved. 



2. OBSERVATIONS AND DATA REDUCTION 

The ASCA satelhte consists of four identical, co-aligned X-ray telescopes (XRTs; 
Serlemitsos et al. 1995). Two solid-state imaging spectrometers (known as SISO and SISl), 
each consisting of four CCD chips, sit at the focus of two of the XRTs, and provide coverage 
over the ~0.4-10 keV band (Burke et al. 1994). Two gas imaging spectrometers (known 
as GIS2 and GIS3) sit at the focus of the other two XRTs, and provide coverage over the 
~0.8-10 keV band (Ohashi et al. 1996 and references therein). Further details on the 
satellite, its instrumentation and performance can be found in Tanaka, Inoue & Holt (1994), 
Makishima et al. (1996), and in the ASCA Data Reduction Guid^ and references therein. 

2.1. The Sample 

Our sample is listed in Table 1 along with the redshift (z) of the source, the axial 
ratio [a/b) of the host galaxy, the effective Galactic hydrogen column density along that 
line-of-sight (Nj^q) as derived from 21 cm measurements, a number of parameters indicative 
of the multiwaveband spectra of the objects, and X-ray luminosity in two bands. The mean 
flux at 220 nm (< /220 >) in the rest-frame of the source, and the mean ratio of the flux at 
125 nm (rest-frame) to /220 for observations carried out within 48 hours of each other from 
the /f/i? observations were calculated from Courvoisier & Paltani (1992). The errors quoted 
for {fi25/ f 220) obs and < /220 > are la or 10%, whichever is the larger, and datasets flagged 
'dubious' by Courvoisier & Paltani were excluded. The observed flux at 2 keV (/2kev) in 
the rest-frame of the source is from the best-fitting model presented in (in all cases the 
statistical error is larger than the differences in /2keV between different models). In the 
case of objects for which we have more than one dataset, /2keV is the exposure-weighted 
mean. The radio-to-optical energy spectral index, aro, is calculated using < /220 > and the 
observed flux at 5 GHz from Veron-Cetty & Veron (1993). The optical-to-X-ray energy 
spectral index, Uox, is calculated using < /220 > and /2kev- It should be stressed that 
these values of aro and aox are based non-simultaneous observations and hence should be 
regarded with appropriate caution. The errors associated aro and aox are rather difficult 
to determine without estimates of the typical amplitude of variability (on timescales of 



^formerly known as the ABC Guide to ASCA Data Analysis and available via 
http : / /heasar c . gsf c . nasa . gov/ 
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months-years) in each band. Unfortunately, such amphtudes are not well determined for 
most of the objects considered here. However even assuming a factor 5 variability in each 
of the three bands (which is likely to be an overestimate in the radio and optical bands), 
for the typical values of aro ~ 0.1 and aox ~ 1-1 derived for these objects, one obtains 
relatively small errors of ~ ±0.1 and ~ ±0.2 respectively. Finally, Table 1 gives the derived 
luminositie^ of the underlying continuum (i.e. correcting for all line-of-sight absorption) 
in the (rest-frame) 0.1-10 keV (Lx) and 2-10 keV (i^2-io) bands. As in the case of /2keV5 
these are calculated from the best-fitting model presented in §H, and are exposure-weighted 
means in the case of objects with multiple datasets. All fluxes and luminosities quoted for 
the sources will be derived using the SISO detector. 



2.2. Data Selection 

The selection criteria and data analysis methods are presented in Paper I. In summary, 
we have used data from the US public archive at NASA/GSFC. Standard criteria were 
applied to reject poor quality data and source spectra extracted using the FTOOLS/XSELECT 
package for each instrument. Background spectra were extracted from source- free regions 
of the instrument. 

In brief, we consider only the gain-corrected (Pulse Invariant or PI) channels for each 
instrument!]. Furthermore, here we consider only the time-averaged spectrum for each 
instrument during each observation in order to maximize the signal-to-noise ratio. The 
raw spectra were grouped such that each resultant channel had at least 20 counts per 
bin, permitting us to use minimization in the spectral analysis. All such fits were 
undertaken using the XSPEC (v9.01) spectral analysis package (Arnaud 1996). Appropriate 
data redistribution matrices were used (those released 1994 Nov 09 for the SIS, and those 
released on 1995 Mar 06 for the GIS). The effective area appropriate for each dataset 
was calculated using ascaarf^ (v2.5). This does not include a parameterization of the 
azimuthal dependence of the point spread function which can cause inaccuracies in the 
normalization, particularly in the case where the source region is not circular. Uncertainties 
in determining the source centroid can also result in normalization errors, and discrepancies 
remain between the absolute flux- calibration of the four instruments. These factors lead 
to cross-calibration uncertainties between the four instruments, but as the point-spread 

^Hq = 50 km s^^ Mpc^^ and go — 0.5 are used throughout 

^using sispi(vl.O) and the gam-history file released on 1994 Jul 28 in the case of the SIS data 
^where, in the case of the SIS we used the Gaussian 'fudge', but not the filter 'fudge' 
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function of the XRT is not strongly dependent on energy, we do not expect them to 
introduce distortions into the individual spectra. We therefore chose to fit the spectra from 
all four instruments simultaneously, but allowed the normalization of that model to be free 
for each detector. For the datasets considered here, the two SIS detectors typically agree in 
flux to within ~l-3%. The normalization of the GIS detectors are currently forced to agree 
with the average of the SIS detectors for the Crab. 

The observing log is presented in Table 2, where the UTC calendar date is for the 
start of the observation. As in Papers I & II, we refer to each observation using the name 
of the source and, in cases where there are multiple observations, the observation number 
in brackets (using the same numbering scheme^ as in Papers I & II). Table 2 contains the 
exposure time (texp), the count rate {Rate) in the 0.6-10 keV band, and the total number of 
photons in the 0.6-5.0 keV and 7.0-10 keV bands, for SISO (only) after all screening criteria 
have been applied (Nphot)- Also listed in Table 2 are the the total number of spectral bins 
{Npts] for all instruments) used in the analysis, and the additional number of good-quality 
spectral bins in SISO and SISl below 0.6 keV (AiVo.e). The latter two quantities enable 
the reader to judge the goodness-of-fit for the results of the spectral analysis presented 
below. Finally in Table 2 we also provide references to work published previously on each 
observation. 



2.3. Calibration Uncertainties 

We wish to fit several spectral features observed in the 0.4-10 keV range. However, 
there is known to be some uncertainty in the low energy calibration of ASCA that must be 
addressed first. While the origin of this problem is currently unclear (Gendreau, p.comm), 
the manifestation of the problem is in the lowest channels of the SIS data. Specifically, the 
ASCA Guest Observer Facility at NASA/GSFC reports that SIS spectra from a variety of 
objects of different astronomical classes consistently indicate uncertainties in the XRT/SIS 
effective area at energies < 0.6 keV, and especially around 0.5 keV. This has been mainly 
due to the inability to interpret soft X-ray features in the spectra from 3C 273 (the SIS/GIS 
cross- calibration target). However, recent simultaneous observations of 3C 273 with ASCA 
and BeppoSAX may improve this situation, and already indicate that the systematic error 
of the SIS relative to the low-energy concentrator system (LEGS) on BeppoSAX may be 
no more than ~ 2 x 10^° cm~^ (Parmar et al 1997). Further details on the low-energy 



^note that observations which failed to meet our selection criteria are still assigned a number, but not 
included here - see Paper I for details 
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calibration of the SIS are given in Dotani et al. (1996) and from the ASCA GOF at 
NASA/GSFCQ 

Thus, here we have chosen to ignore all data below 0.6 keV during the spectral analysis, 
although we do compare the extrapolation of the various best-fitting models below 0.6 keV 
as a diagnostic of the quality of the fit (see § |4.1|) . For the purpose of this paper, we do 
not use the 5.0-7.0 keV band data in the fitting process, the consistency of those data 
with the fits is discussed later. This approach avoids any requirement to include a complex 
parameterization of the iron K-shell line profile in our consideration of the fits to the 
broad-band continuum shape. Thus the spectral analysis was performed on simultaneous 
fits to the 0.6-5.0, 7.0-10.0 keV data. This provides the best method of determining 
the broad-band X-ray continua. (Note that in all figures presented below which display 
data/model ratios, the corresponding ratios are included for the 0.4-0.6 keV and 5.0-7.0 
keV ranges to illustrate where these data lie with respect to the corresponding model). 



3. THE PHOTOIONIZATION MODELS 

The ionized absorber/emitter models used in this paper were generated using the 
photoionization code ION (Netzer 1993, 1996, version I0N95) to calculate the physical state 
of a slab of gas when illuminated by an ionizing continuum. The code includes all important 
excitation and ionization processes, full temperature and radiative transfer solutions, and 
emission, absorption and reflection by the gas are calculated self-consistently assuming 
thermal and ionization equilibrium. Since the gas temperature is well below the equilibrium 
plasma temperature, only a few transitions affect the cooling, all of which are included. We 
assume an ionizing continuum typical of an AGN with a blue bump. Specifically we assume 
the 'weak IR' case of Netzer (1996), with a photon index Fq = 1.5 in the optical/UV band 
from 1.6-40.8 eV, and photon indices in the range 1.5 < F < 2.5 in the X-ray band from 
0.2-50 keV. The continuum in the XUV band is also assumed to be a power law connecting 
the fluxes at 40.8 eV and 0.2 keV, such that the ratio of the fluxes at 250 nm and 2 keV is 
always /25o//2keV = 8.1 x 10^ (corresponding to an optical-to-X-ray energy index aox = 1-5 
- but see below). Following Netzer (1996), we assume undepleted 'cosmic abundances' 
(with (He, C, N, O, Ne, Mg, Si, S, Fe)/H = (10^, 3.7, 1.1, 8, 1.1, 0.37, 0.35, 0.16, 0.4)/10-^) 
and a constant density of uh = 10^° cm~^ throughout the slab. However it should be noted 
that as discussed in Netzer (1996), for the same illuminating continuum, a wide range of 



for up-to-date information see 
http : / /heasarc . gsf c . nasa . gov/docs/asca/ascagof . html 
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densities (10^ ^ uh ^ 10^° cm~^) give rise to similar equilibrium solutions (and hence 
similarly shaped absorption/emission models). Thus under steady-state conditions our 
models are insensitive to the precise value of uh assumed. The ionization state of the gas is 
parameterized by the ionization parameter, Ux, defined as 



where Q{E) is the rate of photon emission at energy E, r the distance from the source to 
the illuminated gas. In the past, ionization parameters (defined in various ways) involving 
an integration of the entire UV-X-ray continuum above 13.6 eV (with an upper limit often 
set to 13.6 keV) have been commonly used in the literature. As pointed out in Netzer 
(1996), the use of Ux is preferable in the study of the effects of ionized gas in the X-ray 
regime as the level of ionization of the dominant species is then directly proportional to Ux- 
Of course the use of Ux has the additional advantage that it is independent of the form of 
the continuum assumed in the 13.6-100 eV band which is unobservable due to absorption 
in our own galaxy for low redshift objects. For convenience, the relationships between Ux 
and such ionization parameters are shown for a number of spectral forms in Fig. 1. It 
should be stressed whilst the conversions between Ux and these more traditional ionization 
parameters are obviously a strong function of the assumed continuum < 0.2 keV, the values 
derived for Ux here are almost completely insensitive to such assumptions. 

A series of such models were run for various Ux and thicknesses through an illuminated 
slab, the latter parameterized by the column density of the gas, Nh^^. Two grids were 
produced for a number of values of F in the above range, and converted to a form able to 
be used within XSPEC. The first grid consisted of only the absorption (with values between 
zero and unity) of the incident continuum taking place within the illuminated slab. The 
assumed ionizing continuum was multiplied by this grid during the spectral analysis, and is 
referred to as the ionized- absorber below. The axes of this grid are spaced logarithmically 
over the ranges 3 x 10^° < Nh,z < 6.3 x lO^^ cm-^ and -3 < \ogUx < 1-0. Such a model 
would be applicable if the gas only covers the small fraction of the solid angle surrounding 
the central source along the line-of-sight (the small clouds case of Netzer 1993). 

The second grid of spectra consisted of the continua Compton-scattered by the slab 
(but including the effects of absorption during its passage through the slab), along with the 
emission lines and recombination continua produced within the gas. The normalizations. 
Ash, of the spectra in this grid are calculated by integrating the total emission emerging 
from the shell. Since to first order this emission is isotropic for the values of Nh,z and Ux 
considered here, Agh can be compared to Api to give an estimate of the fraction of the 
hemisphere required to be present in order to produce the observed emission, and hence 




(1) 
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the solid angle subtended by the emitting gas at the ionizing source fl (assuming spherical 
symmetry). Thus combining these spectra from the 'ionized-emitter' with the effects of the 
ionized-absorber on the underlying continuum enables us to investigate a range of physically 
plausible conditions within, and location and geometry of, the photoionized gas. 

4. MODEL FITTING 

The models considered in the spectral analysis include an underlying continuum which 
we assume to be a powerlaw of photon index F and normalization Api at 1 keV. In all cases 
the spectrum emerging from the source passes through a full screen of neutral material at 
zero redshift. This screen is parameterized by an effective column density Nh^ of hydrogen, 
which is constrained to be > A'^^'^q, the Galactic Hi column density along that line-of-sight 
(Table 1). 

We have performed the analysis assuming a variety of spectral models of varying levels 
of complexity. For the purposes of clarity, we present the results in the following order: 

A The power-law continuum is absorbed by an additional screen of neutral 
material along the line-of-sight at the redshift of the source (parameterized 
by an effective column density Nh,z of hydrogen). This model therefore has 3 
interesting parameters: F, Nh^^ and N^fi- 

B As for model A, but where the absorbing material at the redshift of the source is 
ionized and parameterized by Nh,z and ionization parameter Ux (see This 
model has 4 interesting parameters: F, Nh,z, Ux and Nha- 

C As for model B, but the emission features expected from the ionized gas (with 
Nh,z and Ux) are included, parameterized by fl, the solid angle subtended by 
the emitting gas as seen at the central continuum source (see This model 
has 5 interesting parameters: F, N^^z, Ux, ^ and Nha- 

For each model, the spectral analysis is performed under two assumptions: 

(i) that the absorbing material at the redshift of the source completely covers the 
source, and that no continuum photons are able to reach the observer via any 
other light paths. These models are labelled A(i), B(i) and C(i) throughout. 

(ii) that a fraction -D/ of the power-law continuum is observed without attenuation 
at the redshift of the source, whilst (1 — Dj) of the continuum is attenuated by 
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material along the line-of-sight at the redshift of the source. Of course this is 
mathematically equivalent to either an absorber within the 'cylinder-of-sight' 
which does not completely cover the central source and/or to a geometry where 
the central source is completely covered by the absorber but in which some 
fraction of the radiation emitted by the central source in other directions is 
scattered back into the line-of-sight without any spectral change. These models 
are labelled A(ii), B(ii) and C(ii) throughout, and have Dj as an additional 
interesting parameter. 

We have applied all the models to all the datasets, and tabulate all the results even 
though the additional parameters in the more complex models may not be required by 
some datasets. This enables us to present the limits on all the parameters with a view to 
unification schemes. The results assuming models A(i)-C(ii) are presented in §|[ A number 
of additional models are discussed in §|^, and notes on individual sources are given in the 
Appendix. The ultimate goal is to classify objects according to the categories defined above 
and try to establish the fraction of spectra that require the ionized gas component. 



4.1. Model Acceptability 

We use a number of criteria to establish whether we consider a model to be consistent 
with a given dataset. During the spectral analysis we use the x^-statistic to determine the 
best-fitting model parameters. Thus this statistic is used in conjunction with the number of 
degrees-of-freedom [dof) to determine the goodness-of-fit over the energy range used during 
the analysis (ze 0.6-5, 7-10 keV). To this end, we use the probability, P(x^ | dof), that 
the statistic should be less than the observed value of under the assumption that the 
model is indeed the true representation of the data. Thus I dof) =0.5 corresponds to 

a reduced-x^ value of unity, values of I dof) ~ 1 indicate that the data are poorly 

represented by the model, and values of I dof) ~ indicate that the model is most 

likely an over-parameterization of the data. For the purposes of this paper, we consider 
the model to be an adequate description of the data (in the 0.6-5, 7-10 keV band) if 
Pix^ I dof) < 0.95 (but flag occasions when P(x^ | dof) < 0.05). 

In addition, however, we also require that the best-fitting model extrapolates into the 
bands excluded from the fit (i.e. < 0.6 keV and in the 5-7 keV band). In the light of the 
uncertainties in the calibration of the instrument at energies < 0.6 keV (§ p.3|) , we quote 
the ratio of the increase in the x^-statistic to the increase in the number of data points 
( AAfo6 '^' and the weighted mean of the data/model residuals (-Ro.e); when the best-fitting 
model is extrapolated below 0.6 keV. The rationale behind these two parameters is that 
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whilst the cahbration is suspect at these energies, it is considered unhkely to be in error by 
^ 20%. The use of and i?o.6 therefore allows us to identify models, which are deemed 
acceptable > 0.6 keV, but in which the extrapolation to energies < 0.6 keV is inconsistent 
with the suspected size of the calibrations uncertainties. For the purposes of this paper, we 
consider the model to extrapolate to energies < 0.6 keV in an acceptable manner if either 
^^"■^ < 2.0 or Ro Q lies in the range 0.8 < i?o.6 < 1-2. A number of the models include 
line emission, and clearly we also require that the model does not predict more Fe i^-shell 
emission than observed in the 5-7 keV band. The few occasions where this is the case are 
flagged in the text and tables. We use the F-statistic to determine whether the inclusion of 
an additional free parameter significantly improves the fit. For the number of dof typically 
in these datasets F ^ 3.8 and ^ 6.6 correspond to improvements at > 95% and > 99% 
confidence respectively. 

To summarize, the formal criteria for acceptibility used throughout are that 
P{X^ I dof) < 0.95 and either < 2.0 or lies in the range 0.8 < < 1.2. 

However, with the menagerie of objects under investigation, these formal criteria for 
acceptability were found to be potentially misleading in a number of cases. In such cases 
somewhat more subjective criteria were used to determine whether the model provided 
an acceptable representation of the data, but all such cases are noted in the text. For 
convenience, a summary detailing which models are deemed acceptable for each dataset 
along with various other notes is provided in Table 12. 



5. BASIC MODELS 

5.1. Model A(i): Single power law and neutral absorption 

The results of the analysis assuming a single underlying powerlaw absorbed by a neutral 
gas at the redshift of the source, and a neutral gas at zero redshift are listed in Table 3. 
As well as the best-fitting parameters. Table 3 also lists the total x^-statistic for each fit 



iXAi)} the parameters described in § ^?T| used to determine the acceptability of the 
fits. Only 4 datasets fulfill all the formal criteria for acceptability: 3C 120, NGC 6814(1), 
Mrk 509 and NGC 7469. However inspection of the data/model ratios reveals 4 additional 
datasets (Fairall 9, NGC 4593, Mrk 841(2) and MCG-2-58-22) which appear to be well 
represented by the simple powerlaw model. In the case of NGC 4593 and Mrk 841(2), 
model A(i) provides an adequate description of the data over the 0.6-5, 7-10 keV band, but 
the best-fitting solution predicts a slight excess of counts {Rq.q ~ 1.3) when the best-fitting 
model is extrapolated <0.6 keV. In the case of Fairall 9 and MCG-2-58-22 the data/model 
residuals indicate that the best-fitting models fail our criteria for acceptability primarily due 
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to residuals in the 2-4 keV band beyond that expected from statistical considerations alone. 
This region of the ASCA bandpass contains a number of instrumental features, primarily 
the Au M-edge complex due to the coating on the XRTs, suggesting these residuals may 
arise as a result of slight miscalibration of the Au M-edges and/or gain of the detectors at 
time of these observations (see also Fig. 13). However, since the residuals appear to have 
little effect on the best-fitting parameters, we consider model A(i) most likely does indeed 
provide an adequate description of these datasets. A similar problem occurs when a number 
of the other models are compared to these and other datasets, and all such occurrences are 
also indicated in Table 12. 

It can be seen that the low signal-to-noise ratio of the NGC 6814(1) dataset results 
in an extremely low value of xja — ^60 corresponding to reduced value of xl = 0.78 
for 205 dof. As reflected in the low value of Pai = 0.01 listed in Table 3, this indicates 
that for the observation of NGC 6814(1) even this simple model can be considered an 
over-parameterization of the available data. A similar problem occurs when a number of 
the other models are compared to other datasets, and all such occurrences are indicated 
in Table 12. However, it should be noted that the fact that a number of the datasets are 
adequately described by a simple powerlaw model in the 0.6-6.0 keV band is not purely 
an artifact of a low signal-to-noise ratio of the observations. Whilst a low signal-to-noise 
ratio may provide an explanation in the cases of NGC 6814(1) & Mrk 841(2), it certainly 
does not provide an explanation in the case of the 3C 120 & Mrk 509 datasets which have 
the highest signal-to-noise ratios in our sample (as illustrated by these sources having the 
highest values of A'p/io((SISO) in Table 2). Thus, it is already apparent that some Seyfert 1 
galaxies are indeed adequately described by a simple powerlaw, even though the majority 
are not. (We note, however, that there is evidence for spectral curvature in 3C 120 at 
energies ^ 6 keV - see §|6^.) 

Histograms showing the distributions of the photon index (F) and intrinsic column 
density {Nh,z) for these datasets are shown in Fig 2a and Fig 3a respectively. For these and 
all the others presented here, the histograms are generated by assigning a rectangle of equal 
area to each dataset {not each object). The width of the rectangle is determined from the 
68% confidence range for that parameter as listed in the relevant table (and hence its height 
derived). We chose such a method as datasets for which the parameter is ill-determined 
then make wide, flat contributions to the total, whilst datasets for which the parameter is 
determined very accurately make thin, high contributions (which can visually dominate 
the total). We feel the resultant plot is a better representation of the results. In all cases 
the hatched region shows the histogram for the datasets which satisfy all the criteria for 
acceptability (§ |4.1| ). The open histogram shows the effect of also including the datasets 
which just fail to meet our formal criteria but for which the model is considered likely to be 
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applicable (as described in the text and summarized in Table 12). As can been seen from 
Fig 2a, all the datasets for which a simple powerlaw provides an adequate fit have photon 
indices in the range 1.5 ^ F ^ 2. Furthermore all of these datasets have intrinsic column 
densities Nh,z ^ few xlO^° cm~"^ (Fig 3a) and Galactic column densities {Nh,o) consistent 
with those derived from 21 cm measurements (Table 3), although we cannot distinguish 
between Nh,z and N^fi in low-redshift, low-Nfj^z sources in this sample. 



5.2. Model yl(^zzj:Partial-absorption by neutral material 

The results of the analysis assuming model A(ii) are listed in Table 4. This model 
assumes a fraction Df of the observed continuum is a bare powerlaw whilst a fraction 
(1 — D/) is a powerlaw attenuated by neutral gas at the redshift of the source. As stated 
in such a model is applicable for geometries either where there is an absorber along the 
line-of-sight which does not completely cover the central source or where the central source 
is completely covered by the absorber but in which some fraction of the radiation emitted 
by the central source in other directions is scattered back into the line-of-sight. As can be 
seen such a model provides a vast improvement in the fits as measured by the reduction 
in the x^-statistic for the majority of the datasets. Indeed, from the F-statistic (listed as 
F{^) in Table 4) it can be seen that the addition of Df as a free parameter leads to an 
improvement over the results obtained assuming model A(i) significant at >95% confidence 
for all but 1 dataset (NGC 6814(1)), and at >99% confidence for all but that dataset and 
Mrk 841(1,2). 

A total of 10 datasets satisfy all our formal criteria for acceptability assuming model 
A(ii). These include the 4 datasets for which our criteria for acceptability were also satisfied 
by model A(i), along with the datasets of Mrk 335, NGC 3227, Mrk 766, NGC 4593, 
Mrk 841(2) and MCG-2-58-22. The histogram of the best-fitting photon indices from the 
datasets for which model A(u) provides an adequate description is shown in Fig 2b. It 
can be seen that the distribution is extended over the range 1.5 ^ F ^ 2.5, but as for 
model A(i) contains clear evidence for a peak at F ~ 2. The best-fitting solutions to 
all the datasets have Galactic column densities {Nh,o) consistent with that derived from 
21 cm measurements (Table 3) with the exception of 3C 120 {ANh^q ~ 6 x 10^° cm~^) and 
NGC 7469(2) (AA^^.o - 3 x lO^o cm'^). 

The derived model spectra and data/model ratios for the 8 datasets for which model 
A(ii) is acceptable and offers an improvement at >99% confidence over model A(i) are 
shown in Fig 4. For 2 of the remaining datasets (Fairall 9 and IC 4329A) we find model A(ii) 
to offer a significant improvement over model A(i), yet no significant further improvement 
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in the goodness-of-fit is obtained assuming subsequent models presented in §|^. These 
two datasets are also shown in Fig 4, despite our formal criteria for acceptibility not 
being satisfied. As can be seen from Table 4, for many of the datasets (Fairall 9, 3C 120, 
Mrk 841(2), NGC 6814(2), NGC 7469(2) & MCG-2-58-22) the best fitting solution occupies 
the region of high- NH,zM^g^-Df parameter-space. This is also clearly apparent from the 
histograms of Nh,z and Df for these datasets are shown in Figs 3b and 5a (respectively), 
and gives rise to a relatively subtle upturn in the observed spectrum ^ 7 keV in these 
sources (Fig 4). 

It is important to appreciate that given the ASCA bandpass (and typical signal-to-noise 
ratios) we are only sensitive to and/or can constrain N^^z and Df over a limited range 
of parameter space. Clearly the model is insensitive to Nh,z in all models where ^ 1 
(models with Df = 1 being identical to a powerlaw). The lower limit on the SIS bandpass 
(along with calibration uncertainties) leads to us being insensitive to Nh^z ^ 2 x 10^° cm~^, 
and the upper limit on the ASCA bandpass leads to us being insensitive to Nh,z ^ 10^"^ cm~^. 
Two datasets are consistent with Df = I at 68% confidence: Mrk 841(2) & NGC 6814(1). 
These sources are consistent with a single powerlaw, and thus Nh^z is effectively undefined 
(as illustrated by the errors 'pegging' at the limits in Table 4). In addition it should be 
noted that Fairall 9 is consistent with Df = 1 at 90% confidence and statistically acceptable 
solutions to the Mrk 509 dataset obtained at 90% confidence at both high and low values of 
Nh^z- Thus model A(u) does not provide a unique explanation to the datasets which exhibit 
a subtle upturn at the highest energies (i.e. those with a high- N^^zMs^-Df solution). 
Indeed, this might be an indication of a real curvature in the underlying continuum (see 
§ |6.4| ), although it might also arise as a result of spectral variability during the observation 
or (less likely) the result of incomplete background-subtraction. However, it is worth noting 
that the best-fitting models for 4 datasets (Mrk 335, NGC 3227, Mrk 766 & NGC 4593) 
are in regions of parameter-space where the ASCA data are sensitive to the model. These 
sources have 10^^ ^ Nh,z ^fewxlO^^ cm~^ and 0.3 ^ Df ^ 0.7. Nevertheless yet superior 
fits are obtained with subsequent models in all 4 cases (Mrk 335 - § |6.4.2| ; NGC 3227, 
Mrk 766 & NGC 4593 - 



5.3. Model B(i): Complete covering by ionized material 

The results of the analysis assuming the powerlaw continuum passes through ionized 
material, completely covering the source are listed in Table 5. It can be seen that for the 
majority of the datasets such a model provides a vast improvement over model A(i) in 
which the material responsible for the absorption was assumed to be neutral. Indeed, from 
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the F-statistic (F(ff) in Table 5) it can be seen that the addition of Ux as a free parameter 
leads to an improvement over the results obtained assuming model A(i) (§ |5.1| ) significant at 
>95% confidence for all but Fairall 9 and NGC 6814(1), and at >99% confidence for all but 
Fairall 9, 3C 120, NGC 6814(1) and Mrk 841(2). A direct comparison of Xm in Table 5 with 
Xau in Table 4 (both models having the same number of dof) indicates model B(i) also 
offers a superior description of the data than model A(n) for the majority of the datasets. 
The notable exceptions to this are the observations of 3C 120, IC 4329A and (especially) all 
3 observations of NGC 4151. It should also be noted that -^^^ and i?o.6 (Table 5) are both 
greatly improved compared to those derived assuming model A(ii) (Table 4) for a number 
of datasets, indicating model B(i) provides a far superior description of the data <0.6 keV. 
This is true even for a number of datasets that do not satisfy our criteria for acceptability, 
and in most cases the result of a 'recovery' of the absorbed spectrum at low energies (below 
the oxygen iT-shell edge) due to carbon and nitrogen being fully ionized and hence the gas 
pseudo-transparent at these energies. 

A total of 14 datasets satisfy all our criteria for acceptability assuming model B(i). 
These include the 10 datasets for which our criteria were also satisfied by model A(u), 
along with NGC 3783(1,2), NGC 5548 and Mrk 841(1). The derived model spectra and 
data/model ratios for the 8 datasets for which model B(i) is acceptable and offers an 
improvement at >99% confidence over both models A(i) and A(ii) are shown in Fig 6. We 
note that Psi = 0.03 in the case of Mrk 766, indicating this model (and all subsequent 
models) could be considered an over-parameterization of the data. For 2 of the remaining 
datasets (NGC 3516 and MCG-6-30-15(l)) we find model B(i) to offer a significant 
improvement over both models A(i) and A(ii), yet no significant further improvement in the 
goodness-of-fit obtained assuming subsequent models presented in §||. These two datasets 
are also shown in Fig 6, despite our formal criteria for acceptibility not being satisfied. 

The superior description of the data at low energies afforded by model B(i) compared 
to that of model A(ii) is illustrated by comparison of panels for NGC 3227, Mrk 766, 
NGC 4593 and Mrk 509 in Figs 4 and 6, and attributable to the reduction in opacity of the 
ionized gas at low energies. Inspection of the data/model residuals of the remaining datasets 
reveal that both observations of MCG-6-30-15 fail to satisfy our criteria for acceptability 
primarily due to a deficit of counts around ~1 keV. As noted in Table 12, such a deficit 
is also evident in the data/model residuals for NCG 3516 (along with an apparent excess 
scatter across the entire energy range), and in NGC 3783(1,2). The feature may be the 
result of 'excess' absorption by Ovill, and is further discussed in §|^. 



The distribution of the photon index from these fits, shown in Fig 2c, covers a similar 
range and remains peaked at F ~ 2 as for models A(i) and A(ii). The distribution of the 
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column density of the ionized material shows no clear preference, with a broad variety 
of values over the range 10^^ <^ Nh^z ^ 10^^ cm~^. The distribution for the ionization 
parameter (Ux), plotted in Fig 7a, is extended over the range 10~^ <^ Ux <^ 10, but exhibits 
a clear peak at Ux ~ 0.1. It is worth noting that the datasets with ionization parameters 
in the range 1 <^ Ux ^ 10 lie in regions of Nh,z-Ux parameter-space where the ionized gas 
is largely transparent in the ASCA band. Thus the best-fitting spectra for these datasets 
(e.g. Fairall 9, 3C 120, NGC 7469(2), MCG-2-58-22) are essentially powerlaws, consistent 
with the findings above. In Fig 8a we plot the derived values for Nh,z and Ux for the 11 
datasets for which model B(i) offers an adequate description of the data and leads to an 
improvement at > 99% over model A(i). The corresponding values for NGC 3516 and both 
datasets of MCG-6-30-15 are also shown in Fig 8a for comparison. 

Finally we note that all the datasets for which model B(i) provides an acceptable 
solution have best-fitting Nha consistent with that derived from 21 cm measurements 
(Table 3) with the exception of 3783(1) (AiV^.o ^ 2 x 10^° cm'^) along with 3C 120 and 



NGC 7469(2), both of which have excess column densities similar to those found in §p]2 



5.4. Model B(ii): Partial-covering by an ionized absorber 

The results of the analysis assuming model B(ii) are listed in Table 6. The relationship 
of this model to model B(i) is the same as the relationship between models A(i) and A(ii): 
identical except that now a fraction -D/ of the observed continuum is a bare powerlaw whilst 
a fraction (1 — Df) is attenuated by line-of-sight material at the redshift of the source, 
(see §^). It can be seen that model B(ii) provides a further substantial reduction in the 
X^-statistic for many of the datasets. From the F-statistic {F{^) in Table 6) we find that 
the addition of Ux as a free parameter leads to a improvement over the results obtained 
assuming model A(ii) for 17 (16) datasets at >95% (>99%) confidence[^. 

A total of 15 datasets satisfy all our criteria for acceptability assuming model 
B(ii). These include 13 out of the 14 datasets for which model B(i) was found formally 
acceptable^, along with NGC 4151(2,4). As was the case for model B(i), inspection of 
the data/model residuals reveals that a deficit of counts at ~lkeV is the primary reason 
why several of the remaining datasets fail to satisfy our criteria for acceptability (namely 
NGC 3516, NGC 4151(5) and MCG-6-30-15(l,2)). 



^^See the note to Table 6 for an explanation of the apparent increase in for the observations of Fairall 
9, 3C 120 & IC 4392A. 

^^The exception being Mrk 841(1), which falls outside our formal limit of Pbh < 0.95 by just APbu — 0.001 
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Interestingly, we find that the inclusion of the additional free parameter Df compared 
to model B(i) improves the fit at >99% confidence for 10 of the 23 datasets (as shown by the 
values of the F-statistic F{^) listed in Table 6). Furthermore, 8 of these datasets satisfy 
our criteria for acceptability, indicating the spectrum observed from a significant fraction 
of the datasets may indeed contain a component of the underlying continuum which does 
not experience attenuation by the ionized gas. The best-fitting solutions for all 5 datasets 
for which model B(ii) is acceptable and offers an improvement at >99% confidence over 
both models A(u) and B(i) have Df > 0, as illustrated by the derived model spectra and 
data/model ratios shown in Fig 9. For one of the remaining datasets (NGC 4151(5)) we 
find model B(ii) to offer a significant improvement over previous models, yet no significant 
further improvement in the goodness-of-fit obtained assuming subsequent models presented 
in This dataset is also shown in Fig 6, despite our formal criteria for acceptibility not 
being satisfied. 

In Fig 8b we plot the derived values for Nh.z and Ux for the 10 datasets for which model 
B(ii) offers an adequate description of the data and leads to an improvement at > 99% over 
model A(%%). The corresponding values for NGC 3516, NGC 4151(5), MCG-6-30-15(l,2) 
and Mrk 841(1) are also shown in Fig 8b for comparison. 

Again the limitations imposed on parameter-space by the ASCA bandpass and the 
signal-to-noise ratios are reflected in a number of the fits. Obviously, the comments made 



in § ^.2| regarding solutions with ~ 1, and in § ^.3| regarding the transparency of the 
ionized gas in the Iow-Nh,z, high-?7x region of parameter space are also applicable to model 
B(ii). These effects prevent useful constraints being placed on Nh,z in a number of cases. 
Furthermore, the gas for models in the high-iVj:^^^,low-[/x region of parameter space is close 
to neutral, giving rise to spectra identical to those in model A(ii). As for model A(ii), the 
best-fitting models for a number of datasets exhibit a subtle upturn at the highest energies. 

As might be expected, the distribution of the photon index (Fig 2d) covers a similar 
(but slightly broader) range as that derived from previous models, and is still peaked at 
r ~ 2. As for model B(i), there appears to be no preferred value of the intrinsic column 
density (Fig 3d), with the derived values (if errors are also considered) covering over three 
decades in Nh,z- The datasets for which Df > are gathered in the range 0.3 ^ Df ^ 1.0, 
with the exception of the 3 observations of NGC 4151 which have 0.03 <^ Df <^ OA (Fig 5b). 
The distribution in ionization parameter is extended over the range 10~^ ^ Ux <^ 1.0, but 
as for model B(i), there is a clear concentration around Ux ~ 0.1 (Fig 7b). 

Finally, we note that the best-fitting solutions to 13 of the 15 datasets considered 
acceptable have Galactic column densities (Nnfi) consistent with that derived from 21 cm 
measurements (Table 3). Again the exceptions are 3C 120 and NGC 7469(2), both of which 



have excess column densities similar to those found in 
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5.5. Model C(i): Complete covering by an ionized absorber plus an ionized 

emitter 

The results of assuming model C(i) are listed in Table 7. This model assumes the 
powerlaw continuum passes through ionized material, completely covering the source (i.e. 
as for model B(i)), but in addition contains the emission and scattered spectrum expected 
from the ionized gas assuming it subtends a solid angle fl at the central source (see 
We see that the such a model provides further improvement for many of the datasets. 
Using the F-statistic {F{%) in Table 7) it can be seen that the addition of f2 free 
parameter leads to an improvement over the results obtained assuming model B(i) for 15 
(11) datasets at >95% (>99%) confidence. A direct comparison of Xci Table 7 with Xbu 
in Table 6 (both models having the same number of dof) indicates model C(i) also offers a 
significant improvement over the results obtained assuming model B(ii) for NGC 3783(1,2), 
NGC 4051 and MCG-6-30-15(2). We note, however, that model C(i) provides an inferior fits 
compared to model B(u) for the observation of 3C 120 and all 3 observations of NGC 4151. 
This indicates that consideration of the covering fraction of the ionized absorber is more 
important than its emission in these cases. Model C(i) offers only a modest improvement 
or slightly inferior fit to the remaining datasets. 

A total of 15 datasets satisfy all our criteria for acceptability assuming model C(i). As 
summarized in Table 12, this includes all 14 datasets for which our criteria for acceptability 
were also satisfied by model B(i), along with MCG-6-30-15(2). As was the case for models 
B(i) B(ii), inspection of the data/model residuals reveals that a deficit of counts at 
~lkeV is the primary reason why the datasets of NGC 3516, MCG-6-30-15(l) fail to satisfy 
our criteria for acceptability. The spectrum emitted by the ionized gas includes Fe lines 
in the 6.4-7.1 keV band (rest-frame). Since the bulk of this region is excluded from our 
spectral analysis, it is possible that the best-fitting model will be inconsistent with the 
observations within this band. However for all 15 datasets which satisfying our criteria for 
acceptability, we find the best-fitting solutions lie in a region of Nh,z-Ux-^ parameter-space 
such that the predicted line-strength is less than that observed. The derived model spectra 
and data/model ratios for the 4 datasets for which model C(i) is acceptable and offers an 
improvement at >99% confidence over both models B(i) and B(ii) are shown in Fig 10. 

In Fig 8c we plot the derived values for Nh,z and Ux for the 13 datasets for which 
model C(i) offers an adequate description of the data and leads to an improvement at 
> 99% over model A(i). The corresponding values for NGC 3516 and MCG-6-30-15(l) are 
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also shown in Fig 8c for comparison. 

The inclusion of the emission spectrum does not globally change the distributions 
of the r, Nh,z and Ux for these datasets as shown Figs 2e, 3e, and 7c (respectively). In 
general the parameters cover similar ranges as for models B(i) & B(ii), with the datasets 
exhibiting a range of photon indices (but peaked at F ~ 2), column densities and ionization 
parameter (again peaking at Ux ~ 0.1). As can be seen from Fig 11a, the solid angle of the 
ionized emitter is poorly constrained in most cases, but there is a clear peak at f2 ~ 4-7? (i.e. 
consistent with emission from a complete shell). 

Again, the best-fitting solutions to all the datasets considered acceptable have Galactic 
column densities (A^h,o) consistent with that derived from 21 cm measurements (Table 3), 
with the exceptions of 3C 120 and NGC 7469(2), both of which have excess column densities 



similar to those found in 55.2. 



5.6. Model C(ii): Partial-covering by an ionized absorber plus an ionized 

emitter 

The results of the analysis assuming model C(ii) are listed in Table 8. As described 
in this model is the same as model C(i) except that only a fraction {1 — Df) of the 
observed continuum is attenuated by the ionized gas. (The spectral component due to 
the emission and scattered spectrum expected from the ionized gas is unaffected and 
parameterized by fl, just as in model C(i).) It can be seen that model C(ii) provides a 
significant reduction in the x^-statistic for several datasets. From the F-statistic (F(^) in 
Table 6) we find that the addition of as a free parameter leads to a improvement over 
the results obtained assuming model B(ii) for 14 (8) datasets at >95% (>99%) confidence. 
Alternatively, we find that the addition of as a free parameter leads to model C(ii) 
being preferred over model C(i) for 9 (8) datasets at >95% (>99%) confidence. However, 
we note that model C(ii) offers an improvement at >99% confidence over both models 
B(U) and C(i) for only 4 datasets (NGC 4051, NGC 4151(4), NGC 5548 and Mrk 509). 
The derived model spectra and data/model ratios of these are shown Fig 12. It should be 
noted, however, that in the case of NGC 4051, the best-fitting solution assuming model 
C(ii) fails to meet our formal criteria for acceptability. This is mostly likely due to (at 
least in part) spectral variability within the observation giving rise to spectral curvature 
(see § |6.4|) . In addition, the best-fitting solution to NGC 5548 over-predicts the number of 
counts < 0.6 keV (-Ro.6 ~ 1-3) • Finally it should be noted that the best-fitting solutions for 
NGC 5548 and Mrk 509 lie in the high-A'^^^2,high-[/x,high-iI)j region of parameter space. 
The resultant spectra are therefore essentially powerlaws with subtle upturns at both the 
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very lowest and very highest energies (Fig 12). The meaningfulness of such solutions is 
therefore questionable, as it may simply be indicative of spectral curvature (see §6.4|). 

A total of 16 datasets satisfy all our criteria for acceptability assuming model C(ii). 
These include 14 out of the 15 datasets for which model C(i) was found formally acceptable 
(the exception being NGC 5548, which just fails our criteria for extrapolation <0.6 keV as 
noted above), along with 2 of the 3 observations of NGC 4151. As was the case for previous 
models, inspection of the data/model residuals reveals that a deficit of counts at ~lkeV is 
the primary reason why the datasets from NGC 3516, NGC 4151(5) and MCG-6-30-15(l) 
fail to satisfy our criteria for acceptability. We also note that for 5 of the 16 datasets 
satisfying our criteria for acceptability for this model (Mrk 335, 3C 120, NGC 5548, 
NGC 7469(2) and MCG-2-58-22), the best-fitting solutions over-predict the strength of the 
Fe i^'-shell lines in the 6.4-7.1 keV band. In a number of cases the 90% confidence region in 
Nh,z^Ux-^ parameter-space does include solutions consistent with observed line strengths. 
However, as will be shown in § |6.4| , in all cases the solution obtained assuming model C(ii) 
is most likely an artifact of spectral curvature in the underlying continuum. 

In Fig 8d we plot the derived values for Nh,z and Ux for the 11 datasets for which model 
C(ii) offers an adequate description of the data and leads to an improvement at > 99% over 
model A(ii). The corresponding values for NGC 3516, NGC 4151(5), MCG-6-30-15(l) and 
NGC 5548 are also shown in Fig 8e for comparison. 

As for model C(i), we find the inclusion of the emission spectrum does not globally 
change the distributions of the F, Nh,z, Df and Ux for these datasets as shown in Figs 2d, 
3d, 5c and 7b (respectively). In general the parameters cover similar ranges as for previous 
models. In most cases, the emission from the ionized gas forms a significant fraction of 
the observed spectrum only at low energies, and is dominated by emission lines due to 
i^-shell processes of C, N & O. However, it should be noted that in a number of cases 
(e.g. NGC 4151(4) & Mrk 509) the Fe ii'-shell emission line and scattered continuum do 
contribute to the spectrum observed above ~3 keV. 

Again we find the best-fitting solutions to all the datasets considered acceptable have 
Galactic column densities {Nh,o) consistent with that derived from 21 cm measurements 
(Table 3), with the exceptions of Mrk 509 (with ANh,o ^ 5 x lO^" cm'^), and 3C 120, 
NGC 3783(1) and NGC 7469(2), all of which have excess column densities similar to those 
found previously. 

Finally, in Fig. 13 we show the mean data/model ratios (in the observer's frame) in 
the SIS and CIS for the 16 datasets satisfying our criteria for acceptability assuming model 
C(ii). This plot was constructed by taking the error- weighted average of the individual 
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data/model ratios for each detector from each of the appropriate observations. These 
averages were then rebinned for display purposes. Caution is urged when interpreting such 
plots. They potentially contain artificial features introduced by time-dependent errors in 
the calibration of the instruments (such as slight offsets in the gain of the detectors at the 
time of each observation), along with real features due to errors in the calibration of the 
instrument, as well as real features of astrophysical origin. Furthermore such plots are 
dominated by the datasets with the highest signal-to-noise ratio. Nevertheless, we have 
found Fig. 13 helpful in judging the reality of features/noise in individual datasets. 



6. ADDITIONAL MODELS 

In the previous section we made the premise that the 'primary' X-ray continuum was a 
powerlaw. We then compared the observed spectrum with that expected assuming a number 
of (probably simplistic) assumptions regarding the geometry and physical conditions within 
the circumnuclear material. We found that one or more of this family of models to be 
consistent with 17 of the 23 datasets considered within our formal criteria of acceptability 
(§ [4.1| ). Of the 6 remaining datasets, the primary reason why the criteria of acceptability are 



not achieved is due to excess scatter in the residuals within the 2-4 keV band for Fairall 9 
(see Fig. 4, the observed spectrum otherwise being consistent with model A(ii)), a deficit of 
counts at ~1 keV in the case of NGC 3516, NGC 4151(5) and MCG-6-30-15(l) (see Figs 6 
&; 9), and a more complex spectral form in the case of NGC 4051 and IC 4329A (Figs 12 
& 4 respectively). In this section we extend our spectral analysis to address a number of 
issues raised in the previous section. 



6.1. The redshift of the ionized gas 

In §1 it was assumed that the ionized gas along the line-of-sight was situated close to 
(and not in motion relative to) the active nucleus (i.e. at the redshift of the host galaxy, 
Zgai). The sharp features imprinted by the ionized material (primarily Ovil and/or Ovill 
absorption edges for most datasets, with rest-frame energies of 740 and 871 eV) facilitate 
testing such an hypothesis. The spectral resolution of the SIS is a complex function of 
energy, time and various other instrument effects. However at the epoch of the observations, 
AE/E ~ 0.07-0.09 in the 0.6-1.0 keV band (e.g. Dotani et al 1996). Useful constraints on 
the redshift of the ionized gas Zion are only likely to be obtained for those datasets with the 
strongest absorption and/or emission features. Nevertheless, we have repeated the spectral 
analysis for models B(i)-C(ii) allowing the redshift of the ionized gas (zjon) to be a free 
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parameter. However, in the vast majority of cases we found the 90% confidence range for 
Zion to encompass both Zion = and Zion = Zgai- We find the most stringent constraints on 
Zion to be provided by model C(i) apphed to 3783(1), where Zion = Zgai — 0.002to;oi2 (where 
Zgai = 0.009). The formal limits exclude Zion = and/or Zion = Zgai for a few dataset/model 
combinations, but further investigation reveals that in each case the result is due to fitting 
small statistic fluctuations in the data. Given the recent detection of ionized gas at the 
redshift of the host galaxy in another AGN (the high-redshift quasar PG 1114+445; George 
et al. 1997a) and lack of evidence to the contrary here, we conclude our assumption to be 
valid. 



6.2. Evidence for a second, ionized absorber 

During the analysis presented in §^ a number of datasets were observed to consistently 
exhibit a deficit in their data/model ratio at ~1 keV. This feature was most evident in 
both datasets of NGC 3783 (e.g. Figs 6, 9 & 10) and MCG-6-30-15 (Figs 9 & 10), and 
in the single dataset from NGC 3516 (Fig 6). As noted in the introduction, the ~1 keV 
deficit has also been seen by a number of other workers, and since it appears close to the 
K-sheW edges of O and Ne, it has been suggested that it is the result of absorption by a 
second cloud/screen of ionized gas along the line-of-sight (e.g. Kriss et al. 1996a; Otani et 
al. 1996). 



As will be discussed in § p.8| there are a number of explanations for such a feature, 
and a detailed investigation of these are beyond the scope of the current paper. However, 
in Table 9 we hst the results obtained when an absorption edge of optical depth tonc is 
added to model B(i). The (rest-frame) energy of the edge, EoNe, was allowed to vary 
over the range 0.568-1.362 keV, thus covering the edges due to Ovi-Nex. As can be seen 
from the values of the F-statistic FpNe/Bi, the addition of this component improved the 
goodness-of-fit at > 99% confidence for all 5 datasets. A comparison with Table 5 shows 
that the effect of including such an additional edge is to reduce the best-fitting values of Ux 
and Nh,z- Furthermore, in all but one case (NGC 3783(1)) we find the best-fitting energy 
for the edge to be consistent with that due to Ovill. We also tested the addition of such 
an edge to model C(i) and obtained similar results, except that the edge in NGC 3783(1) is 
now also consistent with that due to Ovill. 

Finally, it should be stressed that such a deficit is not seen in several of the other 
datasets presented here, and hence unlikely to be the result of a calibration error. For 
example, in the case of the high signal-to-noise, powerlaw dataset 3C 120, we find the 
additional of an Ovill edge to model B(i) makes an insignificant improvement in the 
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goodness-of-fit (A^^ < 1), with ToNe ^ 0.02 at 90% confidence (for one interesting 
parameter) . 

We consider these results to be indicative of a second absorbing system in some objects. 
This possibihty, along with the fact that the 1 keV deficit is present in the datasets with 
the strongest absorption features, is discussed in §|8l8. 



6.3. Implications for Fe i^'-shell band 

6.3.1. Fe K -shell emission 

In Paper II we presented an analysis of the spectra of all 23 datasets in the 3-10 keV 
band. The approach used was to approximate the underlying continuum with a powerlaw in 
order to perform a detailed investigation the profile of the Fe Ka line. In the present work, 
we have ignored the 5-7 keV band, which contains the bulk of the Fe Ka line emission. 
Clearly it is inappropriate to repeat here the full analysis presented in Paper II. However 
there are two consistency checks which should be performed. 

First, in §^ we have shown that the 0.6-10 keV continuum in these sources can be well 
represented by a powerlaw imprinted with (primarily absorption) features due to ionized 
gas along the line-of-sight. Comparison of Tables 3-6 with the corresponding results from 
Paper II reveal consistent values of P. Thus the approximation of the local continuum used 
in Paper II is indeed valid, and hence the intensity and profile of the line emission presented 
in Paper II is also valid. This is illustrated in Fig 14 in which we show the mean profiles 
obtained for models B(ii), C(i) and C(ii) compared to that presented in Paper II. 

Second, models C(i) and C(ii) contain the emission spectrum for the ionized gas, which 
includes significant Fe i^-shell features in certain regions of Nh,z^Ux parameter-space. We 
have therefore repeated the analysis for these models, but including the 5-7 keV band 
and an emission line as an additional spectral component. The profile of the emission line 
assumed is that expected from the innermost regions of an accretion disk under the infiuence 
of the strong gravitation effects of the central black hole. For convenience, we have followed 
Paper II and adopted the parameterization of Fabian et al (1989), fixing the inner and outer 
radii of the disk at Ri = 6Rg and Ro = lO^Rg (respectively, where Rg is the gravitational 
radius), the axis of rotation with respect to the line-of-sight of i = 30°, and the (rest-frame) 
energy of the line to be E^a = 6.4 keV. The emissivity of the line is parameterized as a 
powerlaw as a function of radius (oc i?"^'^), leaving only the equivalent width of the line, 
Wxa, as a free parameter. As may be expected from Fig 14, for both models C(i) and C(ii), 
we find results consistent with those presented in Paper II, with mean equivalent widths of 
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the diskline Wxa — 290 ± 30 eV, compared with Wko — 270 ± 14 eV derived from Table 5 
of Paper II. 



6.3.2. Fe K -shell absorption 

Ginga observations indicated a large fraction 40%) of type-1 AGN might possess 
a significant Fe absorption edge in the 7.1-8.9 keV band (Nandra & Pounds 1994). 
However given the spectral resolution of the LAC and the presence of the strong Fe 
K-sheW emission in the 6-7 keV band, the parameters of such a feature could not be 
unambiguously determined by Ginga. Nevertheless, Nandra & Pounds did find a preferred 
energy of ~ 7.9 keV (equivalent to Fexvill) and effective hydrogen column densities 
Nh,z ~fewxlO^^ cm~^. 

Unfortunately the end of the effective bandpass of the SIS and GIS at ~ 10 keV and 
the low signal-to-noise ratio of most datasets in the 7-10 keV band makes the study of 
such a feature extremely difficult with ASGA. We have not attempted such a study here, 
especially since we employed an analysis technique whereby we explicitly excluded the Fe 
iiT-shell emission (in the 5-7 keV band) from the majority of the fits. However as can be 
seen in Figs 6, 9, 10 & 12, our models for the photoionized gas do not predict significant Fe 
i^'-shell absorption edges for the values of Ux and Nh^z typically derived in §^ (the ASGA 
results being primarily based on O i^-shell, Ne i^'-shell and Fe L-shell absorption features). 
If confirmed, the absorption features seen by Ginga imply an additional, large column 
density of highly-ionized gas along the cylinder-of-sight, possibly related to that suggested 
as responsible for the '~ 1 keV deficit' described in §0. However we note that there is no 
obvious evidence for such features in Fig 14. 



6.4. Constraints on a more complex continuum form 

Given the diversity in the observed spectra presented here, the assumption that the 
continuum is well represented by a single powerlaw over the entire 0.4-10 keV band can 
be questioned. Indeed, based on both observational and theoretical arguments, in the 
past a number of more complex spectral forms have been suggested for the underlying 
X-ray continuum in Seyfert 1 galaxies. These can be divided into those in which the 
underlying emission mechanism gives rise to a continuum with a more complex form (a 
simple case case a continuum represented by the sum of two powerlaws), and those in which 
a 'primary' powerlaw continuum is augmented by other sources of emission (both thermal 
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and non-thermal, inside and outside the active nucleus). 

In §^ we found a number of datasets consistently satisfied by best-fitting parameters 
which could imply a gradual flattening of the observed spectrum to higher energies (e.g. 
Mrk 335, NGC 4051, IC 4329A). The data/model residuals of a number of other datasets 
also imply a 'soft excess' of counts at the lowest energies and/or a 'hard tail' at the highest 
energies (e.g. NGC 3516, 3C 120). In the light of this, we have investigated two forms of the 
continuum commonly assumed for Seyfert 1 galaxies: the addition of a 'Compton Reflection' 
component, and an underlying continuum represented by the sum of two powerlaws. We 
have also investigated the inclusion of a thermal emission component at the softest energies. 

It should be remembered, however, that spectral curvature can of course also arise as 
the result of spectral variability within an observation. Throughout this paper we consider 
only the mean spectra. However in Paper I we presented evidence for variability in many of 
the datasets. In most datasets any spectral variability was of relatively low amplitude, but 
a number of cases (primarily the low luminosity sources, and especially NGC 4051), the 
amplitude could be at a level to start affecting our analysis of the mean spectra. 

6.4- i- The effects of including 'Compton Reflection' 

'Compton- reflection' of the 'primary' continuum by optically-thick, neutral material 
out of the line-of-sight is often invoked as an explanation of the high equivalent width of the 
Fe Ka line emission and flattening of the spectrum ^ 10 keV in many AGN (e.g. Nandra & 
Pounds 1994). The presence of an unresolved reflection component can flatten the observed 
continuum, so the underlying continuum could be slightly steeper than observed, which 
could have implications on the best-fitting values of the parameters associated with the 
ionized gas. We have therefore repeated the analysis for model C(ii) including the effects of 
the 'Compton-reflection, which we parameterize by the addition of the enhancement factor, 
J-", of the reflected continuum (only) compared to that expected from a semi-infinite plane 
illuminated by a point source viewed face-on. 

We find fixing JF = 1 does not improve the goodness-of-fit at > 99% confidence for any 
of the datasets. As expected, in most cases the derived spectral index is slightly steepened 
by AF ^ 0.15, and in all cases best-fitting values of Nh,z and Ux are consistent with those 
quoted in Table 8. 

Interestingly, allowing to be a free parameter in the analysis, we do find a significant 
improvement in the goodness-of-fit for 9 datasets (Mrk 335, NGC 3783(2), NGC 4051, 
NGC 4593, MCG-6-30-15(l,2), IC 4329A, NGC 5548, Mrk 509), with -3-5. Such 
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values have been suggested for some of these objects by previous studies using different 
instrumentation (e.g. Nandra & Pounds 1994; Molendi et al, 1997). However, since the 
ASCA bandpass effectively ends at 10 keV, such a component is difficult to constrain 
and weaker Compton-refiection components would be undetectable. For JF <^ 10, the 
predicted equivalent width of the Fe Ka from the reflecting material can be approximated 
by Wko ~ 150J-'( 3 3^^Q_5 ) eV, where Ap^ is the abundance of Fe relative to hydrogen 
(e.g. George & Fabian 1991). The most recent estimations of the Fe abundance are 
Ape ~ 4.7 X 10"^ (Anders & Grevesse 1989) and ~ 3.2 x 10"^ (Feldman 1992), and 
from observations we find Wxa — 290 ± 30 eV ( §6.3. 1|) for these 9 datasets. Thus these 
best-fitting solutions imply a moderate under-abundance of Fe in these Seyferts compared 
to solar values. We do not consider this unreasonable per se, but note these solutions may 
also be indicative of curvature in the underlying continuum shape, rather than the presence 
of a strong reflection component. Indeed alternative explanations are available in the case 
of Mrk 335, NGC 4051, IC 4329A (see § pX^ . In any event, we find the inclusion of the 
reflection component typically only steepens the underlying spectral index by AF ^ 0.2 
in these 9 datasets, and in all cases again find the best-fitting values of Nh,z and Ux are 
within the error ranges quoted in Table 8. 



6.4-2. Double powerlaw continua 

The form of the continuum considered in this section is represented by the sum of two 
powerlaws0, with photon indices F^ and Th (where F^ > Th). The ratio of the normalization 
of the soft powerlaw divided by that of the hard powerlaw at 1 keV (in the rest-frame of 
the source) is denoted by Rs/h- For computational expediency we restrict parameter space 
such that 1.5 > F^ > 5.0 and 0.0 > Th > 2.5. Such a model has been applied to all the 
datasets considered here. Here we concentrate on the results when this additional powerlaw 
component is added to model C(i) as they are representative of the rest. Clearly such a 
model has the same number of interesting parameters as the model discussed in §|6.4.1. 



We find this model leads to an improvement over the results obtained assuming model 
C(ii) for 10 of the 23 datasets at >99% confidence. The best-fitting solutions for 4 of these 
datasets have Rg/h ~ 1 with both powerlaws of similar importance over a sizable fraction 
of the ASCA bandpass. In the case of the remaining 6 datasets, the bulk of the ASCA 



^^i.e. not the somewhat artificial 'broken-powerlaw' model (which has also been applied in the past to 
AGN spectra by some workers) in which the continuum is assumed to be a steep powerlaw below some energy 
at which it sharply breaks to a flatter powerlaw 
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bandpass is dominated by a single powerlaw, with the second powerlaw component making 
a significant contribution only at either the very lowest energies (the 'ultra-soft excess' 
datasets) or highest energies (the 'hard tail' datasets). 

We consider first the 3 datasets (3C 120, NGC 4593, Mrk 509) for which this model 
implies a 'hard tail'. Naturally we find the powerlaw which dominates the observed 
spectrum at the highest energies to be relatively weak {Rg/h >> 1) and flat (although its 
photon index is poorly constrained within the range 0.2 ^ <^ 1.8). In each case the 
steeper powerlaw. dominating the observed spectrum over most of the ASCA bandpass, has 
a photon index steeper by AF ~ 0.2 than that found assuming model C(ii) (in fact all 3 
datasets have Tg — 2.2). However, it should be noted that the underlying continua derived 
for these 3 datasets assuming this model exhibit only a relatively subtle curvature across 
the ASCA bandpass. It is therefore not that surprising that in all cases the best-fitting 
values of the parameters associated with the ionized gas {Ux and Nh,z) are similar to those 
found when the underlying continuum is approximated by a single powerlaw. 

Turning now to the 3 datasets for which this model implies an 'ultra-soft excess' 
(NGC 3516, MCG-6-30-15 (1,2)), in all cases we find the powerlaw which dominates the 
observed spectrum over most of the ASCA bandpass has a photon index fiatter by AF ~ 0.1 
than that found for model C(ii). The best-fitting solutions imply a very steep second 
powerlaw component (pegging at our upper limit of F^ = 5 for all 3 datasets) dominating 
the observed spectrum only at the very lowest energies (with Rg/h ~ 0.3-0.5). However, 
none of the best-fitting models extrapolate <0.6 kcV in a satisfactory manner (even given 
best- fitting values of Nh,o significantly in excess of Nf^^^), and all imply an aox far steeper 
than that observed. Further investigation of such a component is extremely difficult with 
ASCA data. However, again, the best-fitting values of the parameters associated with the 
ionized gas {Ux and Nh,z) are similar to those found when the underlying continuum is 
approximated by a single powerlaw. 

Finally, we consider the remaining 4 datasets for which this analysis indicates both 
powerlaws are of similar importance over a sizeable fraction of the ASCA bandpass 
(Mrk 335, NGC 4051, IC 4329A, NGC 5548). The best-fitting parameters of for these 
datasets are listed in Table 11, and the derived model spectra and data/model ratios shown 
in Fig 15. As further discussed in the Appendix, the form of the continuum implied for 
Mrk 335 is in good agreement with previous observations. However any imprints due to 
ionized gas are weak in Mrk 335, allowing no useful constraints to be placed on A^^^^, 
Ux or Q. The best-fitting solution for the other 3 datasets are somewhat similar, with a 
relatively steep (F^ ^ 3) powerlaw dominating the spectrum ^ 1 keV along with evidence 
for absorption and emission from ionized gas. However, in all 3 cases the reality of such a 
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continuum can be questioned. Significant spectral and flux variability occur on timescales 
much shorter than the observation period in the case of NGC 4051 (Paper I). Thus we 
consider short-timescale spectral variability may be partly responsible for the apparent 
curvature of the underlying continuum in NGC 4051. As is apparent from Table 11, this 
model fails to satisfy our criteria for acceptability in the case of IC 4329A. As further 
discussed in the Appendix, this source appears to have a particularly complex spectral 
form, and hence it is difficult to judge the robustness of the values derived for Nu^^i Ux or 
Vt. In the case of NGC 5548, an acceptable fit is achieved assuming this model. However 
it should be noted that the steep {Tg ~ 3.2), soft component derived does appear to be a 
relatively poor representation of the data < 0.6 keV, implying a neutral column density far 
in excess of > N^^fi (A^h,o ~ 1-5 x cm-^) and giving ^f-^ ~ 14 (though R^q ~ 1). 
Finally, as shown in § |6.4.1| , a single powerlaw with a strong Compton-reflection component 
can provide a comparable flt for all 3 datasets. 

6.4-3. The effects of including a low-temperature, thermal component 

It has long been suggested that Seyfert galaxies are likely to contain several sites which 
give rise to thermal emission at temperatures consistent with soft X-rays. These sites 
occur on different size-scales and arise as a result of processes both directly related to and 
independent of the AGN phenomenon. Examples include regions of active star formation, 
the hot-phase of the narrow emission-line region (NELR), the halo of the galaxy and 
the shock- or X-ray heated ambient interstellar medium (ISM). Thermal emission closely 
associated with the active nucleus is also expected as a result of the accretion process 
itself and as a result of the absorption (followed by thermalization) of the continuum. The 
relative importance of such components probably differs from object to object. In the case 
of NGC 4151, Einstein HRI and ROSAT observations have revealed spatially resolved 
soft X-ray emission (Elvis, Briel & Henry 1983; Morse et al 1995). ASCA observations 
of normal galaxies have also revealed soft X-ray emission with temperatures in the range 
kT ~ 0.3-0.6 keV, and luminosities Lkx ~ 10^^-10^°erg s"^ (e.g. Makishima 1994; 
Serlemitsos, Ptak & Yaqoob 1996). Thus one could question whether the existence of such a 
component might effect our conclusions regarding the photoionized gas - particularly since 
the thermal emission offers an alternative to the decreased opacity and emission features at 
the lowest energies implied in previous models. 

We have investigated this by including a component due to an optically-thin plasma 
(with cosmic abundances, parameterized by a temperature kT and unabsorbed luminosity 
in the 0.1-3.0 keV band) to model C(ii) (but located outside the photoionized 
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absorber). We find that the goodness-of-fit is indeed improved (A^^ ~ 9) in the case of all 
3 datasets of NGC 4151, with kT ~ 0.66 keV and LkT ~ 2^1 x 10^^ erg s'^ consistent with 
the luminosity of the extended gas measured by the ROSAT HRI observation (see §^]2| and 
Appendix). However, it should be noted that in the ASCA bandpass, the dominant features 
from such gas are Fe L-shell lines in the 0.7-1.5 keV regime. Thus whilst L^t is a function 
of Ape, the addition of this optically-thin component does not significantly affect the derived 
values for all other of the parameters (although the size of the allowed parameter space is 
obviously increased - e.g. see Fig 16). 

We find no other dataset is significantly improved by the inclusion of emission from 
thermal gas at such a temperature. This is hardly surprising since only NGC 4151 occupied 
a region of the Nh,z^Ux plane in which emission of similar strength would not be swamped 
by the underlying continuum. Thus we consider the the omission of such a component is 
generally not important in our results. 



6.4-4- Summary 

In summary, we find evidence for significant curvature in the observed continuum in 
approximately half of the datasets considered. We have shown that for the majority of the 
datasets, such curvature can be accounted for by the inclusion of a Compton-Refiection 
and/or a second powerlaw component. However, in the majority of cases, the ionizing 
continuum can be well represented by a single powerlaw (with F ~ 2) close to the 
Ovil and Ovill absorption edges. These edges provide the main diagnostics for our 
photoionization models. Thus we consider that the best-fitting parameters associated with 
the ionized gas derived in models B and C are generally reliable. Allowing the presence 
of a low-temperature, thermal component improved the fit only in the case of NGC 4151. 
However, again we found it to have no significant effect on the best-fitting parameters 
associated with the ionized gas. 



7. OBJECTS WITH MULTIPLE OBSERVATIONS 



Our sample contains more than one observation of NGC 3783, NGC 4151, MCG-6-30-15 
and Mrk 841. We have found evidence for absorption by ionized material in all four sources, 
enabling us to compare and contrast the behaviour of the ionized gas. 
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7.1. NGC 3783 

Two datasets of NGC 3783, separated by 4 days, are contained within our sample 
from observations carried out in 1993 December. As shown in § |5.3| - |5.6| , there is clear 
evidence for ionized gas in this source, with models B(i)-C(ii) providing acceptable fits 
to both datasets. The observed count rate was ~ 25% higher for the second observation 
(e.g. Table 2), which is reflected in an increase in the derived luminosity of the ionizing 
continuum (for model C(i), after correcting for absorption) from Lx — 4.6 x 10^^ erg s^^ to 
^ 5.5 X 10^^ erg s^^. Both observations are consistent with F ~ 1.8, Df = 0, Ux ~ 0.1-0.15 
and Nh,z ~ 2 x 10^^ cm~^ (Table 8). However, as shown in Fig 16, the two datasets actually 
occupy slightly different regions of the Nh,z^Ux plane, consistent with an increase in Ux by 
^ 40% and/or a decrease in Nh,z (by ^ 5 x 10^^ cm~^) between the 2 epochs. The former 
offers the simplest interpretation of these data: a screen of ionized gas completely covers 
the cylinder-of-sight in NGC 3783, becoming more highly ionized as the photoionizing 
flux increases. Both observations are consistent with the presence of emission from the 
ionized material, with Le ~ 2-6 x 10^^ erg s~^ in the 0.1-10 keV band (after correcting for 
absorption). Unfortunately however, the signal-to- noise ratio of these datasets is too low to 
unambiguously determine whether the intensity of this emission responded to the increase 
in the ionizing continuum. 

Interestingly, a '1 keV deficit' was also evident for both datasets in the data/model 
residuals of all the models considered in §|[ This can be modeled by an additional 
absorption edge {tons ~ 0.5) at an energy close to that of Ovill (§ |6.2|) and most likely due 
to an additional screen of more-highly ionized gas covering part, or all, the cylinder-of-sight. 

7.2. NGC 4151 

Three datasets of NGC 4151 fulfilled the criteria defining this 'sample' (Paper I). These 
were an observation performed 1993 Nov 05 (NGC 4151(2)), an observation performed 
exactly 30 days later (NGC 4151(4)), and an observation performed 2 days thereafter 
(NGC 4151(5)). The observed count rate for NGC 4151(4) was only ~ 2% higher than that 
for NGC 4151(2), but that for NGC 4151(5) was ~ 15% higher. As shown in Table 12, we 
find a fit satisfying our criteria for acceptability is achieved for the first two observations 
(NGC 4151(2,4)) only when the absorbing material is allowed to be ionized and a fraction 
of the underlying continuum is allowed to escape without suffering any attenuation (i.e. 
models B(ii) & C(ii)). Such an hypothesis is strongly supported by the third observation 
(NGC 4151(5)), even though none of the models in §^ provide a solution which formally 
satisfies our criteria for acceptability (see Fig. 9). 
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As can been seen from Tables 6 and 8, all 3 datasets are consistent with -D/ ~ 3-4% 
and r ~ 1.5. However the parameters associated with the ionized gas vary in an interesting 
manner (see also Fig 16). In the month between the NGC 4151(2) and NGC 4151(4) 
observations, there is an increase in the absorbing column density {ANh,z ~ 2 x 10^^ cm~^) 
and increase in ionization parameter {AUx ~ 30%), whilst the derived luminosity of the 
ionizing continuum (in the 0.1-10 keV for model C(ii), after correcting for absorption) 
increases by ^ 10% (from Lx — 1-3 to 1.4 x 10^^ erg s~^). Two days later, at the time 
of the NGC 4151(5) observation, Nh,z is similar to NGC 4151(4). However the ionization 
parameter has decreased {AUx ^ 15%) despite an increase in the luminosity of the 
ionizing continuum (to Lx ~ 1.8 x 10^^ erg s~^). Such behaviour is difficult to understand 
in terms of a simple model, unless the Nh,z and/or the number density of the gas uh, 
(but not Df) are allowed to vary on timescales ^ 2 days, or unless the photoionized 
gas is out of equilibrium. A simple explanation of such behaviour is that the absorbing 
gas completely covers the cylinder-of-sight but is inhomogeneous (with differing Nh,z 
and/or nn) on a scale comparable or greater in size to that of the central source, whilst 
a constant fraction -D/ of the underlying continuum follows an alternate path. If the 
transverse velocity of the absorbing material is sufficiently high, then the movement of 
such clumps across the cylinder-of-sight would then give rise to variations in Ux, Nh,z 
seemingly decoupled from variations in the intensity of the underlying source. Assuming 
only Keplerian motion, the length-scale / passing through a given line-of-sight in a time 
t is / ~ 2 X 10^ {fboix / fEdd)^'^Lx\4rl~f'^t cm at a distance rid light-days from a central 
object with a bolometric luminosity Lhoi = fboixLxu emitting at a fraction fEdd of its 
Eddington luminosity, where t is in seconds, Lxm = lO'^'^Lx erg s~^ and Lx defined over 
the 0.1-10 keV band. Equating t with the timescale tabs on which significant variations in 
Nff^z are seen, and putting I > ct^ar, the size of the central source, we obtain 

rid ^ 4.5 X 10'^ {fboix/fEdd)LxM{tabs/tvarf light - days (2) 

The ratio fboix/ fsdd is generally considered to be ~few, and from our analysis we find 
Lx44 — 0.15 for NGC 4151. Considering the variability between the 4151(2) and 4151(4) 
datasets and thus setting tabs ~ 30 days, we obtain the limit rid ^ ^^{fboix/ fEdd) light-days 
assuming t^ar ~ 10^ s. The limit if the same explanation is proposed to explain the 
variations between the 4151(4) and 4151(5) datasets {tabs ~ 2 days) is a factor ~ 200 
smaller. We stress that these limits are based only upon the Keplerian motion, which may 
not be the major component governing the dynamics of the photoionized clouds. 

In § |6.4.3| we found that there was evidence for an additional spectral component, 
which we tentatively identified with the extended soft X-ray emission seen in this source 
and modeled as an optically-thin plasma. We found that all 3 datasets are consistent with 
such gas with kT ~ 0.7 keV and L^t ~ x lO^^erg s~^. Such a component alone (but 
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absorbed by Nf^^) would give rise to a count rate in the range 10"-''-0.07 count s'^ for the 
ROSAT HRI, and hence is consistent with that reported for the extended emission in this 
source (0.06 count s~^, Morse et al 1995). However, the inclusion of this component did 
not significantly affect the values derived for the ionized gas (dashed contours in Fig 16), 
and hence the behavior of the source between the 3 observations still suggests clumps of 
gas of differing N^^z traversing the cylinder-of-sight. These results are further discussed in 
relation to previous studies in the Appendix. 



7.3. MCG-6-30-15 

Two datasets of MCG-6-30-15 are contained within our sample from observations 
carried out in 1993 July, separated by ~ 3 weeks, with the observed count rate ~ 52% 
higher for MCG-6-30-15(l) (Table 2). As shown in § p.3H5.6| , there is clear evidence for 
ionized gas in this source, with models B(i)-C(ii) providing acceptable fits. In all cases 
a '1 keV deficit' was also evident, which can be modeled by an additional absorption 
edge (§16:^). The best-fitting continuum for MCG-6-30-15(l) has T ~ 2.1, whilst that 



for MCG-6-30-15(2) is slightly flatter (AF ~ 0.1). Both observations are consistent 
with Df = 0. As shown in Fig 16, there is clear evidence for an increase in Nh^z (by 
~ 1-6 X 10^^ cm~^) over the 3 weeks between the observations. However, the ionization 
parameter remains constant, or perhaps even increases (but with AUx ^ 50%), between 
the 2 observation despite a decrease in the luminosity derived for the ionizing continuum 
from Lx — 4.5 x 10^^ erg s~^ to ^ 3.0 x 10^^ erg s~^ (for model C(ii), after correcting for 
absorption). Such behaviour has been reported previously for these datasets by Fabian 
et al (1994), and as in the case NGC 4151, is difficult to reconcile with that expected 
in the simplest picture whereby a uniform shell of gas reacts to changes in the intensity 
of the photoionizing continuum. Again appealing to only Keplerian motion to move 
inhomogeneities through the cylinder-of-sight between the 2 observations (i.e. setting 
tabs ^ 3 weeks) , and again setting the variability t^ar ~ 10^ s, from Eqn |^ one obtains the 
rather weak constraint r;^ ^ 5 x lO"^ {fboix / fEdd) light-days. Again we note that these limits 
are based upon only the Keplerian motion, which may not be the major component of the 
dynamics of the ionized gas. 

However, there is evidence in a subsequent ASCA observation that whilst the depth 
of the Ovil edge appears to be constant on long timescales in this source, the depth of 
the Ovill edge varies on a timescale ~ 10^ s (Otani et al 1996). Such behaviour has been 
recently confirmed in BeppoSAX LECS observations (Orr et al 1997). Setting tabs ~ 10'^ s in 
Eqn Hone obtains ria ^ ^{fboix I fEdd) light-hours. However, since the photoionization and 
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recombination timescales for these ions are similar, such behaviour is incompatible with the 
observed Ovil and Ovill edges being produced in the same region of ionized gas. Otani 
et al suggested the most natural solution was a two-zone model in which the gas giving 
rise predominately to the Ovil edge is stable, whilst that responsible for the Ovill edge 
undergoes rapid variations. Such an hypothesis is of course supported by our observations 
of a (relatively small) '1 keV deficit' in the data/model residual of all the models considered 
in In §^]2| we showed that this can be modeled by an additional absorption edge 
{ToNe ~ 0.14) at an energy close to that of Ovill (see § |8.8|) . 

Both observations of MCG-6-30-15 are consistent with the presence of emission from 
the ionized material, with Lg <^ 10^^ erg s~^ in the 0.1-10 keV band (after correcting for 
absorption). However, the signal-to- noise ratio of these datasets is too low to determine 
unambiguously whether the intensity of this emission responded to the decrease in the 
ionizing continuum. 

7.4. Mrk 841 

Two datasets of Mrk 841 are contained within our sample from observations carried 
out in 1993 August and 1994 February. The observed count rate was (~ 33%) higher for 
Mrk 841(1) (e.g. Table 2). In §^we found evidence for ionized gas in this source. For models 
B(i)-C(ii), the underlying continuum for Mrk 841(1) has F ~ 1.9, but the values derived 
for Ux and Nh,z are rather model dependent and often poorly constrained. In contrast, 
the data from Mrk 841(2) are consistent with a simple powerlaw, preventing any strong 
constraints being placed on Ux and Nh^^- In all cases, the region of the Ux-Nh,z plane 
occupied by Mrk 841(1) for a given model is consistent with the limits from Mrk 841(2) 
assuming that model. However the underlying continuum for Mrk 841(2) does appear to be 
significantly flatter (F ~ 1.7). Dramatic spectral variability has been observed previously 
in this source (see Appendix and references therein). The luminosities derived for the 
ionizing continuum are Lx — 1.6 x 10'^^ erg s~^ and ~ 1.1 x 10^^ erg s~^ for Mrk 841(1) 
and Mrk 841(2) respectively (for model C(ii), after correcting for absorption). In view of 
the poor constraints which can be placed on the ionized gas from these data, few further 
comments can be made. 
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8. DISCUSSION & CONCLUSIONS 

As discussed in the introduction, it is becoming increasingly clear that the spectra of a 
number of individual AGN contain features indicative of absorption by ionized gas within 
the cylinder-of-sight. It seems highly likely that this gas is photoionized by the intense 
radiation field produced by the central object. The aim of the work presented here was 
to determine the frequency and characteristics of such gas based upon new, self-consistent 
modelling using the ION photoionization code. 

In the preceeding sections we have presented the results from 23 ASCA observations of 
18 objects. When considered together these objects certainly do not constitute a complete 
sample. Nevertheless, in at least some aspects, we do believe our results to be representative 
of the variety of properties exhibited by the type-1 AGN as a whole. 

In §^ we considered models assuming the underlying continuum in the 0.6-10 keV 
band of our sources is well-represented by a single powerlaw. A number of models were 
considered, starting with the case in which the spectrum emerging from the source passes 
through two screens of neutral material fully covering the source: one at the redshift of the 
source and a second at zero redshift (with column densities Nh,z and Nh^ respectively). 
Models of increasing complexity were then considered by allowing a fraction (-D/) of 
the underlying continuum to be observed without attenuation by Nh,z, allowing the gas 
represented by Nh,z to be photoionized (with an ionization parameter Ux), and finally by 
also including the emission spectrum expected from the ionized gas (assuming it subtends a 
solid angle fl at the central source). 

In §^ the results of further spectral analysis were presented, including models with 
additional spectral components, with the primary goal of testing the robustness of the 
derived characteristics of the ionized gas. 

8.1. General Considerations 

All models considered in this work are single-zone models, and given the assumed gas 
density uh = 10^° cm~^, the gas reacts instantaneously to continuum variations. However, 
the models are intended to represent a wide range of conditions including cases with much 
lower values of uh- In such cases, the level of ionization represents the response of the gas 
to the long-term, average flux of the illuminating source. We have not considered all such 
cases here since we do not consider the data yet warrant this level of discussion. 

In §|5| we found that one or more of the models considered offer an explanation of 
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17/23 datasets within the formal criteria outhned in §0. The criteria used to determine 



the acceptabihty of model as a description of the data is of course a somewhat subjective 
decision. The primary criterion used here that P(x^ | dof) < 0.95 (i.e. that the model is 
a true representation of the data with a probability of better than only 1 chance in 20) 
might be considered relatively lax by some readers. Furthermore, as can been seen from 
Tables 3-8, a 'perfect' parameterisation of the data (-P(x^ I dof) ~ 0.5 corresponding to 
a reduced-x^ value of unity) is only obtained in a small number of cases. Relatively high 
values of P{x^ \ dof) imply that there are residuals in the data not modelled adequately. 
However, whilst some of these residuals might be carrying physical information associated 
with the source, others might be the result of inaccuracies in the calibration of the 
instruments (see also Fig. 13). The latter is currently under intense study by the instrument 
teams in anticipation of being able to assign unmodelled residuals to the former explanation 
with a higher degree of confidence. 

However our intention in this paper is to explore just how well the data from the 
sample can be described by our relatively simple models for the ionized gas. We believe this 
approach compliments that taken by some previous workers where the absorption features 
introduced by the ionized material are modelled as a series of absorption edges. Throughout 
the analysis we have not assumed any systematic errors, and have implicitly assigned the 
relatively poor values of P{x^ \ dof) to be the result of problems with the calibration in the 
case of 4 additional datasets in §^ (Table 12). Furthermore, beyond the limited number 
considered in §|^, we have not explored more complex models in detail. For instance we 
have not considered models in which the ionized material has a significant {'^ 10^ km s~^) 
velocity dispersion along the line-of-sight. Given the spectral resolution of the SIS detectors, 
dispersions smaller than this value are impossible to constrain with the signal-to-noise ratio 
of the data presented here (see also §|6.1|). We have not explored different abundance ratios, 
or multi-component models in those cases where we have indications for their existence (see 
§8.8| ). Our intention in the present work is merely point the reader to the need of such 
analysis and defer it for future discussion. 



8.2. The form of the underlying continuum 

It is well established that many Seyfert 1 galaxies exhibit so-called 'soft-excesses' and 
'Compton-reflection' tails (e.g. Turner & Pounds 1989; Nandra & Pounds 1994). After 
accounting for these effects, and those of the ionized-absorber, Ginga observations showed 
that the underlying X-ray continua can be generally well represented by a powerlaw of 
photon index F ~ 1.9-2.0 (e.g. Nandra & Pounds 1994). In we obtained acceptable fits 
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for most of the datasets by one or more of the models considered, with the distribution of 
the photon index peaked at F ~ 2 (Fig. 2). In §6.4| we found significant improvements in the 



goodness-of-fit were obtained for many of the datasets by the inclusion of additional spectral 
components. However, in the majority of cases, we again found the bulk of 0.6-10 keV 
bandpass to be consistent by an underlying powerlaw with F ~ 2. It should be noted that in 
a number of cases the best-fitting models imply that very little of the underlying continuum 
is observed directly (e.g. see NGC 3783(1,2) in Fig 6; NGC 4151(2,4,5) in Fig 9). In passing 
we note that in such cases, F can still be well constrained by the data presented here as 
the shape of the 'recovery' in the absorption features in the 1-5 keV band is well defined 
in the models and well constrained by ASCA data. For comparison with other work, we 
have calculated the expectation value and dispersion of the distribution of spectral indices, 
employing the method of Maccacaro et al. (1988). Considering those datasets for which we 
obtained acceptable fits, we find < F >= 1.94 ± 0.08, with a significant, intrinsic dispersion 
of 0.18 lb 0.07. These values are in excellent agreement with those found by Ginga and in 
Paper II, when the effects of Compton reflection are taken into consideration. 

In a number of cases, complex continua do seem to be present in the ASCA bandpass. 
Mrk 335 does indeed appear to possess real curvature in the underlying continuum (Fig 15). 
Evidence was also found for curvature in the observed continuum in NGC 4051, IC 4392A 
& NGC 5548. However in at least NGC 5548 we consider the evidence for curvature in 
the underlying continuum to be questionable and/or alternative explanations are available 
(see § |6.4.2|) . The observed spectrum of NGC 4151 has long been considered somewhat 
idiosyncratic amongst Seyfert 1 galaxies. However we do find acceptable solutions and 
consider the underlying continuum to indeed be a powerlaw, albeit flatter (F ~ 1.5) than 
that observed for the majority of sources (see and § ^.7| ). 



8.3. The frequency of neutral absorbers 

Unlike Seyfert 2 galaxies, Seyfert 1 galaxies are generally not believed to contain 
significant column densities of neutral gas detectable in the 0.6-10 keV band (ie. within the 
range 10^° ^ Nh^z ^ 10^^ cm~^). In the case of models A(i) and A(ii), the absorption at the 
redshift of the source (with column density Nh^z) is assumed to be due to neutral material. 
As described in § |5.1| , significant absorption is not detected in any of the datasets for which 
model A(i) offers an acceptable solution (with Nh,z ^ few xlO^° cm~^; see Fig 3a). Model 
A(u), in which a fraction Df of the underlying continuum is able to escape without passing 
through Nh,z, does offer an improvement in the goodness-of-fit for many datasets, and a 
range of derived column densities for those in which an acceptable solution was obtained 
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(Fig 3b). However, as discussed in § ^.2| , these solutions are most likely due to a residual 
problem in the instrument calibration (e.g. 3C 120), lie in a region of the N^^zyDf indicative 
of (relatively-subtle) spectral curvature (e.g. NGC 7469), and/or superior solutions are 
obtained assuming subsequent models (e.g. Mrk 335). We note that it has long been 
suggested that NGC 4151 might contain a significant column density attenuating ~ 95% of 
the continuum. However as shown in we find acceptable solutions only for models in 
which the gas is ionized (see also §0 and §|8l^). 

We have not explicitly tested for additional intrinsic column densities of neutral 
material in any of our models which contain ionized gas (i.e. in models B(i)-C(ii)). Such a 
situation might arise for example in the case where the primary continuum passes through 
a region of ionized gas, and then a further screen of neutral gas during its passage through 
the host galaxy. Nevertheless, all our models do contain a neutral column at zero redshift, 
Nh.q, to account for the Galactic absorption (and hence expected to be ~ Nf^^). At the 
redshift of most of the sources in our sample, any absorption by neutral material intrinsic 
to the source would be indistinguishable from absorption at zero redshift. Thus any such 
absorption would result in a derived value of N^fl » ^f^o- Generally we find no such 
evidence in our sample, with the only reliable detection in the case of IC 4329A (with a 
neutral column density ~ 4 x 10^^ cm~^). The host galaxy is this AGN is observed edge-on 
so this result is not surprizing (see Appendix). 



8.4. The frequency of ionized absorbers 

In §^ we found 16/23 datasets (11/18 objects) are improved at > 99% confidence 
over models A(i) & A(ii) if the absorbing material is assumed to be photoionized. Of 
these, 12/16 datasets (9/11 objects) can be adequately described by one or more of models 
including ionized gas presented in §^ (i.e. models B(i)-C(ii)), under the assumption of 
a single power-law continuum. Of the remaining 4/16 datasets, 3 datasets (only one of 
which being from an additional object: NGC 3516) show strong evidence for ionized gas 
(with the statistical poverty of the fit due to further spectral complexity). Thus we find 
clear evidence for absorption by ionized gas in 15/23 datasets (10/18 objects). However 
the presence of ionized gas is also strongly suspected in 3 additional datasets (3 additional 
objects: NGC 4051, IC 4329A, NGC 7469) from the analysis presented in §§ 

• Thus we conclude there is evidence for ionized gas in 18/23 datasets (13/18 objects) 
in our sample. 
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Reynolds (1997) has also presented results from a sample of 24 type-1 AGN using 
ASCA data, using the same datasets as presented here for 15 sources (Tabic 2). Based on a 
search for Ovil and Ovill edges, Reynolds found strong evidence for ionized gas in 12/24 
objects. These 12 include 9 sources also in our sample (Table 12), along with Mrk 290, 
3C 382 and MR 2251-178. Reynolds also fitted these data with a single-zone photoionization 
model (based on the photoionization code CLOUDY which is similar in approach and scope 
to ION used here), assuming a photoionizing continuum consisting of a single power law over 
the entire 0.0136-13.6 keV band. Converting their quoted ionization parameters (^) to Ux 
using Fig.lb, gives 1.5 ^ T ^ 2.0, lO^^ ^ Nh,z ^ 2 x lO^^ ^.^-2 q q2 ^ Ux ^ 0.09, 
broadly in agreement with the distributions found here (Figs 2, 3 and 7 respectively). A 
detailed source-by-sourcc comparison between our results and those of Reynolds (1997) and 
other workers is provided in the Appendix. 

8.5. Column density and ionization parameter of the ionized gas 

We find no preferred value of Nh,z under any of the models tested, with clear differences 
from object-to-object (e.g. Fig. 3). However, in each of the models including photoionized 

gas, wc find the distribution in the ionization parameter strongly peaked at Ux ~ 0.1 
(Fig. 7). Considering those datasets for which we obtained acceptable fits (and again 
employing the method of Maccacaro ct al. 1988), we find < logUx >= —0.92 ±0.21. While 
the values cluster around this value, there is a highly significant intrinsic dispersion in 
logUx of 0.2llo;o9. Nevertheless this relatively narrow range in Ux is our most interesting 
and unexpected result. As such it deserves some further discussion as the whether it is 
purely the result of selection effects, or has physical significance. 

We considering first values of Ux » 0.1. With the signal-to-noise ratio typically 
afforded by ASCA data, the presence of ionized gas is most readily inferred by deep 
(r ^ 0.2) absorption edges which it introduces into the observed X-ray spectrum. For a 
given value of Nh,z, there is a maximum value of Ux above which the gas is unstable and 
goes to the high electron temperature (Te ~ 10^ K) branch of the two-phase curve (e.g. 
sec Fig. 2 in Netzer 1996). The value of Ux at which this thermal instability occurs is a 
function of the gas density and composition, along with the form of the ionizing continuum. 
Under the assumptions made here, the instability occurs for Ux > 10^^^ A''//.^. Such gas 
is completely transparent, and the observed spectrum will be identical to the underlying 
(powerlaw) continuum, consistent with the data from a minority of objects in our sample. 
Therefore the presence of such material along the cylinder-of-sight is impossible to prove 
(or disprove) using the current ASCA data alone, offering a potential explanation of the 
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apparent lack of objects with Ux » 0.1 in Fig. 7. 

For Ux << 0.1, different ions contribute different continuum opacity, depending on the 
gas composition and level of ionization. For 10^^ < Nh,z < 10^^ cm~^ and 0.03 < Ux < 0.3, 
most of the opacity is due to Ovil and Ovill (Fig 8). For lower values of Ux, carbon, 
nitrogen and lower ionization states of oxygen are more important. However, we note that 
none of the objects presented here lay in the region of Ux^Nh^z parameter-space where 
Ux < 0.03, Nh,z > 10^^ cm~^. There is no reason a priori why there should not be objects 
with ionized material in this region. However, there is a potential selection-effect here 
whereby if such material lays along the cylinder-of-sight to the broad emission line region 
(BELR) and contains embedded dust, the broad emission lines will be absorbed by this 
material. This may lead to such objects being classified as something other than Seyfert 1 
galaxies, and hence being excluded from the sample of sources presented here. Alternatively, 
should this explanation not be the case, then there must be some yet unknown mechanism 
that determines these unique conditions by either fixing the pressure or perhaps the 
dynamical state of the absorbing material. Future studies of a larger sample of objects are 
required to investigate these possibilities. 

We have compared the parameters of the ionized gas with the derived continuum 
luminosity, Lx, in the 0.1-10 keV band (after correcting for absorption). We find no clear 
relation between Nh,z and Lx, but some indication that Ux oc Lx (Fig. 17). If true, from 
the definition of Ux (Eqn. |l|), this implies the ionized gas has similar densities, uh, and 
is at similar radii, r in all sources. However, the majority of the sources considered here 
for which Ux is well constrained lie within a restricted range of Lx (10^^-10^^ erg s~^). A 
recent observation of the quasar PG 1114-1-445 {Lx ~ 6 x 10^^ erg s~^) also exhibits strong 
absorption features due to ionized gas with Nh,z — 2 x 10^^ cm~^ and Ux — 0.1 (George 
et al. 1997b). Furthermore, Nandra et al (1997e) have recently reported the results from 
ASCA observations of the RQQ MR 2251-178 {z = 0.068, Lx ^ 2 x 10^^ erg s'^), finding 
Nh,z ~ 2 X 10^^ cm~^ and Ux ~ 0.07. Inclusion of these objects on Fig. 17 argues against 
any obvious relation between Ux and Lx- 

We found no evidence for any relationship between any other pairs of parameters 
derived from the X-ray results (e.g. Lx vs F; F vs Ux', F vs Nh,z)- In most cases this is 
again primarily due to the limited range of Lx and/or F exhibited by most the objects 
in the sample. Thus, further progress requires high signal-to-noise observations of similar 
sources covering the high and low Lx regimes. 
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8.6. Constraints on ionized emitters 

The dominant effects of the ionized gas on the observed spectrum from these sources 
is absorption by the material along the line-of-sight. However, ionized material out of the 
line-of-sight will give rise to emission lines and recombination continua, along with the 
(absorbed) continuum Compton-scattered back into the line-of-sight. The inclusion of such 
an emission component leads to an improvement at >99% confidence in 11 datasets (9 
objects) in the full covering case (model C(i)), and 8 datasets (7 objects) in the partial 
covering case (model C(ii)). However there is no case for which this component is required 
(i.e. a satisfactory fit only obtained for model C). It should also be noted that the inclusion 
of the emission component generally does not have a significant effect on the values derived 
for Ux and Nh,z- For the value of Ux and Nh,z derived for the majority of sources, the 
emitted spectrum contains significant line emission in the 0.6-2.0 keV band, but constitutes 
only a relatively small fraction (^ 10%) of the total flux observed in this band, even when 
the ionized gas subtends a large solid angle at the central source. The best-fitting values of 
Q are ill-constrained, but in most cases consistent with Q = An, and there is no case where 
r2 > 47r at > 95% confidence. Formally we find < logr2/47r > = 0.23lo[23 the datasets 
for which we obtained acceptable fits. 

In cases where the ionizing continuum source is constant, then the normalization of the 
emission spectrum is a direct measure of the geometry of the emitting gas. However in cases 
where the ionizing continuum varies, then the normalization of the emission spectrum will 
follow the luminosity history of the continuum source, with a lag due to light travel-time 
effects. In the latter case we are therefore offered the opportunity to determine the location 
of the ionized gas via the delay in the intensity of the emission features following variations 
in the ionizing continuum. Unfortunately however, the signal-to-noise ratio is too low 
for the few objects presented here for which there are multiple observations (^) to place 
even crude constraints on whether the intensity of the emission changed in an appropriate 
manner. 

8.7. Evidence for an unattenuated component 

In §^ & §^ we find that allowing a fraction -D/ of the powerlaw continuum to be 
observed without suffering attenuation by the column density Nh,z significantly improves 
the fit (compared to the corresponding model with Df = 0) for a large number of 
dataset/model combinations. Indeed, inspection of (Table 4) reveals all but 3 

datasets are improved in the case of model A(ii), and inspection of F(^) and F{^) 
(Tables 6 and 8) reveals that 8 datasets are improved in the case of models B(ii) and C(ii). 
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However, in many of these cases the best-fitting solutions have Dj ~ 1 (Fig. 5) and as 
described in §|6.4| are an artifact of spectral curvature. Indeed, as shown in Table 12, only 



in the case of 2 datasets (1 object: NGC 4151) are solutions which satisfy our criteria for 
acceptability obtained only if Df > 0. Thus we conclude most of the sample show solutions 
close to the full-covering case and allowing Df a.s a. free parameter does not significantly 
change the properties of the sample. 

As noted earlier, this model is appropriate to a partial covering of the cylinder-of-sight 
and to a geometry in which some fraction of the underlying continuum, initially emitted 
in other directions, is scattered back into the line-of-sight. As discussed in §|^, a possible 
explanation for the spectral variability observed in NGC 4151 consists of clumps of 
photoionized gas (of differing column density) moving through the cylinder-of-sight on 
timescales <^ 1 day. Under such an hypothesis, the fact that Dj is similar for all the 
observations (at ~5%) suggests scattering as the most likely explanation. Circumstantial 
support is provided by the observation of a similar component in many type-2 AGN, 
with implications for unification schemes (e.g. Turner et al 1997a,b). Indeed the optical 
spectrum of NGC 4151 has been observed to exhibit both type-1 and type-2 characteristics 
at various times. It should be noted that the parameterization used here approximates 
the scattered continuum as a simple powerlaw as would be appropriate in the case of pure 
electron scattering in a very highly- ionized medium. In a more realistic treatment, the 
ionization state of the scattering gas would be taken into account, which under a wide range 
of parameter space would lead to emission and absorption features further complicating the 
spectrum ^ 2 keV. It should be noted that the fact that such a component is unambiguously 
observed only in NGC 4151 amongst our sample is simply due to the suppression of the 
transmitted continuum by the high column density of gas within the cylinder-of-sight of 
this source (with the gas having Ux such that there's still significant opacity) 



8.8. Implications of the '1 keV deficit' 

A number of sources have evidence for a deficit of counts at ~1 keV compared to the 
predictions of our photoionization models. As noted above, such a deficit occurs in those 
sources exhibiting the strongest absorption features. In § |6.2| we showed that this could be 
modeled with an additional edge (consistent with Ovill) in the rest frame of the source. 

There are several possible explanations of such a feature. First, it could be indicative 
that the cylinder-of-sight contains at least two absorption systems, each of different density, 
riH, and column density, Nh^^- Such an hypothesis is supported in the case of MCG-6-30-15 
by the short-timescale variability observed in the depth of the Ovill edge whilst the depth 
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of the Ovil edge remains constant (^). Circumstantial evidence for such a possibihty also 
comes from the Ginga observations that a large fraction of Seyfert 1 galaxies (including 
some of those considered here) contain an Fe ii'-shell absorption edge, far deeper than that 
predicted by the models considered here (Nandra & Pounds 1994; § |6.3.2| ). Unfortunately 
the geometry of the absorption systems cannot be constrained by the current data. One 
possible geometry is that both absorbing systems consist of complete screens, but at 
different radii from the central source. A full and correct treatment of such a scenario 
requires the ionization-equilibrium of the second screen of gas to be calculated in the 
same way as the first screen, but with the photoionizing continuum appropriate for that 
which has already passed through the first screen. Such a treatment is beyond the scope 
of the current paper. An alternative geometry is that the absorbing screen is clumpy 
on scale-sizes smaller than the emission region, resulting in the observed spectrum being 
the sum of the spectra transmitted by regions of different density and Nh,z- Yet another 
possibility is that there is indeed only a single absorbing cloud along the line-of-sight 
at any instant, but that different clouds (of different densities & Nh^z) move across the 
line-of-sight on timescales much shorter than the typical ASCA observation. In this case 
our analysis reveals time-averaged parameters of the absorber. Finally, we note that some 
of our detailed assumptions regarding atomic processes may be inapplicable to AGN, such 
as other processes important in determining the ionization equilibrium for Ovil/Ovill. 

Given the uncertainty in the form of the 'primary' continuum, the spectral resolution 
and typical signal-to-noise ratios typically achieved, it appears to be extremely difficult to 
distinguish between these possibilities using ASCA data. 



8.9. The location of the ionized gas 

The location and geometry of this ionized gas is currently unclear, although there are 
prospects for future progress. Of course if one knows the density, nj^, of the ionized gas, Ux 
and the intensity of the photoionizing source, one can use eqn. |l| to derive its radius as 

~ 3 X 105(i^)0-5 light - days (3) 
UxUh 

Assuming hh = 10^° cm~^ (as used for the model calculations), Ux — 0.1 and Lxaa, ~ 0.1-1 
(as appropriate for the majority of the objects presented here), one obtains rid ~few 
light-days. Such a location is comparable to that of the BELR in these objects. However 
as noted in §|^, our models are insensitive to densities in the range 10'^ ^ uh ^ 10^° cm~'^, 
thus values of ^few light-years cannot be excluded on this basis. 

As described in §^ the apparent variations in the column density of the ionized gas 
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seen in the X-ray band can be used to place upper limits on the radius of the gas if one 
assumes these variations are due to inhomogeneities in the matter passing through the 
cylinder-of-sight due to only Keplerian motion. For the few cases observed to-date, such an 
assumption also implies a location close to the BELR and hence is consistent with all/part 
of the same gas being responsible for the absorption features seen in the optical/UV. 
However such estimates are currently extremely crude, and it should be stressed that the 
bulk motion of the ionized gas is highly uncertain and could easily be dominated by other 
kinematics components. Constraints on the location could be placed by applying (under 
certain assumptions) photoionization models to both the UV/optical and X-ray absorption 
features. As will be described in §8.9.2| , the results from current studies are mixed and 



future progress requires simultaneous, high-resolution observations in all three wavebands. 
Reverberation mapping using the emission features produced within the ionized gas is 
another means whereby the location of the material could be determined. However as 
described in §fj.6|, current results are inconclusive. 



However we believe that in all likelihood there is ionized gas throughout the nuclear 
region in AGN. Multiple components, separated in at least velocity-space are seen explicitly 
in the optical/UV absorption features of some objects. At least two components are 
implied spectroscopically and/or from temporal variations in the X-ray absorption features 
of some objects. In addition, highly ionized gas (occupying a different region of Ux, Nh,z 
parameter-space than that responsible for the absorption on Seyfert 1 galaxies) is commonly 
invoked to explain the 'scattered' radiation in Seyfert 2 galaxies (e.g. see Turner et al 
1997b,c). 



8.9.1. The possible effect on UV-IR observations 

Should the volume of ionized gas contain embedded dust, and is located outside the 
appropriate emission regions, then one might expect a correlation between the column 
density Nh,z derived from X-ray observations and the various reddening indicators at longer 
wavebands. Such studies therefore have the potential for further constraining the location 
of absorbing material and/or geometry of the nuclear regions. Brandt, Fabian & Pounds 
(1996) have recently suggested that the infrared-bright quasar IRAS 13349+2438, which 
exhibits a large degree of reddening in the optical, might contain photoionized gas with 
internal dust. Reynolds (1997) noted that MCG-6-30-15, whose X-ray spectrum is clearly 
attenuated by ionized gas, also exhibits a large degree of reddening at longer wavelengths. 
The study of the creation, survival and implications of dust within the nuclear regions of 
AGN has received much attention for several decades (e.g. see Laor & Draine 1993 and 
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references therein). 

Unfortunately the use of the various reddening indicators in order to determine the 
column density of dust along the line-of-sight (and which can then be compared to that 
observed in the X-ray regime) rely upon making one or more assumptions. The primary 
assumption is of course that one knows the intrinsic value of the quantity from which 
the observed reddening is determined. If one knew the form of underlying continuum, 
multiwaveband continuum measurements could be applied in the optical-UV . However such 
a technique often suffers from a lack simultaneous observations (as well as observational 
'contamination' of other spectral components). Broad emission line ratios, in particular 
broad hydrogen lines, are poor reddening indicators since the lines are significantly affected 
by optical depth effects and there is no satisfactory theory to predict their intrinsic ratios 
(Netzer 1990). Clearly all techniques also require assumptions concerning the composition 
(chemical and grain-size distribution) of the dust itself, and well as the dust/gas mass 
ratio. With these points in mind, we consider only a single reddening indicator here 
(resonance-absorption line studies in the UV/optical are also considered in § ^.9.2| , but there 
do not provide a direct diagnostic of the dust-phase material). 

Fig. 18 shows (/i25//22o)obs (defined as the mean ratio of the observed flux at 125 nm 
to that at 220 nm, determined as described in § |2.1j and listed in Table 1) is shown against 
Nh^z for the datasets for which model C(ii) was considered acceptable (§5.6|). It can be seen 
that the majority of sources are consistent with 2 ^ (fub/ f 220) obs ^ 5, but that 3 datasets 
(NGC 3227, MCG-6-30-15(l,2)) have values a factor ~ 10 lower. Thus it is tempting to make 
the assumption that the intrinsic flux ratio0 is in the range 2 ^ {fi25/ f 220) int ^ 5, with the 
{fub/ f 220) obs lower in the case of NGC 3227 & MCG-6-30-15 due to reddening. Also shown 
in Fig. 18 is the predicted {f i2-> / f 220) obs assuming { f 12^ / f 220) int = 3.1 for various values of 
the difference between the optical depths at 125 nm and 220 nm (where Ar = ri25 — T22o)- 
It can be seen that NGC 3227 and MCG-6-30-15(l,2) are consistent with {fi2b/ f 220) int 
similar to the other objects if their fluxes at 125 nm and 220 nm are absorbed by a column 
density similar to that observed in the X-ray band and At = 0.5Nh^z/^^'^^ (dashed line). 
Interestingly, the 'standard Galactic' extinction curve (for a 1:1.12 silicate:graphite mixture 
of dust grains with sizes 5 < a < 250 nm, number density oc a~^'^ and total dust/gas mass 
ratio of 10~^ — see Laor & Draine 1993 and references therein) gives Ar close to this 
value (Ar = 0A2Nh,z/10^^)- Thus in the case of NGC 3227 and MCG-6-30-15(l,2), and 



^''it should be stressed that hne emission may contribute to /125 and (to a lesser extent) /220, and thus 
(/i25//22o)int IS not necessarily a good indication of the underlying UV contimmm. However, the use of the 
color here is only based on noting {f 125 / f 220) obs appears to be crudely constant (within a factor ^ 2) for the 
majority of the objects considered. 
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assuming a Galactic dust/gas mass ratio, this UV color is consistent with reddening by a 
column density of dusty material similar to that inferred by the X-ray observations. There 
are several implications of these results. 

First, we consider the implications on the location of the dust and gas in the case of 
NGC 3227 and MCG-6-30-15. The consistency between the Nh^^ derived from the X-ray 
absorption studies and the UV-reddening in the case of these two sources, assuming a 
'reasonable' dust/gas mass ratio, implies that the dust- and gas-phase material might occupy 
the same volume. If true, this allows us to place constraints on the location of this material 
based upon the survival of the dust. For a typical AGN continuum of bolometric luminosity 
Lhoi = fboixLxu (as defined in §|7.2D, and assuming maximum equilibrium temperature 



of ~ 1750 K, Laor & Draine find sublimation radii of Vsub — '^^ifboixLxu)^^'^ light-days 
and ~ l50{fboixLxuy^'^ light-days for graphite grains with a = 10"^ cm and 5 x 10^^ cm 
respectively. (Silcate grains, with a slightly lower maximum equilibrium temperature 
(~ 1400 K), have slightly larger values of Vsub-) Assuming fboix — 4, Vsub ~ 10 light-days 
and ~ 5 light-days for large graphite grains in MCG-6-30-15 {Lx44 — 0.3) and NGC 3227 
(Lx44 — 2 X 10~^) respectively. The corresponding value in the case of NGC 5548 
(Lx44 — 1) is Tsub ~ 40 light-days for large graphite grains. The latest results from the 
intensive reverberation mapping of the broad emission line region (BELR) in NGC 5548 
have suggested lags of ^ 2 ~14 days, with the lowest lags seen in the highest ionization 
lines (e.g. Korista et al 1995 and references therein). As first pointed out by Netzer & Laor 

1 /2 

(1993), since both the radius of the BELR and Vgub scale as Lf^'^i we conclude dust can only 
survive (i.e. in equilibrium) at radii just beyond the BELR in these objects. If the dust 
covers a substantial fraction of the BELR, then it will attenuate the broad emission lines 
in these objects. For instance, assuming the same 'standard Galactic' extinction curve used 
above, r(Clv) = 9Nh^z/^^'^'^ at the energy of Civ emission. Thus from the values of Nh,z 
found in MCG-6-30-15 and NGC 3227, a broad Cvi emission line is predicted to be totally 
absent in both sources, consistent with observations (see e.g. Courvoisier & Paltani 1992). 
The extinction in the optical is lower (e.g. r(H^) = 3Nh,z/10^'^, r(H„) = 2Nh,z/10'^^) 
but sufficient to attenuate ~80% and ~95% of the broad line, and ~ 65% and ~85% 
of the broad line in NGC 3227 and MCG-6-30-15 respectively (again assuming the 
dusty absorber completely covers the entire BELR). However it is clear from the value 



of rsub above and the discussion in § [7.3| that dusty clouds cannot survive at a sufficiently 
small radius that Keplerian motion alone can offer an explanation of the short-timescale 
variability (~ 10^ s) in the Ovill edge observed in MCG-6-30-15. Thus either the dusty 
absorber is extremely non-uniform and/or its kinematics are extremely non-Keplerian, or 
(far more likely) the gas giving rise to the bulk of the Ovill edge is at much smaller radii 
than that with embedded dust. Interestingly NGC 3227 and MCG-6-30-15 have the lowest 
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axial ratios (a/6, Table 1) in our sample, with the exception of the well-known, edge-on 
galaxy IC 4329A. 

We now consider the majority of the objects on our sample, for which {fi2b I f 220) obs 
provides no evidence for reddening. There are two obvious explanations. First, the column 
density of ionized gas within the cylinder-of-sight may simply not contain any embedded 
dust in these objects. This might arise from differences in the origin of the material, or as 
a result of the material in these objects being within Tgub- Alternatively the dust may be 
present yet unobservable. This will true if the dust is of a different composition in which 
Ar ~ (such as is the case if the grains are large), or if the dust clouds cover only a 
fraction of the UV emission region. 

8.9.2. Relationship to 'associated absorbers' in the UV/ optical 

The presence of highly-ionized gas, instrinsic to the Seyfert 1 nuclei is also obtained 
from observations carried out in the UV. lUE observations revealed only a small fraction 
(~3%) of sources to exhibit narrow absorption lines close to the systemic velocity (e.g. 
Ulrich 1988). However, recent results taking advantage of the higher spectral resolution 
and signal-to- noise ratio available with the FOS and GHRS instruments onboard HSTheive 
revealed that up to 50% of Seyfert 1 galaxies may in fact contain such absorption features 
(e.g. Crenshaw 1997). Such features have also been reported in ifC/T observations of 
NGC 4151 (Kriss et al. 1992) and NGC 3516 (Kriss et al. 1996b). Furthermore variability 
in the UV absorption features have been observed on timescales of weeks-years in NGC 4151 
(Bromage et al. 1985), NGC 3783 (Maran et al 1996), and NGC 3516 (Koratkar et al. 
1996), and multiple radial velocity components have been reported in Mrk 509 (Crenshaw, 
Boggess & Wu 1995) and NGC 3516 (Crenshaw, Maran & Mushotzky 1997). 

As noted by Crenshaw (1997), there is a reassuring trend whereby Seyfert 1 galaxies 
which have evidence for strong absorption features due to ionized gas in the X-ray band 
also tend to have strong absorption lines in the UV/optical. This raises the question as to 
whether the associated absorbers seen the UV/optical are caused by (part or all of) the 
same ionized material which gives rise to the absorption features observed in the X-ray 
band. The UV absorption components are often blueshifted with respect to the emission 
hues, but with velocity shifts (typically ~fewxlO^ km s~-^) too small to be detectable with 
current X-ray observations. Thus direct tests for kinematic consistency are not yet possible. 
However various attempts have been made to relate the UV/optical and X-ray absorption 
features using photoionization calculations (e.g. Mathur 1994; Mathur et al. 1994; Mathur, 
Wilkes & Elvis 1995; Shields & Hamann 1996; Kriss et al 1996a,b). These have met 
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with varying degrees of success, which is hardly surprizing given the uncertainties in the 
location and kinematics of the gas and the spectral form of the photoionizing continuum. 
Nevertheless there have been a number of impressive successes linking the UV/optical and 
X-ray absorbers (e.g. see Mathur 1997). The observed UV/optical absorption features 
are often superimposed on the wings of broad emission lines. This clearly places the 
gas responsible for these 'associated absorbers' at radii greater than that where the bulk 
(of the blue-wing) of that emission line originates (^ 2 ~14 light-days in the case of 
NGC 5548). However ionized gas could exist at a wide range of radii, with only part of the 
gas responsible for the X-ray absorption within the cylinder-of-sight to the UV/optical 
emission line region(s). Indeed the presence of multiple components to the UV/optical 
absorption features (and the indication of additional absorption zones in the X-ray — e.g. 
§ |8.(j| ) clearly indicates a single screen of gas to be an over-simplification. Thus the column 
densities implied by the UV/optical absorption features may be small (and dominated by 
different ionization states) compared to those implied from the X-ray features. Furthermore, 
the emission and absorption features predicted in the UV/optical from the ionized gas are 
somewhat dependent on the form of the photoionizing continuum in the UV-X-ray band. 
This is poorly determined for all objects, leading to further ambiguities. In addition, in at 
least some objects the absorption features in both the UV and X-ray are known to vary 
with time. To date, most comparisons between the absorption features have been performed 
using non-simultaneous data, leading to obvious risk of deception. Finally, there is the 
usual irritation of uncertainties in the abundances within the absorbing material. The best 
hope for progress is provided by intense, simultaneous monitoring of the UV/optical and 
X-ray absorption features. 



9. SUMMARY & OPEN ISSUES 

We have presented the systematic analysis of the 0.6-10 keV ASCA spectra of 18 
Seyfert 1 galaxies in order to examine the characteristics of X-ray absorption by ionized gas. 



9.1. Summary of Results 

In summary, we find: 

• For the majority of sources, an underlying powerlaw with F ~ 2 dominates the 
emission within the ASCA bandpass (§ p.2| and Fig. 2). 

• Absorption by neutral material does not provide a satisfactory model for the X-ray 
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attenuation for most objects (§ |8.3|) . However, there is evidence for a significant 
column density of ionized gas in the hne-of-sight to 13/18 objects (§ ^.4| ). 



The X-ray ionization parameter for the ionized material is strongly peaked at Ux ~ 0.1 
(Fig. 7), which may in part be the result of selection effects. The column densities of 
ionized material are typically in the range N^^z ~ 10^^-10^^cm~^ (Fig. 3), although 
highly ionized (and hence psuedo-transparent) column densities up to 10^^cm~^ 
cannot be excluded in some cases (see also 5B.5I) 



Allowing a fraction of the continuum to be observed without attenuation significantly 
improves the fit to many sources ( §8.7| and Fig. 5). Such a component is required only 
in the case of NGC 4151. 

Inclusion of the emission from the ionized material significantly improves the fit to 
many sources, at an intensity consistent with the material subtending ~ Att sr at the 
central source. However such a component is not required in any individual source 
(Qand Fig. 11). 

A deficit of counts is observed at ~ 1 keV in the sources exhibiting the strongest 
absorption features (§ p.8| ). This is consistent with the presence of a second zone of 
(more highly) ionized gas, which might have been seen previously in the deep Fe 
i^-shell edges observed in some Ginga observations. 

There is evidence that the ionized material in NGC 3227 and MCG-6-30-15 contains 
embedded dust, whilst no such evidence in the other sources (§ ^.9.1| and Fig. 18). 



9.2. Open issues 

While X-ray absorption by ionized gas appears to play a very important role in shaping 
the observed X-ray spectra of AGN, there remain a large number of important questions 
which cannot be answered with the available data, including: 

• What is the frequency of occurrence of warm absorbers at higher luminosities. Is 
there a simple relationship between Ux and Lx ? 

• Are there several distinct warm absorber components attenuating the X-ray bandpass 
? 
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• What is the location, density and distribution of the material comprising the X-ray 
warm absorber ? Sensitive mcasTircments of the emission spectrum from the gas will 
help to constrain the gas distribution. 

• Is the warm absorber flowing into or out from the active nucleus ? Do the UV 
associated-absorbers arise from the same X-ray absorbing gas ? 

• Does the gas to dust ratio differ significantly between objects, and how does this 
affect the observed properties of the source ? How important is this component in the 
unification of AGN ? 

• Is the absorption variability due to bulk motion of material across the line-of-sight, 
variable ionization or a combination of both effects ? 

• How stable is the warm absorbing gas ? 
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APPENDIX: NOTES ON INDIVIDUAL OBJECTS 

The sources presented here are all well-known, bright Seyfert 1 galaxies, detected early 
in the history of X-ray astronomy. Thus, most have been observed by all the major X-ray 
astronomy satellites. The results from, and references to, almost all X-ray observations 
published prior to the end of 1992 are included in the compilation of Malaguti, Bassini & 
Caroli (1994). The following abbreviations are used for papers cited numerous times in 
this section: NP94 - Nandra & Pounds (1994); R97 - Reynolds (1997); TP89 - Turner & 
Pounds (1989); T91 - Turner et al (1991); W95 - Weaver, Arnaud & Mushotzky (1995). As 
stated in §8.4| , R97 also fitted a single-zone photoionization model to a number of sources, 
and quotes the ionization parameter ^ (= L/hhR^-, where hh and R are as defined in 
eqn.|l], but where L is the integrated luminosity) calculated over the 0.0136-13.6 keV band 
assuming a single powerlaw. For convience, below we give the corresponding t/^^'' (defined 
as in eqn.[^ for such a continuum over the 0.1-10.0 keV band) using the conversion factors 
shown in Fig. lb. Ginga observations have shown evidence for a significant Fe absorption 
edge in the 7.1-8.9 keV band for a number of sources (NP94), in some cases confirming 
suggestions in earlier observations. These sources are noted, but as described in § |6.3.2| a 
detailed analysis of such features has not been attempted here. Unless otherwise stated, 
the Seyfert classification is from Whittle (1992) and all limits on a parameter are quoted at 
90% confidence for 1 interesting parameter. 

As a reminder, models A(i)-C(ii) all assume a single powerlaw continuum absorbed 
by a neutral column density at z = (constrained to be > A^"q, the Galactic Hi column 
density along that line-of- sight) and an additional column density, Nh^zi at the redshift 
of the source. In models A (§|5.1|-pl2|) Nh^^ is assumed to be neutral, whilst in models 
B (§p3|- p^ & C (§5.5|-|5^) Nh,z is assumed to be photoionized, with a corresponding 



ionization parameter Ux- Models C also include the emission features expected from the 
ionized gas (assuming it subtends a solid angle VL at the central source). In the models 
denoted (i), the absorbing gas covers the entire cylinder-of-sight to the source, whilst in 
models (ii) a. fraction -D/ of the underlying continuum are observed without attenuation. 
See for further details. 

Mrk 335: 

This Seyfert 1.0 galaxy was first detected in X-rays by UHURU (Tananbaum et al 1978). 
From our analysis of the data from ASCA observation performed in 1993 December, we find 
model A(i) to provide an unsatisfactory description of the spectra (Table 3). Fits satisfying 
our criteria for acceptability are obtained if ~60% of the continuum escapes without 
suffering attenuation by neutral material (model A(ii)) and for all the models in which the 



absorbing gas is assumed to be ionized (models B(i)-C(ii)). However, in § |6.4.2| we obtained 
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a superior fit if tlie continuum is assumed to be tlie sum of a steep powerlaw (r^ = 3.3I0!?) 
dominating at energies ^ 1 keV and a flatter powerlaw (F^ ~ LSlg g) dominating at energies 
^ 3 keV. Tliis model gave x^/dof = 621/642, had Rg/h ~ 1-6 and required no absorption in 
addition to A'^^^q. Such a double-powerlaw continuum is in good agreement with previous 
observations (Turner et al 1993b and references therein). The flux in the 2-10 keV band 
during the ASCA observation was F2-10 = 9.5 x 10~^^ erg cm~^ s~^, similar to that during 
the BBXRT era, and indeed the X-ray spectrum of Mrk 335 appears relatively stable on 
long timescales, though the source exhibits variability on short timescales (Pounds et al 
1987; Lee et al 1988; Paper I). In an independent analysis of the same ASCA dataset, R97 
modeled the continuum in terms of single powerlaw and a blackbody (with kT ~ 0.14 keV 
and luminosity Lbb ~ 10"^^ erg s~^) to account for the steepening to low energies. We 
confirm these results, but find such a model provides an inferior description of the data 
{x^/dof = 629/641) compared to the double-powerlaw model. Explanations for the X-ray 
continuum in Mrk 335 have been presented which involve e^-pair cascades (Turner et al 
1993b) and the Comptonization of the spectrum emitted by an accretion disk (Zheng et al 
1995a). 

We find any features imprinted by any ionized gas to be relatively weak during the 
epoch of the ASCA observation, resulting in poor constraints on Nh,z and Ux in models 
B(i)-C(ii). Turner, George & Mushotzky (1993) found evidence for an absorption edge 
(r ~ 0.3 ± 0.1) in the 0.7-0.9 keV band during a ROSAT PSPC observation of Mrk 335. 
Furthermore, R97 found the addition of an Ovil edge (of optical depth tqj ~ 0.3) 
significantly improved the fit to a single powerlaw model (F ~ 2.1) to these ASCA 
data. However, we find only a marginal improvement (Ax^ = 4) when such an edge 
{to7 = 0.2 ± 0.2) is added to our double-powerlaw model, and no improvement when an 
Ovill edge is also included {tqs ^ 0.1). Some evidence exists for a Fe i^'-shell edge in 
the 7.1-8.9 keV band in two Cinga observations (NP94), and possibly in an EXOSAT 
observation (Turner et al 1993b). The energy of the features could not be constrained in 
any of the datasets, but NP94 inferred Nh,z ~ 2 x 10^^ cm~^. 

There is no evidence for any UV absorption features due to ionized gas at the redshift 
of Mrk 335 (Zheng et al 1995a; Crenshaw 1997). 

Fair all 9: 

This Seyfert 1.0 galaxy was suggested to be an X-ray source based upon its location in 
UHURU and ylne/-F error boxes (Ricker 1978). From our analysis of the data from ASCA 
observation performed in 1993 November, we find no model in §^ which satisfies our formal 
criteria for acceptability. However, this is primarily due to scatter in 2-4 keV band, most 
likely due to an underestimation of systematic errors for this dataset. In view of this. 
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we consider model A(i) to provide an adequate description of the data. However model 
A (a) provides a significant improvement in the goodness-of-fit if ~ 80% of the continuum 
escapes without attenuation by a large column density {Nh.z ~ 3 x 10^'^ cm~^) of neutral 
material which gives rise to a relatively subtle upturn in the spectrum ^ 5 keV (Fig 4). 
No significant improvement in the goodness-of-fit was obtained assuming the more complex 
models described in or §^ For all models we find no evidence for absorption in addition 
to Nf^Q and a derived slope (F ~ 2.0) consistent with the Einstein (T91) and EXOSAT 
(TP89) observations of the source, slightly steeper than that observed by Ginga (F = 1.84, 
NP94) and only marginally consistent with the ROSAT PSPC observations (F = 2.2 ± 0.2, 
Walter & Fink 1993). X-ray variability has been reported on both short and long timescales 
(Morini et al 1986; Paper II). Unfortunately these ASCA do not offer any insight into the 
possible spectral variability suggested by the EXOS'tIT observations (Morini et al 1986). 
However, we note that the flux in the 2-10 keV band observed during both the Ginga 
(NP94) and this ASGA observation (F2_io = 2 x 10~^^ erg cm~^ s~^) indicates that the 
historical decline in the brightness of Fairall 9 in the X-ray band (Morini et al 1986) has 
ceased. 

We find any features imprinted by any ionized gas to be weak during the epoch of the 
ASGA observation, resulting in poor constraints on Nh,z and Ux in models B(i)-G(ii). R97 
also found a lack of evidence for Ovil and Ovill absorption edges in their analysis of the 
ASGA data from Fairall 9. They find a best-fitting photon index consistent with our value, 
although unlike our findings, require a small column density {Nu^^ ~ 10^° cm~^) of intrinsic 
absorption. 

There is no evidence for any UV absorption features due to ionized gas at the redshift 
of Fairall 9 (Zheng et al 1995b; Crenshaw 1997). 

3C 120: 

This bright, flat-spectrum radio source, has Seyfert 1 properties in the optical regime. The 
source was first detected in X-rays by SAS-3 (Schnopper et al 1977), and has since been 
observed by all the major X-ray astronomy satellites except Ginga. From our analysis of 
the data from the ASGA observation performed in 1994 February, we find model A(i) to 
provide an adequate description of the spectra, with some evidence for intrinsic absorption 
{Nh,z ~ 5 X 10^° cm~^). However, we find some evidence for a 'hard tail' in this source, 
with significantly superior fits obtained if ~ 30% of the continuum is attenuated by 
Nh,z ~ 5 X 10^^ cm~^ (model A(ii); Fig. 4), or by a flattening of the underlying continuum 
^ 6 keV (§ |6.4.2| ). We find no compelling evidence for ionized gas in 3C 120 (intense 
emission from ionized gas is implied by the best-fitting solution assuming model G(ii), but 
at a level incompatible with strength of the observed Fe iiT-shell line - see Fig 12). This 
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confirms the finding of R97 wlio also found a lack of evidence for Ovil and Ovill absorption 
edges in their analysis of these ASCA data. They find a best-fitting photon index consistent 
with our value, but a slightly larger absorbing column density instrinsic to the source 
{Nh,z ~ 8 X 10^° cm"^) then we found for model A(i). 

In all cases, the derived slope is F ^ 2.0, and is consistent with a poorly constrained 
iiOS'^r PSPC spectrum (T ~ 1.5l|'g, Boiler et al 1992), and slightly steeper than than 
those derived from HEAO-1, Einstein and EXOSAT {V ~ 1.5-1.9, W95; T91; TP89). 
Spectral variability was apparent in the Einstein observations (T91). Unfortunately the 
ASCA data are unable to constrain the Compton-Refiection component seen in HEAO-1 
observations (W95), although such a component clearly offers an an alternative explanation 
of the hard tail. We also note that many of the fits to the 3C 120 dataset imply a 
local column density > iV^°Q {^Nh^ — 5 x 10^° cm~^). Furthermore this dataset seems 
particularly badly affected by a trough-like deficit of counts at energies < 0.6 keV compared 
to the extrapolated model. It is unclear how these facts may be related to the well-known, 
but poorly-understood calibration problems of the XRT/SIS at these energies (see also 

NGC 3227: 

This Seyfert 1.2 galaxy was first detected in X-rays by Ariel-V (Elvis et al 1978). In our 
analysis of the data from ASCA observation performed in 1993 May, we find model A(i) to 
provide an unsatisfactory description of the spectra (Table 3). Fits satisfying our criteria for 
acceptability are obtained if ~35% of the continuum escapes without suffering attenuation 
by neutral material (model A(ii)] Fig. 4). However yet superior fits are obtained if the 
absorbing gas is assumed to be ionized (model B(i)] Fig. 6). The parameters of the ionized 
gas are fairly well constrained, with Ux — 10~^ and Nh,z ~ 4 x 10^^ cm~^ consistent with 
the findings of Ptak et al (1994) and R97 (f/^^^ ^ 0.02, A^^^,^ ~ 4 x lO^^ cm'^) in their 
independent analysis of this ASCA observation. 

We find that statistically there is no requirement for any of the underlying continuum 
to escape without suffering attenuation by the ionized material, or for any significant 
emission from the ionized gas, (models B(ii)-C(ii)). In all cases there is no requirement for 
additional absorption by neutral gas in excess of Nj^^. The derived photon index (F ~ 1.6) 
found for NGC 3227, is flatter than for the majority of the sources (Fig. 2), consistent 
with previous measurements by HEAO-1 (W95), Einstein (T91), EXOSAT (TP89), Cinga 
(NP94) and ROSAT (Rush et al 1996). However, it is important to note that EXOSAT 
observations of NGC 3227 revealed significant spectral variability on timescale as short as 
10^ s, most likely as a result of changes in absorption (TP89). Furthermore, such variability 
was also present during the ASCA observation presented here (Ptak et al 1994; Paper I). 
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Thus it should be remembered that our results most likely represent the time-averaged 
values. We note evidence was found for a Fe i^-shell edge in the 7.1-8.9 keV in two 
Ginga observations of the source. NP94 found evidence for absorption by ionized gas with 
Nh.z ~ 3 X 10^^ cm~^, higher than that observed here. 

NGC 3516: 

This source, one of the original galaxies classified by Seyfert (1943) as a type 1.0 galaxy, 
was first detected in X-rays by Einstein (Maccacaro, Garilli & Mereghetti 1987). In our 
analysis of the data from an ASCA observation performed in 1994 April, we find no 
model in ^ which satisfies our formal criteria for acceptability. Nevertheless, we find 
model B(i) to provide a significant improvement over models A(i) & A(ii), offering clear 
evidence for ionized gas and confirming the results of R97 and those from a second ASCA 
observation carried out in 1995 March (Kriss et al 1996a). However, whilst a single-zone 
photoionization model clearly reveals a strong '~1 keV deficit in this source, we find that 
despite providing a vast improvement, a statisfactory fit is still not obtained with the 
inclusion of an additional absorption edge (§ |6.2| ). Interestingly, we note that neither R97 nor 
Kriss et al (1996a) were able to find a solution which satisfy our criteria for acceptability. 
We find a derived index (F ~ 1.9), steeper than that which has been claimed from previous 
Einstein (Kruper, Canizares & Urry 1990), EXOSAT (Ghosh & Soundararajaperumal 
1991) and Ginga (Kolman et al 1993) observations in the 2-10 keV band. However all 
these observations required a neutral Nh^ » Nf^^, undoubtedly due to the incorrect 
modelling of the deep Ovil and Ovill edges in this source with neutral gas. A steeper index 
(F ^ 2.4) is implied when a single powerlaw model is applied to the ROSAT PSPC data 
(Boiler et al 1992), again probably the result of the deep Ovil and Ovill edges. For models 
B(i)-G(ii), we find a neutral Nha ~ 5 x 10^° cm"^ and an ionized Nh,z ~ 9 x 10^^ cm~^, 
with Ux ~ 0.1. For comparison, R97 find Nh,z ~ 10^^ cm'^, F ~ 1.7 and Uf^ ~ 0.06 
from an independent analysis of the ASGA dataset obtained in 1994 April, and assuming 
a single-zone model. Kriss et al (1996b) have performed a detailed analysis of the 1995 
March observation (using XSTAR), and for a single- zone model find Ux ~ 0.03 (assuming 
their quoted spectrum), N^^z ~ 10^^ cm~^ and F ~ 2.0. These workers also find adding a 
second zone of ionized gas significantly improves the fit, obtaining a best-fitting solution 
with Ux ~ 0.1, Nh,z ~ 1.4 X 10^^ ^,^^-2 _ q_q2^ jy^^^ _ 7 x iQ^i cm'^ for the two 

zones. Nandra et al (1997c) have recently shown the Fe emission line varied significantly 
between the two ASGA observations. Evidence for absorption by ionized gas was found in 
3 Ginga observations of NGC 3516, with Nh,z ~ 2-3 x 10^^ cm'^ (Kolman et al 1993). 
We also note that ROSAT YiYiY observations show possible evidence for elongation of some 
fraction of the soft X-ray emission in NGC 3516 in a direction consistent with the extended 
narrow- line region (Morse et al 1995). 
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NGC 3516 has long been known to exhibit deep, variable absorption lines in the UV 
(Koratkar et al. 1996; Kriss et al 1996a and references therein) with at least 4 distinct radial 
components visible in Civ (Crenshaw, Maran & Mushotzky 1997). From the combined 
analysis of the simultaneous ASCA and HUT, Kriss et al (1996b) conclude that a range 
of column densities and ionization parameters are required to explain the depths of all the 
X-ray and UV absorption features seen in NGC 3516. 

NGC 3783: 

This Seyfert 1.2 galaxy was first detected in X-rays by Ariel-V (Cooke et al 1976). 
Absorption due to ionized gas in NGC 3783 was first suggested by ROSAT PSPC 
observations (Turner et al 1993). In our analysis of the data from 2 ASCA observations 
separated by 4 days in 1993 December, we find models involving absorption by neutral 
material (model A(i) A(ii)) to provide unsatisfactory descriptions of the spectra. However 
fits satisfying our criteria for acceptability are obtained if the absorbing gas is assumed to 
be ionized (models B(i) & B(ii)), with a significant improvement in the goodness-of-fit with 
the inclusion of the ionized emitter (models C(i) h C(ii)). In all cases, the parameters of the 
ionized gas are fairly well constrained, with Ux ~ 0.1-0.15 and Nh,z ~ 2 x 10^^ cm~^ (e.g. 
Fig. 16). This confirms the findings of George, Turner & Netzer (1995), and R97 who found 
the addition of Ovil and Ovill edges (of optical depths tqj ~ 1.2 & tqs ~ 1-4) significantly 
improved the fit to a single powerlaw model (F ~ 1.4) to these ASCA data. R97 also fitted 
a single-zone photoionization model and found f/^^'' ~ 0.08, Nh,z ~ 2 x 10^^ cm~^, F ~ 1.7. 
The index derived for these ASCA observations (~ 1.8) is somewhat steeper than that 
found by George et al (1995), primarily due to these workers not adequately modeling the 
Fe i^-band. A flatter spectral index has also been suggested by some (but not all) previous 
observations in the 2-10 keV band (e.g. Ghosh et al 1992). This may indicate spectral 
variability or (more likely) be an artifact of incorrectly modeling the deep absorption 
troughs due to the ionized gas in this source. As described in § [7.1| , the behaviour of the 
ionized gas in response to the change in flux seen between the two ASCA observations 
is consistent with expectations. NP94 found evidence for absorption by ionized gas with 
Nh,z ~ 10^^ cm~^, higher than that observed here. 

UV observations of NGC 3783 have revealed absorption due to Lya, Nv and Civ, with 
the latter exhibiting variability on a timescale of ^ 6 months (Reichert et al 1994; Lu et 
al 1994; Maran et al 1996). Shields & Hamann (1997) have compared the implied column 
densities with the photoionization modeling of the X-ray absorption presented in George 
et al (1995). Whilst alternative scenarios could not be excluded, they found the UV and 
X-ray results to be consistent with a single-phase, photoionized plasma. By comparing 
the observed variability timescale in Civ with the radiative recombination, they derived 
a density n ;^ 50 cm~^ for the ionized gas and hence, from the ionization parameter, a 
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location r ^ 30 pc from the ionizing source. Since the parameters associated with the X-ray 
absorption are are generally consistent with those in George et al (1995), we agree with the 
conclusions of Shields & Hamann (1997). 

NGC 4051: 

This Seyfert 1.5 galaxy was first detected in X-rays by Einstein (Marshall et al 1983) and 
has been seen to exhibit rapid variable on timescales as short as fewxlO^ s (e.g. Lawrence 
et al 1987). Indeed, in Paper I we showed NGC 4051 the source to exhibit the most 
pronounced variability of any of the sources considered. Absorption due to ionized gas in 
this object was first suggested by ROSAT PSPC observations (Pounds et al 1994). In our 
analysis of the data from the ASCA observation performed in 1993 April, we find no model 
in §^ which satisfies our formal criteria for acceptability. A satisfactory fit is found however, 
if a strong (JF ~ 4) Compton-refiection component is included in the model (§ |6.4.1D . 

However, such a curvature in the observed spectrum may be due in part to the spectral 
variability observed during the observation (Paper I). Nevertheless, we do find a significant 
improvement in the fits for models including ionized gas, with Ux ~ 0.7, Nh,z ~ 2- 
10 X 10^2 cm-^ F ~ 2.3 for all models B(i)-C(ii), and Ux ~ 0.4, Nh,z ~ 10^^ cm'^, F ~ 2.4 
for the model including Compton-refiection (Table 10). Mihara et al (1994) and R97 have 
both performed independent analyses of this ASCA dataset. Both find the addition of 
two edges to a single powerlaw model significantly improved the goodness-of-fit over the 
0.6-1.0 keV band, confirming the results from the ROSAT PSPC and from a later ASCA 
observation carried out in 1994 July (Guainazzi et al 1996). R97 fixed the edge energies to 
those appropriate for Ovil and Ovill and found tqi ~ ros ~ 0.2. 

Photoionization models were also applied to these ASCA data by Mihara et al and 
R97, as well as to the earlier ROSAT PSPC data (Pounds et al 1994; McHardy et al 1995) 
and data from the later ASCA observation (Guainazzi et al 1996). Unfortunately detailed 
comparisons are difficult due to differing assumptions regarding the form of the observed 
and photoionized spectra. Mihara et al find ^ ~ 50 (but for an unspecified XUV continuum) 
and Nh^z ~ 2 x 10^^ cm^^ when the observed continuum is modeled by a (relately) fiat 
powerlaw (F ~ 1.9) plus a low temperature blackbody {kT ~ 0.1 keV and luminosity 

~ 3 X 10^1 erg s-^). R97 find Uf' ~ 0.02, Nh,z ~ 10^^ cm'^ and F ~ 1.9, but also find 
evidence for a similar blackbody component. From the analysis of combined ROSAT-Ginga 
data Pounds et al (1994) also prefer the inclusion of a blackbody [kT ~ 0.2 keV and 

~ 5 X 10"^^ erg s~^) and powerlaw (F ~ 2.1) finding Nh,z ~ 7 x 10^^ cm"^ and 
U ~ 0.8 (but for an unspecified XUV continuum). From subsequent ROSAT data 
McHardy et al (1995) find Ux ~ 0.5 (given their continuum), Nh,z ~ 8 x 10^^ cm"^ 
and F ~ 2.2, with no requirement for a cool thermal component. Finally, Guainazzi et 
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al find Nh,z ~ 2 x 10^^ cm~^, T ~ 2.1 and Ux ~ 0.2 (given their continuum). These 
workers found evidence for additional emission ^ 1 keV, which could be modeled by 
either a blackbody or a series of emission lines. However given the differential variability 
between the depths of the Ovil and Ovill edges seen on a timescale ~ 10^ s during the 
1994 observation, Guainazzi et al considered neither explanation to be satisfactory. We 
find no strong requirement for a low temperature blackbody component in our analysis 
of NGC 4051 (with our models for having 1.0 ^ i?o.6 ^ 1-2) • For our model including 
Compton- reflection we find ^ 3 x 10^^ erg s~^ for /cT = 0.1 keV. Instead our findings 
do imply significant emission from the photoionized gas in NGC 4051 at energies ^ 1 keV 
(with a luminosity of ~ 10^^ erg s~^ in the 0.1-1 keV band), emission which was ascribed 
to the blackbody component by previous workers. However we do note that the derived 
intensity of this emission is up to a factor 3 times that expected given the intensity of the 
continuum during the observations. It is currently unclear whether these results indicate 
that the source was brighter in the past, that there is indeed an additional component at 
the softest energies, or the physical conditions with NGC 4051 are much more complex. We 
consider it unlikely that significant progress in decouphng the complex spectral variability 
exhibited by NGC 4051 ^ 1 keV will not be made until high signal-to-noise observations 
are performed with high spectral and temporal resolution. Nevertheless, the time-averaged 
results presented here are useful for comparison with those from those sources which appear 
to exhibit less dramatic variability. 

NP94 found evidence for absorption by ionized gas in Ginga observations, with 
Nh,z ~ 6 X 10^2 cm-2. 

NGC 4151: 

This Seyfert 1.5 galaxy was first detected in X-rays by UHURU (Giacconi et al 1974), and 
has been extensively studied ever since at at all wavebands. Interestingly, despite being one 
of the brightest Seyfert 1 galaxies in the X-ray band, the source appears to be somewhat 
atypical (e.g see Warwick et al 1996 and references therein). Specifically, it has long 
been known that NGC 4151 appears to have a flat powerlaw (F ~ 1.3-1.7), a substantial 
column density of absorbing gas instrinsic to the source {Nh,z ~ 10^^-10^^ cm^^) and a 
strong 'soft-excess' at energies <^ 2 keV. At present, the consensus is that the usually hard 
underlying continuum is the result of different physical conditions associated with the 
emission mechanism compared to most other objects in its class (but see Zdziarski, Johnson 
& Magdziarz 1996). The soft-excess has been interpreted as either due to separate emission 
component (e.g. Warwick, Done & Smith 1995) and/or the result of the absorbing gas only 
partially covering the source (e.g. Holt et al 1980). The suggestion that the absorbing gas 
in NGC 4151 might be photoionized was first made by Yaqoob, Warwick & Pounds (1989). 
However subsequent observations, in particular the different variability characteristics 
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observed below and above ~ 1 keV, have suggested the situation may be more complex (e.g. 
Weaver et al 1994a,b; Warwick et al 1996). Emstein HRI (Elvis, Briel & Henry 1983) and 
ROSAT HRI (Morse et al 1995) observations have shown spatially resolved X-ray emission, 
probably associated with the extended narrow-line region, at an intensity likely to also 
contribute to the observed spectrum ^ 1 keV. 

In our analysis of the data from 6 ASCA observations carried out prior to 1994 
May, we found 3 datasets fulfilling the criteria defining this 'sample' (see also Paper I). 
As summarized in Table 12, we find that a fit satisfying our criteria for acceptability is 
achieved for the first two observations (NGC 4151(2,4)) only when the absorbing material 
is allowed to be ionized and ~5% of the underlying continuum is allowed to escape without 
suffering any attenuation (i.e. models B(u) & C(ii)). Such an hypothesis is strongly 
supported by the third observation (NGC 4151(5)), even though none of the models in §^ 
provide a solution which formally satisfies our criteria for acceptability (see Fig. 9). In 
§ |6.4.3| we found evidence in all 3 datasets for additional emission due to an optically-thin 
plasma at a temperature and intensity consistent with the extended emission. In all 
cases we find best-fitting values of F ~ 1.5, Ux ~ 0.07, Nh,z ~ 4-7 x 10^^ cm~^, and 
Lx ~ 1.4-2.0 X 10"^^ erg s~^. However, as described in §7.2| , the ionized gas appears to 
respond in an interesting way in response to changes in the intensity of the illuminating 
continuum between observations. A possible explanation of such behaviour was proposed 
in which the gas responsible for the absorption is clumped into regions of differing N^^^ 
with a characteristic scale-size greater than that for the central source. It was shown that 
Keplarian motion would allow such Nh,z perturbations to traverse the cylinder-of-sight 
on a timescale consistent with observations for reasonable range of radii, giving rise to an 
apparent decoupling between variations in Ux-, Nh,z and the intensity of the underlying 
source. 

These NGC 4151 datasets have been independently analysed by other workers. 
Weaver et al (1994b) considered the 4151(2) dataset and found a satisfactory fits with 
F ~ 1.6, U ~ 0.3 (for an unspecified XUV spectrum), Nh,z ~ 5 x 10^^ cm'^, Df ~ 0.03. 
However when a similar model was applied to the 4151(1) dataset, these workers found 
U (and Nh,z) to be higher when the source was fainter. (We confirm this general 
behaviour, finding F ~ 1.4, Ux ~ 0.13, Nh,z ~ 10^^ cm^^ Dj ~ 0.09, n/Air ~ 4 and 
Lx ~ 0.9 X 10^^ erg s~^ when model C(ii) is applied to the 4151(1) dataset.) Weaver et al 
also suggested an explanation is terms of inhomgeneities in the ionized material trasversing 
the cylinder-of-sight. 

Zdziarski, Johnson & Magdziarz (1996) have considered the 4151(1) ASCA dataset 
along with contemporaneous CGRO OSSE data. They find the spectra ;^ 3 keV to be 
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better described by thermal Comptonization, and dual, cold absorber (with twice cosmic 
abundance of Fe), rather than by an ionized absorber. Furthermore, they show that if 
a Compton-reflection component is added, the underlying spectral index of this source 
is steepened to F ~ 1.8 and hence similar to that observed in other Seyfert 1 galaxies. 
However, since we have not included 4151(1) in our general analysis we feel unable to 
comment further on this result here (although see above, and also note that in § |6.4.1| we did 
not find a significant improvement when Compton-reflection was included in our analysis of 
the NGC 4151(2,4,5) datasets). Nevertheless, we note Zdziarski, Johnson & Magdziarz also 
consider several archival X- and 7-ray datasets, and make a a number of points concerning 
the emission mechanism in this source. 

Warwick et al (1996) considered the 4151(4,5) datasets (as well as 4151(3,6), and 
psuedo-simultaneous ROSAT PSPC and CGRO OSSE observations). They confirmed 
previous findings that the time-variability patterns differ in the soft and hard X-ray 
bands, with the ROSAT PSPC count rate in the 0.1-1 keV band remaining constant at 
~ 0.46 count s~^ whilst the 1-2 keV count rate increased a factor ~2 (to ~ 0.32 count s~^) 
on a timescale ~ 2 days (around the time of the NGC 4151(3) ASCA observation). 
Unfortunately the i?0S'j4r observations ended ~ 2 days later (and at the time of the 4151(4) 
ASCA observation), but the subsequent ASCA observations showed that 1-10 keV count 
rate remained relatively constant (^ 30% variability) for a period of ~ 6 days thereafter. 
We do not attempt a detailed comparison of our results with these findings here. However 
we note that our best-fitting model to the 4151(4) dataset does predict a count rate in the 
1.0-2.0 keV band of the ROSAT V^VC of 0.30 count s~^, consistent with that observed. 
The corresponding count rate prediction for the 0.1-1.0 keV is clearly a strong function 
of the spectrum assumed < 0.6 keV especially the value assumed for Nu^. This band is 
clearly poorly constrained with ASCA which results in a range of predicted PSPC count 
rates ~ 0.2-0.8 count s~^ at this epoch, but which encompasses that observed. 



In § |6.4.3| we suggested that the spectrum in ^ 1 keV was composed of a fraction 
Df ~5% of the underlying continuum which is allowed to escape without suffering any 
attenuation, the emission spectrum of the photoionized ionized gas and thermal emission 
from an optically-thin plasma. At energies 1.5 keV, the observed spectrum becomes 
increasingly dominated by the underlying continuum transmitted by the ionized absorber 
(e.g. Fig. 9). The emission from the optically-thin plasma, identified with the extended 
emission observed in the source, is assumed to be constant with time and under the above 
assumptions is of relatively minor importance to the observed flux (corresponding to a 
count rate in the 0.1-1 keV band of the ROSAT PSPC of 10^^-0.15 count s^^). The 
variability timescales for remaining spectral components differ, with that for the transmitted 
continuum determinded by the the variability timescale of the central source and/or the 
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timescale over which significant changes in the absorbing column density occur within 
the cyhnder-of-sight (above). The timescales on which variations in the intensity of the 
scattered continuum and emission from the photoionized gas occur are determined by the 
hght travel-times to the corresponding regions. If we assume that the rise in the 1-2 keV 
count rate observed by ROSAT is due to a brightening in the central source, then the 
lack of a corresponding increase < 1 keV seen in NGC 4151(4,5) implies these regions lie 
;^ 4 light-days. However if the rise seen by ROSAT is simply the result of a change in the 
absorbing column, then no such constraints can be placed on the location of these regions. 

We note that Weaver et al (1994a) presented the results of a BBXRT observation of 
NGC 4151, and considered a model similar to our model C(ii) (but with Q = Att). They 
found (see Weaver et al §4.3) T - 1.5, Nh,z ~ 10^^ cm'^ Df - 0.04 (from their Fig.l4), and 
Ux ~ 10~^ for their 'high-state' dataset, decreasing by a factor ~ 2 during the 'low-state' 
1 day later. 

NGC 4151 has long been known to exhibit variable absorption features in UV (e.g. 
Penston et al 1981). However the wide range of ionization states observed, including both 
low (Lyman series. Cm) and high (Nv, Ovi) ionization UV lines, exclude a single-zone 
photoionization model accounting for the absorption in both the UV and X-ray bands 
(Kriss et al 1995). 

Mrk 766: 

This Seyfert 1.5 galaxy was first detected in X-rays by Einstein (Kriss, Canizares & Ricker 
1980) and has since been revealed to be one of the most highly variable AGN (e.g. Ghosh &: 
Soundararajaperumal 1992b). In our analysis of the data from ASCA observation perfomed 
in 1993 December, we found Mrk 766 to be second most variable source in our sample, with 
evidence that the amplitude was larger at ^ 2 keV (Paper I). Leighly et al (1996b) have 
presented a more detailed analysis of this dataset and came to the same conclusion, finding 
the amplitude of the variability to be highest at ~ 1 keV. They suggest this is due to the 
intensity of a soft component decreasing slightly whilst the intensity and steepness of the 
underlying continuum ^ 1 keV intensity increased (from P ~ 1.6 to P ~ 2.0). A similar 
change in the form of the continuum with flux state was previously suggested on the basis 
of ROSAT PSPC results for this source (Molendi & Maccacaro 1994; Netzer, Turner & 
George 1994). 

Prom the analysis of the mean spectrum presented here, we find model A(i) to provide 
an unsatisfactory description of the spectra. Pits satisfying our criteria for acceptability 
are obtained if ~50% of the continuum escapes without suffering attenuation by neutral 
material (model A(ii)). However yet superior fits are obtained if the absorbing gas is 
assumed to be ionized (model B(ii)) giving P ~ 2, Ux — 0.1, Nh,z ~ 1-6 x 10^^ cm~^ and 
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Df 0.5. We note that such a model requires no additional soft component (eg. Fig 9). 
Indeed, we find models C(i) h C(ii) to offer only small improvements in the goodness-of-fit 
and thus find no compelling evidence for significant emission from the ionized gas. In no 
case is there a requirement for any addition absorption by neutral gas in excess of Nf^Q. It 
should be noted however, that that models B(i)-C(ii) all have I dof) < 0.05. 

Leighly et al also report the presence of an absorption feature, identified with Ovil, 
and presented the results of photoionization calculations assuming a model similar to our 
model B(i). They find a best-fitting column density Nh,z ~ 6 x 10^^ cm~^, in agreement 
with our findings (Table 5). Unfortunately Leighly et al do not specify the limits on their 
assumed continuum, thus we are unable to compare their values to Ux directly, however 
they do find an increase in ionization parameter with intensity. 

NGC 4593: 

This Seyfert 1.0 galaxy was first detected in X-rays by UHURU (Forman et al 1978). In 
our analysis of the data from ASCA observation performed in 1994 January, we find model 
A(i) provides a satisfactory description of the data analysed, but fails our criteria for 
extrapolation < 0.6 keV. Fits satisfying our criteria for acceptability are obtained if ~70% 
of the continuum escapes without suffering attenuation by neutral material (model A(ii)). 
However yet superior fits are obtained if the absorbing gas is assumed to cover the entire 
source and is ionized (model B(i)), giving F ~ 2, Ux — 0.1, Nh^z ~ 2 x 10^^ cm~^. This 
confirms the findings of R97 who found the addition of Ovil and Ovill edges {toj ~ 0.3 
& To8 ~ 0.1) significantly improved the fit to a single powerlaw model (F ~ 2) to these 
ASCA data. R97 also fitted a single-zone photoionization model and found U§^'^ ~ 0.05, 
Nh,z ~ 3 X 10^1 cm-2, F ~ 1.9. 

We find that statistically there is no strong requirement for any of the underlying 
continuum to escape without suffering attenuation by the ionized material, or for any 
significant emission from the ionized gas, (models B(ii)-C(ii)). In all cases there is no 
required for any addition absorption by neutral gas in excess of N^^q. We note that a 
substantially fiatter spectral index has been suggested in the past (eg. F ~ 1.1, T91). 

MCG-6-30-15: 

This Seyfert 1.0 galaxy was first detected in X-rays by HEAO-1 A-2 (Marshall et al. 1979). 
Absorption due to ionized gas first suggested by Ginga observations, with Nh,z ~ 10^'^ cm~^ 
(Nandra, Pounds & Stewart 1990). Our analysis of the data from two ASCA observations 
separated by ~ 3 weeks in 1993 July, we confirm the presence of ionized gas. However 
whilst none of the models presented in §^ provide a fit which formally satisfies our criteria 
for MGC-6-30-15(l), this is primarily due to the presence of a deficit at ~ 1 keV (Fig.6). 
An acceptable solution is obtained for MGC-6-30-15(2) for model C(i), but again some 
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evidence for a relatively small deficit at ~ 1 keV (Fig. 10). 

However, as further discussed in §[7.3| there is a significant increase in column density 



(from Nh^z ~ 6 X 10^^ cm~^ to ~ 9 x 10^^ cm~^), along with a slight increase in ionization 
parameter (from Ux ~ 0.11 to ~ 0.15) despite a factor ~ 2 decrease in the intensity of the 
illuminating continuum between the 2 observations (e.g. Fig 16). This, the presence of the 
~ 1 keV deficit, and a decrease in the depth of the Ovill edge on a timescale ~ 10^ s (but 
not in the depth of the Ovil edge) observed in an ASCA observation performed in 1994 
July (Otani et al 1996) are impossible to reconcile with that expected from a uniform shell 
of gas reacting to changes in the intensity of the photoionizing continuum. 

A single-zone photoionization models have also been applied to these data by Fabian 
et al (1994) and Reynolds et al (1995), assuming an ionizing continuum consisting of 
a single powerlaw of F = 2. Using Fig lb to convert from their quoted values of ^ we 
obtain results in good agreement with those found here {Nh.z ~ 6 x 10^^ cm~^, Ux ~ 0.09 
and Nh,z ~ 13 x lO^^ cm'^, Ux ~ 0.11 for the MCG-6-30-15(l,2) datasets respectively). 
Photoionization models for the extended (4 day) ASCA observation carried out in 1994 
July have been presented by Otani et al (1996) and R97. The source exhibited a factor 
~ 4 variability in the observed fiux during this observation (covering the two fiux states 
represented by MCG-6-30-15(l,2)) leading Otani et al to divide the observation into 2 
parts. They found z ~ 5 x 10^^ cm~^, Ux ~ 0.04 for the first 300 ks of the observation, 
but Nh.z ~ 10 X 10^^ cm~^, Ux ~ 0.18 for the last 60 ks. For reference, R97 analysed the 
time-averaged spectrum and obtained t/^^^ ~ 0.05 and Nh,z ~ 6 x 10^^ cm~^) The physical 
implications of these results have been further explored by Reynolds & Fabian (1995). 

Curiously (but see Reynolds & Fabian 1995) MCG-6-30-15 appears to be heavily 
reddened in the optical/UV thus little is known regarding the presence of UV absorption 
features. 

IC 4329A: 

This Seyfert 1.0 galaxy was first detected in X-rays by Ariel-V (Elvis et al 1978). In 
our analysis of the data from ASCA observation performed in 1993 August, we find no 
model in which satisfies our formal criteria for acceptability. However we do find a 
significant improvement in the fits for models including ionized gas, with Ux ~ 0.02-0.04, 
Nh,z ~ 2 X 10^^ cm~^ and F ~ 1.9 for all models B(i)-C(ii). A substantial column density 
of neutral gas is also required for this source Nh^q ~ 4 x 10 cm" , far in excess of NfQ, but 
consistent with previous findings and the edge-on orientation of the host galaxy. Cappi et al 
(1996) and R97 have both performed independent analyses of this ASCA dataset. Both find 
the addition of two edges to a single powerlaw model (F ~ 1.9) significantly improved the 
goodness-of-fit over the 0.6-1.0 keV band. R97 fixed the edge energies to those appropriate 
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for Ovil and Ovill and found tq^ ~ 0.6 and tqs ~ 0.1. Cappi et al allowed the line energies 
to be free and obtained a different ratio of optical depths, but otherwise similar results. 
Both Cappi et al and R97 also fitted the data with a similar photoionization model and 
obtained ^ ~ 10 (= f/f^^ ~ 0.02), Nh,z ~ 3 x lO^i cm'^ and Nha ~ 3 x 10^^ cm'^. These 
results are consistent with those found in a ROSAT PSPC observation (Madejski et al 
1995). 

We also found evidence for a strong Compton-reflection component in this source 
with ^ ~ 4 (§ |6.4.1| ), similar to that found by Cappi et al. This hard tail was also found 
in HEAO-1 (W95) and Ginga (Miyoshi et al 1988; Piro, Yamauchi & Matsuoka 1990) 
observations, and offers an explanation the CGRO OSSE spectrum of this source (Madejski 
et al 1995). 

None of the models quoted by Cappi et al (1996) and R97 formally satisfy our criteria 
for acceptability (though some of the fits described by Cappi et al are considered acceptable 
to those workers). It consider it likely that the poverty of the various fits to the ASCA data 
are a combination of unaccounted for systematics errors in the calibration and additional 
spectral complexity in the soft X-ray band. As can be seen from Table 10, even we applied 
the Compton-reflection fit, an excess of counts remains < 0.6 keV. A similar effect was 
reported by Cappi et al. This may be the result of a real steepening of the intrinsic emission 
from IC 4329A, or the result of contamination of the ASCA data by emission from the 
nearby elliptical galaxy, IC 4329 (see Madejski et al 1995). 

NP94 found evidence for absorption by ionized gas in Ginga observations, with 
Nh,z ~ 10^2 cm-2. 

NGC 5548: 

This Seyfert 1.2 galaxy was first detected in X-rays by OSO-7 (Hayes et al 1980). 
Absorption by ionized gas in NGC 5548 was first suggested by Ginga observations (Nandra 
et al 1991) and confirmed by the ROSAT PSPC (Nandra et al 1993). In our analysis of 
the data from ASCA observation perfomed in 1993 July, we find acceptable solutions with 
Ux ~ 0.1, Nh,z ~ 3 X 10^^ cm-2 and P ~ 1.9 for all models B(i) & C(i). Consistent results 
are also obtained for model B(ii), with no requirement for an unattenuated {Df = 0). This 
confirms the results of R97 who found the addition of Ovil and Ovill edges (ror ~ 0.25 & 
To8 ~ 0.16) significantly improved their fit to a single powerlaw model (P ~ 1.9) to these 
ASCA data. R97 also fitted a single-zone photoionization model and found f/^^'' ~ 0.08, 
^^^^,2 ~ 5 X 10^^ cm~^, assuming a single optical-X-ray continuum P ~ 1.9. Mathur, Elvis 
& Wilkes (1995) have also presented photoionization calculations based on previously 
published ASCA results and obtain Ux ~ 0.06 (for their quoted optical-X-ray continuum) 
assuming Nh.z ~ 4 x 10^^ cm~^. 
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We find a significant improvement in the goodness-of-fit is achieved assuming model 
C(ii). However the best-fitting solution occupies the region of high-A^j:^ ^, high-[/x, high-Dj 
parameter-space and does not extrapolate well < 0.6 keV (Fig 12). A yet superior 
improvement was found in § H.4.2| consisting of the combination of a steep powerlaw 



{Tg ~ 3.2) dominating the lowest energies and a fiatter powerlaw {Th ~ 1-5) dominating at 
highest energies. The ionized gas in this case has Ux ~ 0.2, Nh^^ ~ 7 x 10^^ cm~^ but also 
requires a neutral column density {Nfjfi ~ 1-7 x 10^^ cm~^) far in excess of A'^^^q. 

A spectral fiattening towards higher energies also seen in HEAO-1 (W95) and Ginga 
(NP94) along with evidence for absorption by ionized gas, with Nh,z ~ 5 x 10^^ cm^^. A 
'soft excess' has been proposed in NGC 5548 from the analysis of EXOSAT (TP89) and 
ROSAT (Nandra et al 1993; Done et al 1995) datasets, but we find any such component to 
be insignificant > 0.6 keV from the analysis presented here. 

NGC 5548 has been observed to have variable absorption lines of CiV, Nv and Lya in 
the UV (Shull & Sachs 1993; Mathur, Elvis & Wilkes 1995). A detailed comparison between 
the UV and X-ray absorption systems indicates the same gas is probably responsible 
(Mathur, Elvis & Wilkes 1995). 

Mrk 841: 

This Seyfert 1.5 galaxy was first detected in X-rays by Einstein (Tananbaum et al 1986) 
and has since been observed to exhibit dramatic spectral variability (eg. George et al 1993 
and references therein). In our analysis of the data from two ASCA observations performed 
in 1993 August and 1994 February, we find model evidence for ionized gas in Mrk 841(1), 
with model B(i) offering an acceptable description of the data with F ~ 1.9, Ux ~ 0.16, 
Nh,z ~ 3 X 10^^ cm~^. However the features are relatively weak (Fig 6) and not required by 
the Mrk 841(2) dataset, presumably due to its lower signal-to-noise ratio. A superior fit is 
obtained for Mrk 841(2) assuming model B(ii). However the solution occupies the region of 
h.i^-NH,z) high-f/x, high-Dj parameter-space and the parameters are poorly constrained, 
primarily since any derivations from a single powerlaw are rather subtle (e.g. Fig 9). R97 
also found a lack of evidence for Ovil and Ovill absorption edges in their analysis of the 
ASCA data from Mrk 841(2). Interestingly, all models suggest the underlying continuum 
was fiatter during Mrk 841(2) (typically F ~ 1.6-1.7) although the observed fiuxes in the 
2-10 keV band were similar (-F2-10 ~ 1-2 and 1.1 x 10^^^ erg cm~^ s~^ respectively). 

No significant improvement in the goodness-of-fit are obtained for models C(i) h 
C(ii) with the high values for the intensity of the ionized emitter derived in these cases 
(n/47r ~ 50) most likely being an artifact of a separate soft component. The existence of 
a separate soft component has been suggested by some (but not all) previous observations 
(e.g. see George et al 1993). Most recently, the spectrum from an extended ROSAT PSPC 
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observation of Mrk 841, performed over a ~ 6 day period in 1992 can be well described by 
a powerlaw with F ~ 2.2 (i.e. steeper than that found here) along with a kT ~ 0.08 keV 
blackbody (Nandra et al 1995). However the data were equally well described by a yet 
steeper powerlaw (F ~ 2.4) and a photoionized absorber with Nh,z ~ 2 x 10^^ cm~^ 
Ux ~ 0.04 (assuming a powerlaw optical-X-ray continuum). Unfortunately, as described in 



7^, given the signal-to-noise ratio of all the datasets and the dramatic spectral variability 
exhibited by the source, a further interpretation of these results is somewhat problematic. 
However we do note that the extrapolation of the underlying powerlaw implied by the 
PSPC spectrum does agree remarkably well with the UV data, implying the luminosity of 
the component dominating the 13.6-100 eV band may have been grossly overestimated in 
the past (e.g. see fig 3 of Nandra et al 1995). 

NGC 6814: 

In the early 1990s, NGC 6814 had its 15 months of fame due to the apparent detection of a 
periodicity in its X-ray lightcurve (Mittaz & Branduardi-Raymont 1989; Done et al 1992) 
spawning a frenzy of theoretical activity. However, subsequent observations by the ROSAT 
PSPC revealed the periodicity to be the result of strong contamination by a Galactic source 
within the field-of-view of previous instruments (Madejski et al 1993). It now appears that 
the X-ray emission from this Seyfert 1.2 galaxy is relatively weak. ASCA observed NGC 
6814 on 3 occasions prior to 1994 May, but only one dataset (NGC 6814(1)) met the criteria 
for inclusion in our sample (Paper I). The count rate for this observation, performed in 
1993 May, is a factor ten lower than any of the other sources in our sample. Given the 
signal-to-noise ratio, it is perhaps not surprizing that we find the dataset to be consistent 
with a simple powerlaw (F ~ 1.7) absorbed by Nf^^. R97 also found a lack of evidence for 
Ovil and Ovill absorption edges in their analysis of this dataset. 

Mrk 509: 



This Seyfert 1.2 galaxy was first detected in X-rays by Ariel-V (Cooke et al 1978). In 
our analysis of the data from an ASCA observation performed in 1994 April, we find a 
simple powerlaw (model A(i) with F ~ 1.9) provides a reasonable description of the data 
(§ ^.1|) . Whilst a significant improvement in the goodness-of-fit is achieved assuming model 
A(ii), ~ 80% of the continuum is unattenuated giving rise to a relatively subtle change in 
the spectrum ^ 1 keV (Fig 4). A further, significant improvement in the goodness-of-fit 
is obtained for model B(i), with F ~ 2.0, Ux ~ 10~^, Nh,z ~ 8 x 10^° cm~^, but again 
the features imprinted on the continuum are relatively weak (Fig 6). This confirms the 
results of R97 who found the addition of shallow Ovil and Ovill edges (r^y ~ 0.11 & 
To8 ~ 0.04) significantly improved their fit to a single powerlaw model (F ~ 2) to these 
ASCA data. R97 also fitted a single-zone photoionization model and found Ux^'^ ~ 0.1, 
Nh,z ~ 3 X 10^^ cm~^, the discrepancy with our results possibly arising as a result of their 
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flatter continuum (F ~ 1.9). 

No significant improvement in the goodness-of-fit are obtained assuming the subsequent 
models in §P, despite the suggestion of line-like soft-excess by Einstein (T91) and EXOSAT 
(Morini, Lipani, Molteni 1987) observations. However a significant improvement is achieved 
if a moderate (1 ^ JF ^ 5) Compton-reflection component is included (§ |6.4.1| ). Such a 
component is also suggested by HEAO-1 (W95) and Ginga (NP94) observations of the 
source and its inclusion steepens slightly the underlying powerlaw (AF ~ 0.2) and gives rise 
to a larger ionized column density {Nh,z ~ 2 x 10^^ cm~^) whilst leaving Ux as for model 
B(^). 

NGC 7469: 

This Seyfert 1.2 galaxy was first detected in X-rays by UHURU (Forman et al 1978), and 
observed by ASCA on 3 occasions over a period of ~ 9 days in 1993 November. However, 
only one observation (NGC 7469(2)) met the criteria for inclusion in our sample (Paper I). 
We find a simple powerlaw (model A(i) with F ~ 2.0) provides an acceptable description 
of the spectrum. However there is evidence for some spectral curvature, leading to a 
significant improvement in the goodness-of-fit for models A(ii), B(ii) and C(ii) if the bulk 
of continuum is unattenuated {Df ^ 0.6) whilst the remainder is absorbed by a large 
column density {Nh,z ;^fewxlO^^ cm~^), or if there is a strong (JF ~ 7) Compton-reflection 
component (§ |6.4.1|) . 

Guainazzi et al (1994) have presented a combined analysis of the NGC 7469(1,2) 
datasets. They found a ~ 40% decrease in intensity between the two epochs, with some 
evidence for a steeper spectrum in the soft X-ray band (^ 0.8 keV) when the source is 
brighter. No such spectral variability was seen ^ 2 keV, with both epochs consistent with a 
powerlaw with F ~ 1.9. We see no reliable evidence for a soft-excess in our analysis of the 
NGC 7469(2) dataset, however we note that the results of Guainazzi et al (their Fig. 3) 
indicate such a component would be almost indistinguishable at that epoch. Spectral 
complexity in NGC 7469 ^ 1 keV has been suggested by previous Einstein, EXOSAT and 
ROSAT observations although its form remains illusive (e.g. Brandt et al 1993; Leighly et 
al 1996a and references therein). We note some soft flux is provided in our models which 
include emission but generally this is small (except in the case of our model including a 
Compton-reflection component). 

In their analysis of the mean spectra from the NGC 7469(1,2) datasets, Guainazzi et 
al flnd the putative soft-excess to be consistent with a blackbody with kT ~ 0.08 keV. 
Guainazzi et al flnd no evidence for Ovil and/or Ovill absorption edges, whilst R97 flnd 
the addition of an Ovil edge (of optical depth to7 ~ 0.2) signiflcantly improved the flt to a 
single powerlaw model (F ~ 2.1) to the NGC 7469(1) dataset. 
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In § 6.4.1 we found that a strong Compton-reflection component is allowed by the data. 
Such a result is also obtained from previous HEAO-1 (W95) and Ginga (NP94; Leighly et 
al 1996a) observations. 

Finally we note that NGC 7469 has recently been the subject of an intense monitoring 
campaign in the optical, UV and X-ray bands. Most interestingly, preliminary results from 
the RossiXTE observations reveal the X-rays to have exhibited large-amplitude variability 
apparently uncorrelated with the variations seen in the UV and optical (Nandra et al 
1997d). 

MCG-2-58-22: 

This Seyfert 1.2 galaxy was first detected in X-rays by Ariel-V {Cooke et al 1978). In our 
analysis of the data from an ASCA observation perfomed in 1993 May, we find a simple 
powerlaw (model A(i) with F ~ 1.7) provides a reasonable description of the spectrum 
(§ p.l| ). A similar result was obtained by Weaver et al (1995) from an independent analysis of 
these data, who suggested the source to be a 'bare' Seyfert 1 continuum with weak (if any) 
additional absorption or emission features. These workers also noted that the continuum 
was relatively flat compared to the canonical F 1.9-2.0 for most Seyfert 1 galaxies 
popularized by the Ginga results of NP94. However, as shown in we do find significantly 
better fits when the more complex models are assumed, and in some cases obtain a derived 
index F ~ 1.9 (e.g. models B(ii) & G(ii)). Some of these models obviously contain ionized 
gas, but the constraints of Ux and Nh,z are poor. We note that previous Einstein (T91), 
EXOSAT (Ghosh & Soundararajaperumal 1992a) and ROSAT V'&VC (Turner, George & 
Mushotzky 1993) observations of MCG-2-58-22 have shown evidence for additional emission 
(features) <^ 1 keV. 

In § |6.4.1| we found a solution including Compton-refiection, but its strength is poorly 
constrained (JF ~ 2-40), and all but the lower values difficult to reconcile with the observed 
equivalent width of the Fe fiorescent line. However, we note that NP94 found no evidence 
for a Compton-refiection component in the Ginga data (JF ~ 0.1), but did find evidence for 
absorption by ionized gas in Ginga observations, with Nh,z ~ 10^^ cm~^. 
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TABLE 1 

Ttie Sample 



Name 




z 


a 


( 


10" cm-2) 


< /220 > 
(mJy) 




/2keV 








a 


(erg s^^ 


log L2-10 
(erg s-i) 


Mrk 335 





025 





80 





37 


2.35±0 


94 


3.55 ± 





40 


1.16±0 


06 


-0 


04 


1 


23 


44 


2 


43.5 


Fairall 9 





046 





80 





32 


4.64 ±2 


19 


5.19± 


1 


14 


2.37 ±0 


12 






1 


29 


44 


8 


44.3 


3C 120 





033 





76 


1 


23 


0.31±0 


20 


6.40 ± 


3 


44 


5.29 ±0 


26 


+0 


72 





69 


44 


9 


44.4 


NGC 3227 





003 





45 





21 


0.76±0 


08 


0.19± 





02 


2.15±0 


11 


+0 


25 


1 


00 


42 


3 


42.0 


NGC 3516 





009 





75 





29 


1.49 ±0 


55 


1.75 ± 





45 


7.49 ±0 


37 


+0 


13 





90 


43 


8 


43.4 


NGC 3783 





009 





80 





96 


2.35± 1 


00 


3.30 ± 





72 


3.74 ±0 


19 


+0 


10 


1 


10 


43 


7 


43.4 


NGC 4051 





002 





68 





13 


1.46 ±0 


15 


1.88 ± 





19 


2.67 ±0 


13 


+0 


22 


1 


07 


42 


7 


41.6 


NGC 4151 





003 





74 





22 


6.45 ±2 


97 


2.21 ± 





99 


5.33 ±0 


93 


+0 


20 


1 


21 


43 


2 


43.0 


Mrk 766 





012 









18 












2.26 ±0 


11 










43 


7 


43.1 


NGC 4593 





009 





78 





20 


1.20±0 


36 


1.98 ± 





52 


3.93 ±0 


20 


-0 


02 





98 


43 


5 


43.1 


MCG-6-30-15 





008 





58 





41 


0.56±0 


06 


0.29 ± 





03 


4.70 ±0 


82 


-0 


01 





81 


43 


6 


43.1 


IC 4329A 





016 





16 





46 


0.05±0 


02 








7.45 ±0 


37 


+0 


44 





32 


44 


3 


43.9 


NGC 5548 





017 





83 





16 


2.43 ±0 


71 


3.25 ± 





60 


5.14±0 


26 


+0 


09 


1 


05 


44 


2 


43.8 


Mrk 841 





036 









22 


1.14±0 


11 


3.60 ± 





36 


1.21 ±0 


20 


+0 


07 


1 


17 


44 


2 


43.8 


NGC 6814 





006 





95 





98 












0.13±0 


01 










41 


6 


41.3 


Mrk 509 





035 





85 





44 


3.75 ±0 


85 


3.29 ± 





85 


5.22 ±0 


26 


-0 


05 


1 


12 


44 


8 


44.4 


NGC 7469 





017 





58 





48 


2.62±0 


65 


3.15± 





18 


3.27 ±0 


16 


+0 


21 


1 


05 


44 


2 


43.7 


MCG-2-58-22 





047 





67 





35 


1.55± 1 


21 


4.22 ± 


1 


07 


1.36 ±0 


07 


+0 


08 


1 


20 


44 


6 


44.2 



NoTK . — Columns- (.1) Sourer name; (2) redshift from Veron-Cetty \^rron (1993); (3) axial ratio from De Zotti & Gaskell (1985); 
(4) Galarlir TTt column density from the 21 cm measurements (from Murphy et al 1996; Dickey & Lockman 1990: Elvis. Lockman 
t Wilkes 1989: Stark et al 1992); (5) mean flux at 220 nm, < /220 > (from Courvoisier t Paltani 1992); (6) moan ratio of flux at 
125 nm to tbal at 250 nm (from Courvoisier & Paltani 1992); (7) observed flux al 2 VcX, f^u^y from Ihr rrsiills prrsmird in §5; (8) 
radio-to-optiral spectral index a^o (= log(/5GHz//22o)/5.38) using /gcH? from ^"^ron-C^tTy K' A'rroii (199^1): (0) optiral-To-X-ray 
spectral index n^.r (= l'^2;(./'25o//2kev)/2.548) ; (10) luminosity, Lx- of The iniderlyina; eoiiTiinnim over tlie O.I-IO ke\^ hand from 
the results presoiiTed in §5 (using 17^ — 50 km s~^ Mpc~^ and qci = O.o): (I I) corresponding Inminosity. L-^-io- of the nnderlying 
continuum over The 2-10 ke^" hand. All flnxes are measured at frequencies appropriate for the rest-frame of the source. For objects 
with more than one clataset presented here. f'>\^^\-, J-x ^H'^ -^2-10 ^^re exposure-weighted means. It should be remembered that the 
radio, UV and X-ray measurements are non-simultaneous. 



TABLE 2 
Observation Log 



Observation 


Date 


iea^p(SlSO) 
(10 s) 


(count s ) 


.Vp,,„,(SlSO) A'p„(all) 
(excl 5-7 keV) 




-^^2-10 
/in-12 ^ -2 -1\ 
(10 erg cm s ) 


Keierences 


Mrk 335 


1993 Dec 09 


19.3 


0.528 ± 0.005 


9837 


650 


14 


9.5 ± 0.5 


Re97 


Fairall 9 


1993 Nov 21 


22.0 


0.945 ±0.007 


19983 


1146 


26 


19. 0± 1.0 


Ro97 


3C 120 


1994 Feb 17 


47.4 


1.828 ± 0.006 


83155 


1150 


26 


43.8 ± 2.2 


Re97 


NGC 3227 


1993 May 08 


33.4 


0.644 ± 0.004 


20028 


1162 


II 


25.5 ± 1.3 


Pt94, Re97 


NGC 3516 


1994 Apr 02 


29.4 


2.353 ±0.009 


65410 


1161 


26 


73.5 ± 3.7 


Re97 


NGC 3783(1) 


1993 Dec 19 


17.5 


1.032 ± 0.008 


16442 


1065 


14 


•Ifi.'l ± 2.3 


Ge95, Re97, Ge97 


NGC 3783(2) 


1993 Dec 23 


14.9 


1.267 ± 0.010 


17270 


1034 


14 


55.9 ± 2.8 


Ge95, Re97, Ge97 


NGC 4051 


1993 Apr 25 


27.6 


0.939 ±0.006 


25131 


1015 


12 


21.2± 1.1 


Mi94, Re97 


NGC 4151(2) 


1993 Nov 05 


13.4 


2.126±0.013 


23726 


816 


22 


182.5± 9.1 


We94 


NGC 4151(4) 


1993 Dec 05 


10.8 


2.167±0.015 


18626 


1483 


19 


198.7 ± 9.9 


Wa96 


NGC 4151(5) 


1993 Dec 07 


11.7 


2.455 ±0.015 


23096 


1263 


14 


223.4± 11.2 


Wa96 


Mrk 766 


1993 Dec 18 


33.7 


0.828 ±0.005 


36761 


1207 


26 


18.1 ±0.9 


Le96 


NGC 4593 


1994 Jan 09 


21.4 


1.463 ±0.009 


30079 


1337 


16 


31.8± 1.6 


Re97 


MCG-6-30-15(l) 


1993 Jul 09 


29.8 


1.806 ±0.008 


51748 


1110 


14 


43.2 ± 2.2 


Fa94, Re95 


MCG-6-30-15(2) 


1993 Jul 31 


31.4 


1.192±0.007 


35628 


1088 


14 


35.6± 1.8 


Fa94, Re95 


IC 4329A 


1993 Aug 15 


29.2 


2.075 ±0.009 


56677 


1325 


14 


77.3 ± 3.9 


Ca96 


NGC 5548 


1993 Jul 27 


20.4 


1.795 ±0.010 


34919 


1188 


14 


46.1 ±2.3 


Re97 


Mrk 841(1) 


1993 Aug 22 


30.5 


0.440 ±0.001 


12897 


759 


14 


12.2 ± 0.6 




Mrk 841(2) 


1994 Feb 21 


20.9 


0.332 ±0.0-12 


6623 


420 


14 


10.3 ± 0.5 


Re97 


NGC 6814(1) 


1993 May 04 


40.3 


0.028 ±0.001 


1001 


211 


5 


1.2±0.1 


Ro97 


Mrk 509 


1994 Apr 29 


40.1 


2.008 ±0.008 


77504 


1318 


26 


42.9±2.1 


Ro97 


NGC 7469(2) 


1993 Nov 26 


14.5 


1.052 ±0.009 


15202 


1003 


26 


27.5± 1.4 


Gu94 


MCG-2-58-22 


1993 May 25 


13.3 


0.423 ±0.006 


5309 


860 


17 


14.0 ±0.7 


We95, Re97 



Note. — Columns- (1) Source name; (2) observation date (start of exposure): {?,) exposure time in SJSO of selected data: (4) mean SISO count 
rate in the 0.6-10.0 keV band; (5) number of source photons detected by STSO in tlio 0.6-5.0, 7.0-f 0.0 koV liand: (6) numbor of spectral data points 
(all instruments; 0.6-5.0, 7.0-10.0 keV band) included in the spectral analysis; (7) number of good-quality spectral data points (SISO & STSl) below 
0.6 keV; (8) observed flux in the 2-10 keV band (excluding any Fe line emission) with 5% errors; (9) reforoucos (o other papers discussing Llic same 
ASCA dataset: Fag5 - Fabian et al 1994; Geg5 - George, Turner, Netzer 1995; Ge97 - George ct al 1flQ7b: Gu94 - Guainazzi et al 1994; Le96 - 
Leighly et al 1996b; Mi94 - Mihara et al 199-1; Pt94 - Ptak et al i90-1; Re95 - Reynolds et al 1995; Re97 - Reynolds 1997; Wa96 - Warwick et al 
1996; We94 - Weaver et al 1994b; We95 - Weaver et al 1995. All datasets are also discussed in Papers I Kl II. 



TABLE 3 

Model A({]: Results for a neutral absorber (with Df = - see §5.1) 





(10^' cm" 


"") 


r 


(10^' cm-' 


) 




PAi 


AJV0.6 


W7 


Mrk 335 




(P) 


1 O+0.03 
'i--l-O_0.02 


00+° °^ , 

"■""-0.00 ( 




719 


0.98 


8.1 


1.3 


Fairall 9 


30+° "' 


(P) 


l-98l°o:S? 


00+° °* , 

"■""-0.00 ( 




1230 


0.97 


2.6 


0.9 


3C 120 


1 00 + 0.62 


(P) 


2.oolS:S^ 


46+°°', 

"■ "-0.46 ( 




1186 


0.82 


9.6 


0.8 


NGC 3227 


1 1 +0.18 
^•-'-'-1.89 


(P) 


1 cc+o.oa 

-^•'J'J-0.02 


00+^ °' , 

"■""-0.00 ( 




1268 


0.99 


20.3 


2.4 


NGC 3516 


i.6o+';°« 


(P) 


'■^^-0.02 


00+° '^° , 

"■""-0.00 ( 




3214 


1.00 


128.2 


3.1 


NGC 3783(1) 


5.62iSii 




'•^^-0.03 


00+^ " , 

"■""-0.00 ( 




2121 


1.00 


17.9 


15.1 


NGC 3783(2) 


4 04+" 

'-2.66 




, O7+0.03 
^■^ ' -0.03 


00+^ '^^ , 

"■""-0.00 ( 




2114 


1.00 


22.7 


12.6 


NGC 4051 


13+°" 


(P) 




00+° °' , 

"■""-0.00 ( 




1703 


1.00 


60.7 


1.7 


NGC 4151(2) 


"•^-'-0.00 


(P) 




lo qq+0-66 


(p) 


1995 


1.00 


18.7 


8.8 


NGC 4151(4) 


21+^ '"^ 


(P) 




2(1 10+°™ 


(p) 


3235 


1.00 


16.1 


9.5 


NGC 4151(5) 


21 +^ ''' 


(P) 


-'-•"-'^-0,03 


99 70+0.69 


(p) 


3518 


1.00 


23.9 


8.8 


Mrk 766 


,-. -. ^ij-fi 1 n 
^'■-^"-0.00 


(P) 


l-°^-o.oi 


00 . 

"■""-0.00 ( 


1415 


1.00 


16.4 


1.4 


NGC 4593 


0.16^°-^; 


(P) 


i.ao_Q Q2 


49+°-" 
"■^''-0.49 




1402 


0.92 


17.1 


1.3 


MCG-6-30-15(l) 




{P) 
(P) 


i.88lS:g? 


51+° , 

"■"-0.51 ( 




1942 


1.00 


J 36.0 


2.0 


MCG-6-30-15(2) 


1 04+°" 

" -0.63 


1 fi7+°°2 


00+° '^'^ , 
"■""-0.00 ( 




2051 


1.00 


127.3 


2.5 


IC 4329A 


"^^-0.00 


(P) 


1 84+°°' 
'■*'^-o.oi 


4 14+0-11 




1710 


1.00 


17.6 


2.2 


NGC 5548 


32+° " 

"••'^-0.16 


(P) 


1 O1+0.02 

l-''l-o.oi 


00+° " , 

"■""-0.00 ( 




1377 


1.00 


26.8 


1.3 


Mrk 841(1) 


22+° " 

""^-^-O.OO 


(P) 


1 81+°°^ 

'■*''-0.02 


00+° " , 
"■""-0.00 ( 




834 


0.98 


9.6 


1.3 


Mrk 841(2) 


22+° °' 

"•^^-0.00 


(P) 


1.63lS:S^ 


00+° " , 

"■""-0.00 ( 




407 


0.43 


3.3 


1.3 


NGC 6814(1) 


1 00+°™ 

-^■""-0.02 


(P) 








160 


0.01 


0.0 


1.0 


Mrk 509 


42+°-^' 


(P) 


1 09+0.02 

i-y'j-o.oi 


o^25!2:°^ ( 




1335 


0.69 


•'1.5 


I.I 


NGC 7469(2) 


62+°-^'^ 


(P) 


2 (11+° °' 

^•^.^-0.02 


00+° =* , 
"■""-0.00 ( 




998 


0.54 


1.8 


1.1 


MCG-2-58-22 


34+° !^ 

"••'^-O.OO 


(P) 


1 74+0.03 

'■ '^-0.02 


00+° " , 

"■""-0.00 ( 




931 


0.97 


1.1 


1.0 



Note. — Fits undertaken in the 0.6-10.0 keV band, excluding the 5.0-7.0 keV band. Error 
bars on the fit parameters are 68 per cent confidence limits for 3 interesting parameters, (p) 
indicates that the parameter 'pegged' at the specified limit. Columns- (1) Source name; (2) 
the column density of neutral gas at 2 — (constrained to be > Nf^^ given in Table 1); (3) 
photon index; (4) the column density at the redshift of the source; (5) x^"Statistic for the fit 
for the Npts data points given in Table 2; (6) probability that the x^-statistic should be less 
than for the number of degrees of freedom, dof, given by Npts — 8; (7) ratio of the increase 
in the x^-statistic to the increase in the number of data points when the best-fitting model is 
extrapolated below 0.6 keV; (8) mean data/model ratio for data points below 0.6 keV 



TABLE 4 

Model A(i{]: Results for a neutral absorber, with Df free (see §5.2) 



O 

00 



Name 


Nh.o 
(10"' cm- 


J 


r 




I lU Llll J 


log 13/ 


2 

Xau 


Pau 


^Xo a 
AJV0.6 


Rq 


6 




) 


Mrk 335 


0.40^2 


12 

00 


(P) 


2.40ig 


OS 
07 




n ,1 r + 0.050 

U.^10_Q Q49 


636 


0.46 


2.7 


1 


1 


83 


7 


Fairall 9 


0.38+2 


17 

OS 


(P) 


2.02+2 


06 
04 


305.65l™5;f 


-0.09712:2?^ 


1222 


0.96 


2.0 





9 


7 


4 


3C 120 


1.7912 


OS 
09 




2.04+2 


02 
03 


538.961?™:^= 


170 + 0.062 
U.-L / '>_o.095 


1 112 


0..')2 


9.0 





8 


44 





NGC 3227 


0.22+2 


35 

00 


(P) 


1.66+2 


04 
04 






LI 89 


0.78 


3.1 


1 


2 


76 


7 


NGC 3516 


0.29+2 


01 

00 


(P) 


1.99t2 


04 
03 






1962 


1.00 


31.0 


1 





736 





NGC 3783(1) 


0.87+2 


05 
00 


(P) 


1.89^2 


05 
06 


14.7712:21 


n 700+0.037 


1338 


1.00 


6.4 


2 


1 


618 


5 


NGC 3783(2) 


0.87+2 


04 
00 
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1.84t2 


06 


13.98i;;2J 


-o-668l2:2'38 


1429 


1.00 


10.2 


2 


5 


491 


7 


NGC 4051 


0.13+2 


01 

00 


(P) 


ri r 1 +0 


07 
03 


27.60l2;o2 


r\ rinn + 040 

-0.29912:241 


1439 


1.00 


12.6 


1 


2 


184 


7 


NGC 4151(2) 


0.21+2 


12 

00 


(P) 


1.4012 


06 
06 


22.57i;;2^ 


1 i-.ifi-l-0 041 


911 


1.00 


3.9 


8 


8 


961 


8 


NGC 4151(4) 


0.21+2 


22 
00 


(P) 


i.3oi2 


05 


33.2715;=? 


-1.19912:21 


1577 


0.97 


3.8 


1 


6 


1551 





INLtLj 4101(0) 


0.21+2 


29 
00 


(P) 


l-4UIo 


04 
04 


QQ 1 1+1.30 


1 0^9 + 0.032 
1.0U^_Q Q29 


1439 


1.00 


4.4 


1 


4 


1814 


1 


Mrk 766 


0.16+2 


04 

no 


■pj 


2.26^2 


07 
06 


i2.76i;s 


n o-\ 7+0.045 
^■■^^ ' -0,044 


1201 


0.53 


■iA 


1 


9 
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8 


NGC 4593 


0.20+2 


12 
00 


(P) 


2.01+2 


05 
02 


4.4oi;i5 


u.iUi_Q Q32 
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0.77 


3.4 


1 


1 


35 





MCG-6-30-15(l) 


' ' -n no 


(P) 


2.25+2 


05 
05 


9.5912:?^ 


^ O'>o + 0.031 

U.OOO_Q Q3Q 


1494 


1.00 


28.1 


1 


3 


330 


4 


MCG-6-30-15(2) 


0.4 1+2 on 


P) 


2.1512 


05 
06 


12.05l2J^ 


-n 41 t^ + '^36 


1510 


1.00 


28.7 


1 


5 


386 


7 


IC 4329A 


0.45^2™ 


.P) 
(P) 


1.89+2 


02 
03 


4.481212 


_1 -|'>4+0.136 


1658 


1.00 


1.1 


1 





41 


3 


NGC 5548 


0.17+2 


03 
00 


1.9112 


04 
04 


7.09i;?2 


-0.10912:2^1 


1322 


1.00 


15.0 


1 


2 


49 


1 


Mrk 841(1) 


0.22+2 


14 

00 


(P) 


1.9212 


08 
11 


11 4Q+1S.21 


-o.o76i2:2S 


829 


0.98 


7.6 


1 


2 


4 


5 


Mrk 841(2) 


0.22+2 


IS 

00 


(P) 


1.7312 


14 
12 


/M.ld_29 43 (j,) 


U.U/1_Q Q92 
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0.40 


2.2 


1 


9 


4 


1 


NGC 6814(1) 


0.98+J 


14 
00 


(p) 


1.76+2 


46 
16 


in 70+620.02 (p) 
lU.I^8_j„ 45 , ) 


-o.o59i2:°J?S 
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0.01 


0.1 





9 





5 


Mrk 509 


0.44+2 


16 

(v) 


1.9712 


03 
04 


9 0^ + 0.40 


-0.08812:2^1 


1326 


0.63 


2.7 


1 





8 


9 


NGC 7469(2) 


0.80+2 


19 
19 




2.07+2 


04 
05 


5I9.49I237 92 


-0 244+° "° 


979 


0.37 


2.8 


1 


2 


19 


5 


MCG-2-58-22 


0.53+2 


46 

19 (p) 


1.8412 


14 
08 


201.301^^3^ 


-0 iso+° °^^ 


917 


0.94 


1.1 


1 





13 






Note. — As for Table 3, except error bars on the fit parameters are 68 per cent confidence limits for 4 interesting parameters 
(dof = Npts — 9), Df is the fraction of the observed continuum unattenuated by the column Nh^^ , and _F(^) is the _F-statistic 
for model A[ii) versus A[i). 



TABLE 5 



00 

I 









Model B{i): 


Results for an ionized absorber (with 


= - SEE 


Name 


Nh.o 
(10"i cm-") 


r 


(10"i cm-2) 







Mrk 335 


40+° 


Fairall 9 


0.30+2 


3C 120 


1.68^2 


NGC 3227 


44+° 

-0 


NGC 3516 


48+2 


NGC 3783(1) 


1.14« 


NGC 3783(2) 


0.87+2 


NGC 4051 


0.13+2 


NGC 4151(2) 


2.5112 


NGC 4151(4) 


2.89+2 


NGC 4151(5) 


3.47+2 


Mrk 766 


0.18+2 


NGC 4593 


0.5112 


MCG-6-30-15(l) 


0.4112 


MCG-6-30-15(2) 


0.4112 


IC 4329A 


3.63+2 


NGC 5548 


0.1712 


Mrk 841(1) 


0.2212 


Mrk 841(2) 


0.2212 


NGC 6814(1) 


0.9812 


Mrk 509 


0.4212 


NGC 7469(2) 


0.6912 


MCG-2-58-22 


0.3412 



00 (p 

11 

00 {p 



00 (p 
02 

00 {p 



26 
24 
07 

00 (p 
16 

31 {p 

09 

00 ip) 
06 

00 (p 
35 
41 
05 

00 (p 
16 

00 ip) 
12 

00 {p 
97 

U(p 

00 (p 
17 
16 
28 

00 (p) 



2.20+2 


04 

03 


40 51+?^ E^ 


1.9661J 


048 
660 


641 


0.51 


2.9 


1 


1 


78.1 


1 98+° 
r.yo_Q 


03 

02 




9.828I9 


418 (p) 
698 


1228 


0.97 


2.9 


Q 


g 


1.9 


2 00+2 


02 
02 


(.-,+22.09 
^^■^-^-19,69 (p) 


9.997I2 


003 (p) 
^53 


1182 


81 


11 3 


Q 




3 9 


1.6312 


03 
03 


— 1.00 


0.01-1+2 ™ 


1175 


0.68 


1.3 


1 





91.3 


1 92+S 


03 

02 


8 71+S ?? 

' -^-0.65 


0.10612 21, 


1400 


1.00 


2.7 


1 





1494.2 


1.8112 


07 
05 


23.401^5* 


0.14612 


027 
013 


1134 


0.95 


1.1 





8 


914.7 


1.7812 


06 
04 


22,05 2 97 


0.15812 


018 
014 


1087 


0.91 


1.4 





8 


969.8 


2.29+2 


04 
03 


20.92+1^? 


0.68612 


070 
105 


1283 


1.00 


3.8 


1 


1 


329.7 


1.4012 


05 
06 


49.441^;?^ 


0.14412 


008 
006 


1032 


1.00 


1.8 


1 


2 


754.1 


1.38+2 


05 
04 


95.191123 


O.212I2 


007 
006 


1721 


1.00 


2.3 


1 


3 


1297.3 


1.5212 


04 
04 


96.7411=^ 


0.20812 


006 
005 


1648 


1.00 


4.6 


1 


3 


1424.2 


2.0912 


04 
03 


6.05lJ°2 


0.16312 


027 
026 


1104 


0.03 


1.6 


1 





337.8 


2.02+2 


04 
04 


1 qt; + 0.99 


0.089+2 


061 
079 


1344 


0.63 


3.1 


1 


1 


57.4 


2.1412 


03 

02 


5.85l2:S 


O.IIII2 


019 
017 


1 188 


0.97 


7.6 


1 





699.4 


2.0212 


03 
03 


9.87121; 


O.I71I2 


016 
017 


1161 


0.96 


5.5 


1 





827.9 


1.89+2 


02 
03 


2 14+°" 


0.037+2 


049 
023 


1672 


1.00 


5.4 


1 


3 


29.9 


1.9412 


02 
03 


r, ,,+0.77 


0.11512 


027 
028 


1232 


0.86 


5.5 


1 





138.9 


1.9212 


07 
06 




0.16312 


164 
079 


816 


0.95 


3.8 


1 


1 


16.1 


1.7012 


09 
OS 


3.00l?^4%. 


„ „K„+10.988 (p) 
".^•^"-0 204 


402 


0.38 


1.6 


1 


1 


5.1 


1.7112 


24 
11 


2.51+f>2' , 

-2. it' (p) 


0.307_|j 3j,g , , 


160 


0.01 


0.1 





9 


0.3 


1.9712 


03 
03 


0.8212S 


O.OIOI2 


039 
035 


1312 


0.53 


4.4 





9 


23.0 


2.02+2 


04 
04 


109.09«?*„f 


, .n+2.S34 (p) 
■ -'-^^ — 4 242 


990 


0.47 


2.1 


1 


1 


8.1 


1.7412 


06 
03 


127.79l?«1? 


7 >;o>;+^^^*i 

1 .OUO_g 


921 


0.95 


1.1 


1 





9.1 



Note. — As for Tables 3 & 4, except error bars on the fit parameters are 68 per cent confidence limits for 4 interesting 
parameters {dof = Npts — 9), Ux is the ionization parameter of the ionized gas, and is the i^-statistic for model B[i) 

versus A[i). 



TABLE 6 

Model B(i{]: Results for an ionized absorber, with Df free (see §5.4) 



Name 


Nh.o 
(10"i cm-=) 


r 




(10^1 cm-2) 






log 13/ 


Xbu 


Pbu 




-^0.6 






Mrk 335 


40+° 


19 

00 (p) 


2.36^2 


15 

14 


11/^ n'>+46S 47 
140.9315582 


0.613lo42o 


-0 


1 <-i£-i+0 090 


631 


0.41 


4.2 


0.9 


5.1 


10.0 


Fairall 9 


U .OD_Q 


27 

00 (p) 


2.oii2 


04 


4-102 ^1 ln\ 


4-11 241 in\ 
U.UU3_„ (,„4 


-0 


086l2;2S 


1223 


0.96 


2.4 


0.9 


t 


4.7 


3C 120 


1 Qn+0 
1.80^0 


10 
10 


2.04+2 


03 

02 


630.96i;2°8°„<|' 


n ns')+°-^" (''J 

"•"°^-0.081 (p) 


-0 


1 f,o-|-0.007 

-'-'^^-0,029 


1146 


0.56 


8.7 


0.8 


t 


34.7 


NGC 3227 


0.38+1; 


90 

16 (p) 


1.64t2 


04 
04 


5.03lf== 


n n-| + 0.001 


-0 


^-^'^-5.085 (p) 


1175 


0.69 


1.4 


0.9 


13.7 


0.0 


NGC 3516 


0.4812 


14 
13 


1.9212 


03 
03 


o 71 +0.76 
B.a_o Y2 


0.10612:2;^ 


-6 


nn(.+0.000 

"J'J'-'-o.ooo (p) 


1400 


1.00 


2.7 


1.0 


462.6 


0.0 


NGC 3783(1) 


0.87+2 


17 

00 (p) 


1.8812 


OS 
07 


25.421^5* 


0.09612:2?^ 


-1 




1121 


0.92 


1.4 


0.8 


204.3 


12.1 


NGC 3783(2) 


0.87+2 


00 (p) 


1.78i2 


04 


22.651^;?? 


U.100_o 


-6 


1^,^^,-1-0.000 

"^"J'-'-o.ooo (p) 


1087 


0.92 


1.4 


0.9 


322.6 


0.0 


NGC 4051 


0.13+2 


00 (p) 


2.35_Q 


04 
03 


100.531^;^;^^ 


0-79210. 133 


-0 


fi/>n+0 024 
260lo.O36 


1178 


1 .00 


3.6 


0.9 


223.9 


83.2 




0.21+2 


00 (p) 


1.49Io 


08 

09 


43-7413 01 


n ni^o+0.005 


— 1 


Qcc+o.oeo 

'^'J'J-0.019 


868 


0-94 


2.0 


1 .0 


88-1 


128. 1 




0.21+2 


00 (p) 


1 49 + ° 


05 
06 


''6.U'-3.S4 


U.U80_o 012 


1 

— 1 


Qoo+0.048 
'*^*'-0.042 


1 4 Sfi 


n fin 

U.DU 




Z.-l 


u.y 


on n 
yu.u 


on 1 /I 


NGC 4151(5) 


0.21+2 


48 

00 (p) 


1.52t2 


05 


63.50+*°,^ 


0.06712:2;? 


-1 


4-19+0.028 

— 055 


1344 


0.96 


2.9 


1.0 


88.7 


231.5 


Mrk 766 


0.18i2 


00 (p) 
30 (p) 


2.1812 


08 
07 


]6.07lt;r. 


o-ii;il2.S33 


-0 


■lia+o. 044 

^"-0.09G 


1088 


0.01 


2.1 


1.0 


124.4 


17.4 


NGC 4593 


0.50+2 


2.02+2 


04 
04 


1 07+1.17 
l.»'_0.75 


o.o89l2:2fo 


-6 


000+°°°° 

"""-0.000 (p) 


1344 


0.63 


3.1 


1.1 


21.8 


0.0 


MCG-6-30-15(l) 


0.41+2 


00 (p) 


2.1412 


03 
03 


^•^^^ 0.67 


lll + ° "22 

U.44±_Q 


-6 


000+°°°° 

"""-0.000 (p) 


1188 


0.97 


7.6 


1.0 


283.6 


0.0 


MCG-6-30-15(2) 


0.41+2 


00 (p) 


2.0212 


04 
03 


9 87+' °' 
^•°'-1.01 


o-i7il2:2;^ 


-6 


000+°°°° 
"""-0 000 (p) 


1161 


0.96 


5.5 


1.0 


324.3 


0.0 


IC 4329A 


3.57+2 


43 

39 


1.8912 


02 
03 


07+0.55 
^•^ ' -0.S7 


0.032l2i2 


-6 


000+"°"° 

"""-0,000 (p) 


1672 


1.00 


2.8 


1.3 


t 


0.0 


NGC 5548 


0.17^2 


06 

00 (p) 


1.94+2 


03 
03 


01+0.87 
■'■-^ ' -0 8S 


n 1ir, + 0.030 
w. I 1 ■-'_o.032 


-6 


000+°°°° 

"""-0.000 (p) 


1 232 


0.87 


.5.-5 


1.0 


86.1 


0.0 


Mrk 841(1) 


0.22+2 


00 ( p) 


1.9312 


07 
07 




0.168l2;o99 


-0 


479+0.409 
^'^-5.528 (p) 


816 


0.96 


3.6 


1.1 


11.7 


0.0 


Mrk 841(2) 


22+° ''° 


1.7812 


IS 
15 


48 85+'*'-" 


A ot;c, + 1.036 


-0 


-1 99+0.030 
-^■^^-0.087 


396 


0-32 


0-7 


1-0 


6-9 


6-2 


NGC 6814(1) 


n as+if4 

"■""-O.OO (p) 


1.7112 


33 

22 


„ „O + 630.64 (p) 
"••''^-0.00 (p) 


n nni+^-^^^ (p) 

"•""-1-0.000 (p) 


-0 


OC7+0-143 (p) 
°'''-5.143 (p) 


160 


0.01 


0.1 


1.0 


0.0 


0.0 


Mrk 509 


0.42^2 


20 

00 (p) 


1.9712 


03 
03 


0.8212°^ 


010+° °" 

U.U1U_Q Q07 


-2 


4t;t;+2.322 
^■'■'-3.545 (p) 


1312 


0.54 


4.4 


0.9 


14.0 


0.0 


NGC 7469(2) 


0.82+2 21 


2.0812 


06 
05 


630.96l2,2°„tr' 


094+° °" , . 

"-"■"-0.093 (p) 


-0 


-104+0.106 
-lo4-0.065 


979 


0.38 


2.9 


1.2 


0.0 


1 I.I 


MCG-2-58-22 




1.9012 


21 
13 


630.09+;i°8°if' 


1 477+0.462 
'•^"-1.476 (p) 


-0 


-|4q+0.081 
-1^^-0.100 


913 


0.93 


1.3 


0.9 


3.7 


7.4 



Note. — As for Tables 3-5, except error bars on the fit parameters are 68 per cent confidence limits for 5 interesting parameters [dof — Npt^ — 10), -f'(^ff) 
is the i^-statistic for model B{ii) versus A{ii), and F{^i) is that for model B{ii) versus B{i). The increase in x'bh compared to X^u for the datasets marked 
with a f is an artifact of the limited range in Nh,z~Ux parameter-space considered in our lOW models. For these 3 datasets a superior fit is obtained assuming 
model A(ii) in which the absorbing material is assumed to be neutral (i.e. Ux — 0). The grids used for models B(i)-C(ii} have a minimuum Ux — 10~^ 
(§3). Due to the interplay between the data and model during the spectral analysis, the solutions for these 3 datasets most likely represent local minima in 
X^-space. However we note that in all 3 cases yet superior fits are obtained for subsequent models presented here, and none of our conclusions rely upon the 
results from these datasets. 



TABLE 7 

Model C(i): Results for an ionized absorber (with Df = 0) and an ionized emitter (see §5.5) 



Name 


Nn.a 


r 






log(Q/47r) 


Xci 


AJVo.6 


-Ro.6 -P(S-) 




(10"i cm-=) 




(10" cm-") 













Mrk 335 


0.40+° 

— 


Fairall 9 


31+S 


3C 120 


1 69+° 


NGC 3227 


0.22+° 


NGC 3516 


0.50+° 


NGC 3783(1) 


1.62+2 


NGC 3783(2) 


1.16+2 


NGC 4051 


0.13+2 


NGC 4151(2) 


2.97+2 


NGC 4151(4) 


3.15^2 


NGC 4151(5) 


3.62^2 


Mrk 766 


0.20+2 


NGC 4593 


0.51^2 


MCG-6-30-15(l) 


0.41^2 


MCG-6-30-15(2) 


0.41^2 


IC 4329A 


3.76^2 


NGC 5548 


0.17+2 


Mrk 841(1) 


0.22^2 


Mrk 841(2) 


0.22+2 


NGC 6814(1) 


2.23+; 


Mrk 509 


0.42+2 


NGC 7469(2) 


0.72^2 


MCG-2-58-22 


0.42+2 



00 (p) 
19 

01 ip) 



29 ip) 
06 

00 (p) 



31 
31 
27 

31 (p) 

00 Ipj 
34 ■ " 
00 ip) 
54 
53 
12 

00 (p^ 
28 

37 

00 (p) 
77 

25 ip) 
20 

00 (p) 
25 
19 
45 

OS ip) 



2.17ig 


06 


51-50 23 79 


2.291IJ 


559 


-0 


034l2;54s 


633 


0.43 


3.0 


1.1 


8.1 


1 97+2 


04 
03 


Afi 17+156.04 


5.802+5 


444 {p) 

642 





732+°-*^'* , > 

'■'^-4.732 ip) 


1225 


0.97 


2.5 


0.9 


2.8 


1 q8 + ° 


02 
02 


68 c,S+45.24 


9 847+2 


153 ip) 
7SS 


Q 


.,7r,+0,252 


1168 


72 


10 5 


8 


13 7 


1.64^2 


04 


. 90+0.50 


0.00912 


037 





828+2'5?i 


1173 


0.67 


6.0 


0.8 


2.0 


1.92+2 


03 
03 


8.74l2-?i 


0.10712 


017 
013 


-1 


063l2ii? 


1400 


1.00 


2.5 


1.0 


0.0 


1.83+2 


OS 
05 


23.77l*'2? 


0.14912 


027 
OlS 





I80l2;igo 


1104 


0.86 


0.5 


0.9 


28.7 


1.80^2 


07 
07 


no 9t;+2.94 


0.16212 


019 
027 


-0 


11412,211 


1076 


0.87 


0.8 


0.9 


10.5 


2.21_„ 


05 
04 


1 f? nfl+2.54 
iu.ay_3 gg 


u.uyo_ Q 


157 
111 


— u 


117+0.105 

' -0.094 


110-1 
1 lz4 


u.yy 


ft 1 
U. 1 


1 1 
1 . 1 


1 -10 Q 
14z .0 


1.40+2 


07 
06 


48.90+152 




009 
006 





1 10+0. 163 
l-^°-0.230 


992 


1.00 


1.8 


1.2 


33.6 


1.39+2 


05 
05 


94.9711^3 


0.21ltg 


008 
006 


-0 


1170+0.253 
^'^-0.556 


1709 


1.00 


2.4 


1.4 


10.3 


1.53t2 


04 
05 


96.89l|ii 


0.207tg 


007 
005 


-0 


C77+0.336 
'-3.423 Ip) 


1643 


1.00 


4.7 


1.3 


3.8 


2.07t2 


06 
04 


5.36i;io 


0.136tg 


050 
044 





105l2i?5 


1085 


0.01 


2.5 


0.9 


20.9 


2.02+2 


04 
04 


1.95lJi^ 


O.OQllg 


067 
0S4 


-4 


OOO^^.ooo 


1344 


0.63 


3.2 


1.1 


0.0 


2.1412 


03 
03 


5.83+2-:^ 


O.lOQlg 


025 
022 


-0 


988^3,0'' 


1 188 


0.97 


7.4 


1.0 


0.0 


2.0012 


06 
04 


8.89i;:2^ 


0.144lg 


040 
023 





nqq+9266 

uyy_o 289 


1135 


0.89 


5.6 


0.9 


24.7 


1.8812 


02 
03 


1 69+°" 


0.020^2 


042 
Oil 





OC2+0-434 
°''^-4.852 (p) 


1666 


1.00 


5.0 


1.3 


4.7 


1.9312 


03 
03 


3.0912:2? 


0.104+S 


034 

046 


-0 


ino+0.455 
-'"'^-3.898 Ip) 


1228 


0.85 


5.5 


0.9 


3.8 


1.88+2 


06 
07 


n qq+2.46 

0-681 J;|2 (p) 


0.023lg 


056 
020 


1 


OQQ+0.679 
^'^'J-1.183 


810 


0.94 


1.7 


1.0 


5.6 


1.6512 


OS 
07 


0.044lg 


OSl 
037 


1 


349l2;fJ 


390 


0.25 


1.4 


0.9 


12.7 


1.7812 


29 
16 


„ „„ + 624.71 (p) 
"•■'^-0.00 (p) 




945 ip) 
054 (p) 


2 


910l°o15?(p) 


154 


0.01 


0.1 


0.9 


7.3 


1.9712 


03 
01 


0.82121^ 


0.010+^ 


057 

009 (p) 


-4 


000+'^°'^^ 

"""-0.000 (p) 


1312 


0.54 


4.4 


0.9 


0.0 


2.0212 


04 
05 


129.69i;™;^2 


7.952^4 


048 Ip) 

732 


-0 




988 


0.46 


2.3 


1.1 


2.1 


1.7312 


09 
07 


105.56llf,f 


5.732^3 


268 {p) 
7SS 


-0 


08412"9;? ip) 


920 


0.95 


1.1 


1.0 


1.0 



Note. — As for Tables 3-6, except error bars on the fit parameters are 68 per cent confidence limits for 5 interesting parameters {dof — 
Npts — 10), Q, is solid angle subtended by the emitting gas at the central source, and F[^) is the F-statistic for model C{i) versus B{i). 



TABLE 8 

Model C(ii): Results for an ionized absorber and an ionized emitter (see §5.3) 



Name 




Nh,o 


r 






Ux 


logD 


/ 




log(Q/47r) 




Pea 




.R0.6 








(10 = 1 g„j-2) 






(1021 ^j^-2) 
























Mrk 335 





40+° 


62 

00 (p) 


ri J 1 +0 

2.41Io 


37 

42 


7A 07+556.69 (p) 
'^■^'-42,94 


U./!I(_u 216 (p) 


-0 104+°"* 

U.1U4_Q IgQ 


(P) 


1 740+'^-^2® 
1. ^4^ 5 ^42 


629 


0.40 


5.6 


0.9 


2.0 


3.4 


Fairall 9 


Q 


39+° 


24 

09 (p) 


2.02tg 


07 




4-10 172 ln\ 
1 048+ 

.l-"^o_l,047 (p) 




KP) 


(«) 


1221 


0.96 


3 7 


0.9 


1.9 


3.7 


3C 120 


^ 


7Q+0 


12 
10 


2.04+g 


04 
03 


630.96l^^fr! 


n 1^0+1.002 
U.l i„ jjjj 


-0.057lo 


n?T 


P) 


n 777+" ^°* 

U. 1 1 1 _4 777 (p) 


1139 


0.51 


9.6 


0.8 


6.5 


28.9 


NGC 3227 





22+1 


09 

00 (rt 


1.651S 


05 
05 


5.00l?;°r8 


no9+° °^^ 
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Note. — As for Tables 3-6, except error bars on the fit parameters are 68 per cent confidence limits for 6 interesting parameters, (dof — Npt^ — 11) , 0, is solid angle subtended 
by the emitting gas at the central source, is the i^-statistic for model C{ii) versus B{ii), and F[^j) is that for model C{ii) versus C{i). (u) indicates that meaningful 

constraints on Q, could not be obtained (0 < Q/Att < 10). 



TABLE 9 
Additional edge fits (see §6.2) 



Name JVs,o T JV^,, Ux -BojVa rojva XojVe PoNe -P(Tr) 

(10=1 rm-=) (10" cm-") {eV) 

NGC 3516 0'29+o:oo l-86lg;2^ 4.97l2;^^ OMSt'omo 853^11 0.32l° J? 1283 1.00 2.2 1.0 105.0 

NGC 3783(1) T-O^mi (p) ^-'StS.Se 16-43t|:;^ 0.119t°;2^? 953t^^ 0.50^° ;^ 1065 0.60 0.6 0.9 68.4 

NGC 3783(2) (l.87t°a^^l l-1'i-t°°ol 14.74^1;^^ 0.129t|5;2" 888t|^ 0.57t|5;^? 1025 0.52 0.8 1.0 61.9 

MCG-6-30-15(l) 0.41+°°*p) 2.13tg°i^ 4.43^J;™ 0.086^|5;2^* S62tlf O.Uto°l 1166 0.92 7.1 1.0 20.8 

MCG-6-30-15(2) 0.41+° °^, , 2.02tg °* 7.90^5;*^ 0-14lt|!;226 947^^39 O.Utoil 1138 0.90 5.4 1.0 21.8 



Note. — As for Table 5, except error bars on the fit parameters are 68 per cent confidence limits for 6 interesting parameters [dof = Npts — 11), 
"^ONe is the optical depth of the additional absorption edge required at a (rest-frame) energy of EoNey and F[ ^^J^ ) is the i^-statistic for XoNe 
and xii. 



TABLE 10 

POWERLAW AND PARTIALLY-COVERED IONIZED ABSORBER + EMITTER FITS WITH REFLECTION (SEE §6.4.1) 



Name Nh.o T N'„ ^ Ux log log(Q/4,r) JF xl Per Ra.6 -F(^) 

(10=1 ^^-2) (JO" cm-=) 
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+ 0.41 
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+ 0,21 
-0,16 



2.80 
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+ 0.53 
0.78 
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-0.255+° 



.082 
001 

.278 (p) 
088 
.323 
070 
.060 
541 
.107 
336 

.267 (p) 
573 

.071 (p) 
909 (p) 



5,745 (p) 



-0.00512 
0.10312 
0.22112 
-4.000l^ 
-4.000l^ 
-0.585l^ 
-0.993l^ 
0.088i; 
0.716l^ 
0.2071° 
0.181+^ 



485 
986 
307 
546 
249 
151 
068 

000 (p) 
154 



007 (p) 
363 
066 
203 

716 (p) 
015 



4.181 (p) 





619 


0.30 


2.5 


I.I 


10.6 


2 4+^-^ 


1069 


0.85 


1.2 


I.I 


17.22 




1077 


0.95 


1.6 


1.0 


23.3 


4.5+?-^ 


1299 


0.31 


4.4 


1.1 


44.9 


2.61] S 


1176 


0.95 


7.6 


1.0 


11.2 




1123 


0.84 


7.7 


0.9 


11.5 




1568 


1.00 


4.0 


1.2 


166.8 


4.9+?« 


1159 


0.37 


8.2 


1.0 


47.7 




1285 


0.34 


11.2 


0.9 


12.2 


6.9l2j 


966 


0.29 


2.6 


1.1 






904 


0.91 


1.8 


0.9 


9.4 



Note. — As for Tables 3-8, except error bars on the fit parameters are 68 per cent confidence limits for 7 interesting parameters {dof = Npts — 12, J- is the 
enhancement of the reflected continuum relative to that expected from a semi-infinite plane illuminated by a point source viewed face-on, and F[-^) is the i^-statistic 
for xlr ^iid Xcii ■ 



TABLE 11 

Results for a double-powerlaw plus ionized absorber and an ionized emitter (see §6.4.2) 



Name Nh.o T. R,,^ Nh,. Ux log(fi/47r) xlp P'lp ^ Ra.6 F(^) 

(10" cm-2) (10" cm-2) 

Mrk 335 OAOtlil,,, 2.61+?-'^' 2.9 1.57^5;- {^j 16.03^-;?? 0.218l-« 

NGC 4051 0.74lJi? 3Alt\fs^'^ 1.5 1.8515=^,^) ^-^Ottfs 0.31212;^ 0.3421^;^?^ 

IC 4329A 5.29+1-;,-^ 3.22+° «^ 2.9 1.55+2"^^ 7.16+^79 0.161+ro76 -0.290+° ?:;j , , 

NGC 5548 1.74i;:;-^ 3.18tlil 3.1 IM^ait O.WStiHl O.IOO^SS 
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0.29 


2.1 


1.0 


11.3 
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0.95 


9.0 


0.9 


19.5 


1573 


1.00 
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1.2 


77.7 


1163 


0.40 


13.5 


1.0 


43.5 



NoTP. — As for Tables ;l-8. rxrrpi error hars on ihr fil parameters are 68 per rnii roiifirlmrr limits for 7 interesting parameters [dof = .Vpf., — 12), f., and T 
are the photon indices of the two powerlaws (constrained such that 1.5 > > 5.0 and 0.0 > T/i > 2.5), Rs/h is the ratio of the normalization of the soft powerlaw 
divided by that of the hard powerlaw at 1 keV (in the rest-frame of the source), -F-statistic for model this double-powerlaw model compared to 

model C{ii) . 



TABLE 12 
Summary of Fits 









Acceptability of Models 






Ionized 


Notes & Comments 


Name 


A{i) 


A{ii) 


B{{) 


B{ii) 


C{i) 


G{ii) 


gas ? 






















Mrk 335 


no 


|YESt| 


YES 


YESt 


YESt 


YES1 




complex continuum (§6.4.2) 


Fairall 9 


yes" 


1 yes" 1 


yes' 


yes' 


yes' 


yes' 




see Appendix 


3C 120 


YES 


YES 


YES 


YES 


YES 


YES 






NGC 3227 


no 


YES 


1 YES 


YES 


YES 


YES 


V 




NGC 3516 


no 


no 




yes^''^ 


yes"-' 


yes'^'' 


V 




NGC 3783(1) 


no 


no 


YES 


YES 


YES 


YES 




see also §7.1 


NGC 3783(2) 


no 


no 


YES 


YES 


YES 


YES 


\/ 


see also §7.1 




no 


no 






no 


[J23 


V 


spectral variability 




no 


no 


no 


1 YES 1 


no 


YES 




see also g r .2 


IVnP 41 fil (A] 








YES 




YES 


1 V 


see also §7.2 




no 


no 




1 yes'^''^ 1 


no 


yes"'" 


V 


see also §7.2 


Mrk 766 


no 


YES 


YESt 


YESt 


YESt 


YES: 






NGC 4593 




YES 


1 YES 


YES 


YES 


YES 






MCG-6-30-15(l) 


no 


no 


1 yei' 1 


yes' 


yc.f 


yes' 


V 


see also §7.3 


MCG-6-30-15(2) 


no 


no 


yes'^ 


yes' 


YES 1 


YES 




see also §7.3 


IC 4329A 


no 




no 


no 


no 


no 


V 


complex continuum (§6.4) 


NGC 5548 


no 


no 


YES 


YES 


YESt 


yes'' 






Mrk 84 1(1) 


no 


no 


1 YES 


yes'-'' 


YES 


YES 




see also §7.4 


Mrk 84 1(2) 


ft 

yes 


YESt 


YESt 


YESt 


YESt 


YESt V 


see also §7.4 


NGC 6814(1) 


YES| 


YESt 


YESt 


YESt 


YESt 


YESt 




Mrk 509 


YES 


YES 


1 YES 


YES 


YES 


YES 


V 




NGC 7469(2) 


YES 


1 YESf 1 


YESt 


YESt 


YESt 


YES 


V 




MCG-2-58-22 


yes"" 


YES 


YES 


YES 


YES 


YES 






TOTAL 


4 + 4 


10 + 1 


14+1 


15 + 6 


15 + 3 


16 + 5 18 





Note. — DatasoTs marked hy YES satisfy al] our criteria for arreptah]ity (sec §4.1] for tlie given model, and tlie hest-fittina; values arc used in the construction of the 
hatched hislo2;rams in Ti.ajs "2, ?>. 5. 7 S.: II. Dalasels donode hy yes do uol satisfy all the criteria dne lo: " llie presence of (random) noise in the data/model residuals beyond 
that expected from pnrely statistical considerations; failnre To extrapolate < 0.6 keV; the presence of a deficit in connts at ^1 keV (see §6.2). Nevertheless, in many cases 
the model is possih]y stih a resonah]o description of tlie data, and tlie effect of including the best-fitTinc;; valncs shown as the open histograms in Tigs "2. ?>, 5. 7 II. The 
given model is not considered to offer an acceptable description of the datasets marked 'no'. For each dataset, the boxed model indicates which of the models considered in 
§5 is not improved at > 999? confidence by the inclusion of additional parameters. Datasets flagged with f have P{x^ \ dof) < 0.5 (i.e. reduced values less than unity) and 
those flagged with | have P{x^ \ dof) < 0.05. Further notes on individual objects are given in the Appendix. 
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Figure Captions 

Figure 1: (a) The solid curves give the ratio of the 'X-ray ionization parameter' Ux 
(eqn. |I]) to the more traditional ionization parameter U (defined over the entire Lyman 
continuum) for the form of the optical-X-ray continua used here (see §§) assuming aox = 1-0 
(upper curve) and aox = 1-5 (lower curve). The dashed curves give Ux/U assuming the 
continuum is a single powerlaw of photon index F from 13.6 eV to 13.6 keV (upper curve) 
and to infinite energies (lower curve), (b) The solid curves give the ratio of Ux to C, (in units 
of erg cm s~^), defined as ^ = L/uhR"^, where nn and R are as in eqn.[^, but where L is the 
integrated luminosity over the 13.6 eV to 13.6 keV band for the form of the optical-X-ray 
continua used here, assuming aox = 1-0 (upper curve) and aox = 1-5 (lower curve). The 
dashed curve gives Ux/C, assuming the continuum is a single powerlaw of photon index F 
from 13.6 eV to 13.6 keV (as used by Reynolds 1997 amongst others). In both panels, the 
open and filled symbols indicate the corresponding ratios for the optical-X-ray continua 
given in Mathews & Ferland (1987) and Krolik & Kriss (1995) respectively. 

Figure 2: Histograms showing the range of photon index (F) of the underlying 
continua derived assuming the various models described in §|[ The histograms are generated 
by assigning a rectangle of equal area to each dataset, of width determined from the la 
error bars for that parameter as listed in the relevant table. In all cases the hatched 
region shows the histogram which formally satisfy all the criteria given in §^^. The open 



histogram shows the effect of also including the datasets which just fail to meet our formal 
criteria but for which the model is considered likely to be applicable (as described in the 
text and summarized in Table 12). 

Figure 3: As for Fig. 2, but for the column density {Nh,z) of the material at the 
redshift of the source. It should be remembered that in models A(i) & A(ii), the material 
is assumed to be neutral, whilst in models B(i)-C(ii) it is photoionized (see §^). 

Figure 4: The derived spectrum and the corresponding rebinned, mean data/model 
ratios for the 8 datasets for which model A(u) (§ |5.2| ; Table 4) is acceptable and offers an 
improvement at >99% confidence over model A(i). Also shown are the corresponding plots 
for Fairall 9 and IC 4329A as, although the fits to these 2 datasets fail to satisfy our formal 
criteria for acceptability, in neither case is the goodness-of-fit significantly improved by 
the inclusion of the additional parameters in subsequent models presented in §^ (see text). 
The objects are presented in order of increasing Right Ascension. The observed spectrum, 
corrected for the effects of Galactic absorption, is shown in the upper panel of each plot as 
the solid, bold curve. The individual spectral components which contribute to the observed 
spectrum are also shown, as is the powerlaw (dashed) that would be visible if there was 
no attenuation along the line-of-sight at the redshift of the source (i.e. if Nh,z = 0). As 
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described in §0, for each observation the spectral analysis was performed by simultaneously 
fitting the mean spectra from all 4 instruments over the 0.6-5.0, 7.0-10.0 keV band. The 
data/model ratios plotted as filled triangles in the lower panel are the (error-weighted) 
means of the ratios from the individual detectors within this energy range, rebinned in 
energy-space for clarity. The open triangles show the corresponding rebinned, mean ratios 
when the best-fitting model is extrapolated < 0.6 keV and into the 5-7 keV band. It should 
be noted that significantly superior fits are obtained for more complex models in the case 
of Mrk 355 (§|X|), NGC 3227 & NGC 4593 (§pD, Mrk 766 (§0) and Mrk 509 



Figure 5: As for Fig. 2, but for the fraction of the underlying continuum (-D/) which 
does not suffer attenuation by material at the redshift of the source. 

Figure 6: As for Fig 4, but for the 8 datasets for which model B(i) (§ |5.3| ; Table 5) is 
acceptable and offers an improvement at >99% confidence over both models A(i) and A(ii). 
Also shown are the corresponding plots for NGC 3516 and MCG-6-30-15(l) as, although the 
fits to these 2 datasets fail to satisfy our formal criteria for acceptability, in neither case is 
the goodness-of-fit significantly improved by the inclusion of the additional parameters in 
subsequent models presented in §^ (see text). It should be noted that significantly superior 
fits are obtained for more complex models in the case of NGC 3516 ( §6.2| ), NGC 3783(1,2) 
(§U & §13), Mrk 766 & Mrk 841(2) (§U) and Mrk 509 (§p|). 



Figure 7: As for Fig. 2, but for the ionization parameter {Ux) of the photoionized 
material, as defined by eqn. [|. 

Figure 8: The best-fitting values for Nh,z and Ux for each of the models including 
ionized gas considered in Shown are only those datasets for which the model offers 
an adequate description of the data and leads to an improvement at > 99% over the 
corresponding model with neutral gas (i.e. model A(i) or A(iiJ). The solid symbols 
represent the datasets for which the model satisfies all our criteria for acceptability. The 
open symbols represent the datasets which do not satisfy our formal criteria, but for which 
we consider the measurements of N^^^ and Ux to be reasonable estimates. Overlaid on 
each panel are curves showing the r = 0.1 (solid) and r = 0.3 (dotted) for Ovil and 
Ovill. As futher discussed in § |8.5| , ionized gas lying in the region of Ux-Nh,z plane where 
Ux > W~'^'^Nh,z has an optical depth r ^ 0.1 for all ions with edges within the ASCA 
bandpass. For the signal-to-noise ratio typical of ASCA observations the ionized gas will 
thus be almost transparent. 

Figure 9: As for Fig 4, but for the 5 datasets for which model B(ii) (§ |5.4| ; Table 6) 
is acceptable and offers an improvement at >99% confidence over both models A(ii) and 
B(i). Also shown are the corresponding plots for NGC 4151(5) as, although the fit to 
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this dataset fails to satisfy our formal criteria for acceptability, the goodness-of-fit is not 
significantly improved by the inclusion of the additional parameters in subsequent models 
presented in §^ (see text). It should be noted that significantly superior fits are obtained 
for more complex models in the case of NGC 3783(1) (§^ & §§3), NGC 4151(4) {^^) 
and Mrk 509 (Q). 

Figure 10: As for Fig 4, but for the 4 datasets for which model C(i) (§ |5.5| ; Table 7) is 
acceptable and offers an improvement at >99% confidence over both models B(i) and B(ii). 

Figure 11: As for Fig. 2, but for the solid angle (fi) subtended by the ionized gas at 
the origin of the illuminating continuum. 

Figure 12: As for Fig 4, but for the 4 datasets for which model C(ii) ( §5.6| ; Table 8) 
is acceptable and offers an improvement at >99% confidence over both models B(ii) and 
C(i). It should be noted that NGC 4051 and NGC 5548 fail to satisfy all our formal criteria 
for acceptability assuming model C(ii). Furthermore it should be noted that in the case of 
NGC 5548 the model predicts a stronger Fe Ka emission than observed. Caution is urged 
in when interpreting these results for NGC 4051, NGC 5548 and Mrk 509 as the best-fitting 
solutions are really only indicative of relatively subtle curvature in the observed spectrum. 
Such curvature can arise as a result of curvature in the underlying continuum, spectral 
variability during the observation, and/or additional spectral components (see text and 

Figure 13: The mean data/model ratios (in the observer's frame) in the SIS and 
CIS for the 16 datasets satisfying our criteria for acceptability assuming model C(ii). This 
plot was constructed by taking the error- weighted average of the individual data/model 
ratios for each detector from each of the appropriate observations. These averages were 
then rebinned for display purposes. Caution is urged when intrepreting such plots. They 
potentally contain artificial features introduced by time-dependent errors in the calibration 
of the instruments (such as slight offsets in the gain of the detectors at the time of each 
observation), along with real features due to errors in the calibration of the instrument, as 
well as of astrophysical origin. Furthermore such plots are dominated by the datasets with 
the highest signal-to-noise ratio. 

Figure 14: The mean data/model ratios in the 3-10 keV band (source frame) obtained 
in Paper II, and for models B(ii), C(i) and C(ii) presented here. The corresponding profile 
assuming model B(i) is almost identical to that assuming model B(ii) and hence is not 
shown. The ?/-axis for each panel covers the same dynamic range. As described in the text, 
the 5-7 keV band is excluded from the spectral analysis presented here. The profiles were 
constructed from those datasets which satisfy our criteria for acceptability (§[4.1|) assuming 
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the best-fitting model and using the SIS data (only). The vertical, dotted line shows the 
location of the Ka flourescent line from neutral Fe (6.4 keV). There are naturally slight 
differences between the profiles, in particular in the case of models C(i) Sz C(ii) (which 
include Fe line-emission for the photoionized-emitter) . However it can be seen that the 
profiles obtained from the analysis presented here have the same characteristics of the profile 
found in Paper-II Specifically, the peak of the line is consistent with that expected from 
Fel, the wings of the profile show a significant breath (the FWHM of the instrument at al 
resolution is ~ 0.12 keV at these energies), the profile appears assymmetric with more flux 
in the low-energy wing. The implications are of such a profile in the context of models in 
which it is due to emission from the innermost regions of an accretion disk are discussed in 
detail in Paper-II. 

Figure 15: As for Fig 4, but for the 4 datasets for which the double-powerlaw model 
(§ |6.4.2 ; Table 11) provides a significant improvement over the fit obtained assuming model 



C(ii) and in which both powerlaws of similar importance over a sizable fraction of the 
ASCA bandpass. As described in the text, the goodness-of-fit of 6 other datasets are also 
significantly improved assuming such a model. However in these cases the second powerlaw 
component makes a significant contribution only at either the very lowest or very highest 
energies. It should be noted that best-fitting solution to IC 4329A fails to satisfy our formal 
criteria for acceptability. 

Figure 16: The solid contours show the 90% confidence limits (for 6 interesting 
parameters) in the Nh,z-Ux plane for the individual observations of NGC 3783, NGC 4151 
and MCG-6-30-15 assuming model C(ii). The dashed contours show the corresponding 
limits for NGC 4151 when the optically thin plasma with kT ~ 0.7 keV is also included 
in the spectral model (§ |6.4.3| ). The dotted curves show the optical depths of the Ovil and 
Ovill edges in this region of parameter-space. The implications of the changes in the state 
of the ionized gas in each source, particularly with regard to the corresponding behaviour 
of the illuminating continuum, is discused in §|^. 

Figure 17: The derived continuum luminosity, Lx-, in the 0.1-10 keV band (after 
correcting for absorption) versus the parameters of the ionized gas assuming model C(ii). 
The solid symbols represent the datasets for which the model satisfies all our criteria for 
acceptability. The open symbols represent the datasets which do not satisfy our formal 
criteria, but for which we consider the measurements of N^^z and Ux to be reasonable 
estimates. As can be seen from the upper panel, there is no clear relation between Nh,z 
and Lx- However from the lower panel it can be seen that there is some indication that 
Ux oc Lx, which would have implications for the distribution of gas in these sources. 
However it should be noted that the majority of the sources considered here in which Ux 
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is well constrained cover only a decade in luminosity (10^^-10^'* erg s~^). The addition of 
the two quasars PG 1114+445 and MR 2251-178 to the plot argues against any obvious 
relation between Ux and Lx (see §^.51). 

Figure 18: The column density, Nh,z, derived assuming model C(ii) versus the mean 
ratio of the flux observed at 125 nm to that at 220 nm, {fun/ f 220) obs (from Table 1). 
It can be seen that the majority of sources are consistent with 2 ^ {fun/ f 220) obs ^ 5, 
whilst 3 datasets (NGC 3227, MCG-6-30-15(l,2)) have values a factor ~ 10 lower. Clearly 
{fi25/ f 220) obs = {fu5/f22o)int c"'^^, whcrc {fub/f22o)int IS the intrinsic flux ratio, and Ar 
(= ^125 — T220) the difference in optical depths at 125 nm and 220 nm. We stress that 
{ f 12b / f 220) int is not necessarily a good indication of the underlying UV continuum (see 
text). The curves show {f ub/ f 220) obs assuming {fub/ f 220) int = 3.1 and (from top to bottom) 
Ar = 0.0, 0.005 A/'2i, 0.05A/'2i (dotted) and Ar = 0.5Ar2i (dashed), where N21 = A^h.^/IO^^. 
It can be seen that NGC 3227 and MCG-6-30-15(l,2) are consistent with { f ub / f 220) int 
similar to the other sources if their fluxes at 125 nm and 220 nm are absorbed by a column 
density similar to that observed in the X-ray band and Ar = 0.5A^2i- A comparison of such 
a value of Ar with that expected from 'standard' dust models is provided in the text. 
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